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Abstract 

Q_ij We review the phenomenology of models with flat, compactified extra dimensions where all of 

the Standard Model fields are allowed to propagate in the bulk, known as Universal Extra Di- 
mensions (UED). UED make for an interesting TeV-scale physics scenario, featuring a tower of 
Kaluza-Klein (KK) states approximately degenerate in mass at the scale set by the inverse size 
of the compactification radius. KK parity, the four-dimensional remnant of momentum conser- 
^ vation in the extra dimensions, implies two basic consequences: (1) contributions to Standard 

Model observables arise only at loop level, and KK states can only be pair-produced at collid- 
ers, and (2) the lightest KK particle (LKP) is stable, providing a suitable particle dark matter 
candidate. After a theoretical overview on extra dimensional models, and on UED in particu- 
lar, we introduce the model particle spectrum and the constraints from precision electroweak 
tests and current colliders data. We then give a detailed overview of the LKP dark matter 
phenomenology, including the LKP relic abundance, and direct and indirect searches. We then 
discuss the physics of UED at colliders, with particular emphasis on the signatures predicted 
for the Large Hadron Collider and at a future Linear Collider, as well as on the problem of 
discriminating between UED and other TeV-scale new physics scenarios, particularly super- 
symmetry. We propose a set of reference benchmark models, representative of different viable 
UED realizations. Finally, we collect in the Appendix all the relevant UED Feynman rules, the 
scattering cross sections for annihilation and coannihilation processes in the early universe and 
the production cross section for strongly interacting KK states at hadron colliders. 
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1 Introduction 

In 1919, Theodor Kaluza extended general relativity to five dimensional space-time in an effort 
to unify gravity with the force of electromagnetism [1] . Seven years later, Oskar Klein proposed 
that this fifth dimension could be hidden from our observations by being compactified around a 
circle of very small radius [2j. Particles moving through the extra dimension would appear very 
massive, as a result of their extra-dimensional momentum being perceived as a contribution to 
its four-dimensional rest mass. These excited states of ordinary particles have become to be 
known as Kaluza-Klein modes. 

The original proposal of Kaluza and Klein suffered from a number of problems and attracted 
only marginal interest for many decades. In the 1980s, however, the relevance of the work of 
Kaluza and Klein grew dramatically, bolstered by the conclusion that string theory could not 
be formulated in a four dimensional space-time. 

More recently, interest in theories including extra spatial dimensions has undergone a resur- 
gence. In particular, a great deal of attention has been given to two classes of extra dimensional 
theories over the past decade. First, scenarios featuring one or more large (millimeter-scale), flat 
extra dimensions were proposed by Arkani-Hamad, Dimopoulos and Dvali (ADD) in 1998 OS]. 
In the ADD model, all of the fields of the Standard Model are confined to a four dimensional 
brane, leaving only gravity to propagate in the space of the extra dimensions (the bulk). As a 
result, the fundamental Planck scale as perceived in the bulk could be similar in magnitude to 
the electroweak scale, providing a potential connection between those scales and the hierarchy 
problem0 

In addition to the possibility of large, flat extra dimensions, Randall and Sundrum proposed 
that the fundamental Planck scale could be lowered to around the electroweak scale by intro- 
ducing an additional small dimension of space with a large degree of spatial curvature [8]. The 
presence of an exponential "warp factor" in the metric of the Randall-Sundrum model leads to 
the fundamental Planck scale (which is near the electroweak scale) to appear to observers on 
our brane to be much higher, thus addressing the hierarchy problem. 

A further class of extra-dimensional models, which goes by the name of Universal Extra 
Dimensions (UED), was proposed in 2001 by Appelquist, Cheng and Dobrescu [SJ, and is the 
focus of this review. By universal, it is meant that all of the Standard Model (SM) fields are 
free to propagate through all of the dimensions of space, rather than being confined to a brane 
as some or all of them are in the ADD and Randall-Sundrum models. Models featuring some or 
all of the SM gauge bosons and fermions in the bulk have appeared in the literature in various 
forms over the past decade, in particular in the work of Dienes, Dudas and Gherghetta [10] 
(see also Ref. [II]). Appelquist, Cheng and Dobrescu were the first, however, to recognize that 
the introduction of KK-parity considerably relaxes the constraints from electroweak precision 
observables, allowing for much lower scales of compactification. 

For a number of reasons, which we will discuss later, extra dimensions of size R ~ TeV -1 are 
particularly well motivated within the context of UED, much smaller than those found in the 
ADD model. In many ways the UED model is more similar to the original proposals of Kaluza 
and Klein than either the ADD or Randall-Sundrum scenarios. Although the motivations for 
UED are very different from those of Kaluza and Klein, each introduce extra flat and compact 
dimensions of space in which all types of particles are free to propag ate@ 

One aspect of UED models which has received quite some attention is their ability to provide 
a natural candidate for the dark matter of our universe. In particular, the lightest Kaluza-Klein 
state can be stable in UED models and, as we will see, naturally colorless, electrically neutral 
and produced in the early universe with an abundance similar to that of the measured dark 



1 For earlier work along similar lines, see Refs. [5][6][7]- 

2 For a recent review comparing and contrasting the ADD, Randall-Sundrum and UED models, see Ref. |12] , 
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matter density. For these and other reasons, Kaluza-Klein dark matter has become a popular 
alternative to neutralinos |13j and other commonly studied dark matter candidates |14j . 

In this review, we will discuss the dark matter and collider phenomenology of UED models. 
With the Large Hadron Collider soon to begin its operation, and numerous dark matter exper- 
iments being developed and carried out, we hope this article will be found useful in the coming 
years. 



2. Universal Extra Dimensions 
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2 Universal Extra Dimensions 

In this section, we introduce the basic UED model and discuss its motivations, the predicted 
mass spectrum of Kaluza-Klein states, and the current experimental constraints. 

The basic UED model is remarkably simple: a flat metric with one or more compact ex- 
tra dimensions. There is one complication to this otherwise very simple picture, however. 
Fermions in a number of dimensions larger than four are non-chiral in the case of the simple 
compactification on a circle. If the extra dimensions are orbifolded, however, this problem can 
be avoided. In the five dimensional case, an S 1 /Z% orbifold removes the unwanted fermionic 
degrees of freedom, allowing for the existence of chiral fermions. The orbifold compactification 
consists of projecting a circular extra dimension onto a line with two fixed points. In six or 
more dimensions, many orbifolding schemes become possible. The case of a T 2 /Z% orbifold in 
six dimensions has particularly attractive features, for example. 

Whenever possible, we will not commit ourselves to a particular number of extra dimensions, 
or choice of orbifold. In some cases, when conclusions can only be drawn in the case of specific 
scenarios, we will specify which class of UED models we will consider. 

2.1 Motivations for UED 

A wide range of motivations for the UED model can be found in the literature. We summarize 
below a few of the arguments we find especially compelling. 

Proton Stability 

In the Standard Model, there exist six-dimensional baryon and lepton number violating opera- 
tors which cause the proton to decay with a lifetime on the order of 



where M is the cutoff scale of the theory where quantum gravity effects set in. To evade the 
constraints on proton stability from Super-Kamionkande (r p > 10 33 s) [15j . the cutoff scale must 
be at least ~ 10 16 GeV. In many extra dimensional scenarios, the cutoff scale is well below this 
value. Although there have been several possible solutions proposed [31 [TBI [17], i n the original 
versions of the ADD and Randall-Sundrum models, protons are generically predicted to decay 
at a rate much higher than is experimentally acceptable. 

In the case of UED, the conclusion which is reached is rather different [18J. The global 
symmetries in UED can limit the higher dimensional operators which lead to proton decay. 
Considering the case of UED in six dimensions, for example, invariance under standard model 
gauge transformations and under the SO(l,5) Lorentz symmetry implies that the lowest dimen- 
sion operator which can lead to proton decay is 9-dimensional (in the 4-D effective theory). 
Higher dimensional operators can be shown to be generically suppressed |18| . and the leading 
proton decay mode to be such as p — > e~~ir + ir + vv, taking place with a lifetime of 



where R is the size of the extra dimensions and A is the cutoff of the theory. In this particular six- 
dimensional case, the proton decay rate is well below the limits imposed by Super-Kamiokande. 
More generally speaking, many of the operators leading to rapid proton decay are forbidden or 
suppressed in UED models, thus alleviating or at least relaxing some of the problems faced by 
the ADD and Randall-Sundrum models. 




(2) 
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2.1. Motivations for UED 



Similar arguments can be made to show that interactions leading to neutron-anti-neutron 
oscillations, lepton number violating interactions leading to large Majorana neutrino masses, 
and interactions leading to large flavor changing neutral currents (FCNC), can each be highly 
suppressed in UED. 

Anomaly Cancellation With Three Generations 

Motivation for UED has also been found from the requirements of anomaly cancellation. In 
particular, it was demonstrated in Ref. [19] that gauge anomalies can be avoided in UED models 
only if there exist three generations of fermions. In particular, in the case of six-dimensional 
UED, an SU(2)l global gauge anomaly exists unless the number of doublets with positive and 
negative chirality differ by an integer multiple of six. For each generation, this difference is either 
two or four, thus with three generations the global gauge anomaly is canceled in six-dimensional 
UED. 

Dark Matter 

As we will discuss at length in this review, the UED model naturally encompasses a stable, 
electrically neutral, and non-colored state which can constitute a viable candidate for dark 
matter. Although this was not among the first motivations for UED, it is one which has later 
received a great deal of attention. In fact, the UED model is the extra-dimensional scenario 
which most naturally provides a suitable particle dark matter candidate, to-date one the most 
compelling pieces of evidence for particle content beyond the Standard Model. 

Collider Searches 

Unlike other extra-dimensional scenarios, including models where only gauge bosons propagate 
in the bulk [20\ 12 lj . in UED KK parity suppresses the constraints on the size of the extra 
dimension(s) from electroweak precision tests, FCNC and other processes. As a result, KK states 
are generically light, making UED, in principle, a fully testable theory at present (Tevatron) 
and future (LHC, LC) colliders. This is especially true in the regions of parameter space favored 
by the abundance of KK dark matter produced thermally in the early universe. 

"Bosonic Supersymmetry" 

The UED model shares striking similarities as well as several relatively subtle differences with 
supersymmetry, the by far most studied TeV-scale new physics scenario. This led some to 
dub the UED scenario "bosonic supersymmetry", referring to the spin of the heavier, yet-to- 
be-discovered states, which, unlike in supersymmetry, feature the same spin as their Standard 
Model counterparts. Both in the realm of dark matter detection and collider searches, this 
has triggered a lot of work on how to distinguish between different new physics scenarios. For 
instance, the cascade decays allowed by KK number conservation have led to a better under- 
standing of the similar signals in supersymmetry, and have generated interest and a number of 
studies on the possibility of measuring spins at the LHC. 

A Non-Confining Strongly Coupled Theory 

The UED model provides an interesting example of an effective field theory which behaves, at 
the cut-off scale, as a non-confining strongly coupled theory. In particular this may lead to a 
dynamical origin for the SM Higgs doublet, along the lines of the analysis of Ref. [22]. The main 
ingredient to achieve this is the fact that fermionic fields propagate in the extra dimensions [23J . 
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2.2 The UED Model 

This section is devoted to a theoretical overview of the UED model. We start in Sec. !2.2.T1 with 
an outline of the Lagrangian of the theory and of its four-dimensional counterpart. Sec. 12.2.21 in- 
troduces Kaluza-Klein parity, a key feature for the model phenomenology. Finally, in Sec. 12.2.31 
we present the particle mass spectrum including the crucial ingredient of radiative corrections. 



2.2.1 The UED Lagrangian 



We consider the Standard Model (SM) in 4 + D space-time dimensions, where all the SM fields 
are allowed to propagate in the extra (universal) dimensions, compactified at a scale 1/R. 

Following Ref. [9], we indicate the usual non-compact space-time coordinates with the no- 
tation x^, fj, = 0,1,2,3, and the coordinates of the extra-dimensions with y a , a = 1,...,D. 
Indexes running over all space-time dimensions are indicated with capitalized letters. The 4- 
dimensional Lagrangian density is obtained by dimensional reduction of the (4+.D)-dimensional 
theory integrating over the compactified extra-dimensions as follows: 



£<7?) 
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i 



1 "-Hi 



rTr 



F t AB {x^y a )F iAB (x^y a ) 



+ 



+\{D, + D 3+a )H(x^y a )\ 2 + f j, 2 H*(x»,y a )H(x»,y a ) - \ [H*(x^,y a )H(x^,y a )] 2 + 
+i (Q, u, d, L, e) (a^, y a ) (T^ + T 3+a D 3+a ) (Q, u, d, L, e) (x», y a ) + 



Q(x^y a ) (\ u u{x^y a )io 2 H*{x^y a ) + \ d d{x^y a )H{x^y a )) + H.c 



+ 



L(x^,y a )X e e{x tM ,y a )H{x^,y a ) + H.c. 



(3) 



In Eq. (J3j) 5 the summation over fermion generations has been suppressed, and we indicate 
with Ff- B the (4 + Z))-dimensional gauge field strengths associated with the SU (3) c x SU (2)w x 
U(l)y gauge group. D M = d/dx^ —A^ and D 3+a = d/dy a — A 3 + a are the covariant derivatives, 
with Aa = —i Yli=i di^Ai^l being the (4 + D)-dimensional gauge fields, and gi the (4 + D)- 
dimensional gauge couplings. The latter, as well as the Yukawa matrices, \ Uj d,ei have dimension 
(mass)" / 2 . The symbols Q,u,d,L,e describe the (4 + D)-dimensional fermions, whose zero 
modes correspond to the SM fermions. Capitalized letters indicate SU{2)\\r doublets, while 
lower case letters indicate SU(2)w singlets. The (4 + L>)-dimensional gamma matrices, T A , are 
anticommuting 2 K+2 x 2 K+2 matrices, where D = 2K if D is even, and D = 2K + 1 for odd 
D, satisfying the (4 + d)-dimensional Clifford algebra {P 4 ,r ! ' B } = 2g AB . In particular, for the 
case of D = 1, one can set = 7 M , and T4 = 275. 

The last step needed to extract the effective four-dimensional Lagrangian from Eq. ([3]) is 
to specify the compactification of the extra-dimensions. The requirement of producing chiral 
fermions in four dimensions forces one, on general grounds, to consider orbifold compactifica- 
tions [9]. We consider here an [(S 1 x S l )/Z2\ K orbifold for D even, and an [(S* 1 x S 1 )/Z2\ K x 
(S 1 /Z 2 ) orbifold for D odd. At this point one must specify how the fields transform under the 
orbifold projection. Choosing these transformation appropriately allows one to eliminate the 
extra, phenomenologically unwanted, degrees of freedom appearing at the zero modes level. For 
definiteness, we shall focus below on the case of D = 1. 

Requiring that the A^ components of the gauge fields are even under the y — > —y transfor- 
mation, and that A§ is odd, will leave the zero modes corresponding to the physical SM gauge 
fields and will eliminate the extra zero modes (which would appear as massless scalar fields after 
dimensional reduction). Analogously, to retain the (zero mode) SM Higgs, one has to require 



8 



2.2. The UED Model 



that H is even under the orbifold transformation. With these properties, we can derive the 
following decomposition of the gauge and scalar fields in KK modes: 



j oo 
\/7T hi L 




n=l 



n=l 



(4) 
(5) 



The normalization is chosen once we assume the range for the y variable to be from to nR. 
The case of the fermions requires taking into account a subtlety related to the problem of chiral 
fields in 5 dimensions. Since one cannot construct, in 5 dimensions, a matrix with the properties 
of 75 in 4 dimensions (i.e. anti-commuting with all T A , and whose square is the identity), there 
is no chirality in 5 dimensions. Practically , this means that bilinears like V'T^TsV' are n °t 
invariant under 5D Lorentz transformations, so they cannot appear in the 5D Lagrangian. 
Consequently, one cannot have the left and right-handed components of the zero components 
of the zero excitation, ipQ, couple differently to the gauge fields, as in the SM. Therefore, one 
cannot get a chiral SM fermion using a single 5D fermion field, and is forced to introduce two 
5D fermionic fields for each Dirac fermion field [24J. Given a SM Dirac fermion field, ifj SM , 
one introduces two 5D fermion fields, ipi,R, each with the quantum numbers of left and right 
handed spinors, V'f^R- O ne can then set tp SM = Pl^l,o + PR^Rfi, with Pl,r = (1 T 7s)/2 the 
4-dimensional chiral projection operators. The generic KK decomposition will then be cast, for 
any given quark and lepton field, as 



/irR 



V SM K) + V2 ]T P L ^ L , n (x») cos (^) + P R ^ R , n (x») sin (!| 



n=l 



(6) 



As a result, the higher KK modes are 4-dimensional vector-like fermions, while the zero modes 
are chiral. This will induce a set of special Feynman rules, which we review and list in Appendix 

EU 



2.2.2 Kaluza-Klein Parity 

As the number of the KK level of a particle is a measure of its momentum in an extra di- 
mension of space, one might expect KK-number to be a conserved quantity by the virtue 
of extra-dimensional momentum conservation. In UED models, however, the introduction of 
orbifold compactifications breaks this translational symmetry along the extra dimension and, 
through Noether's theorem, leads to KK-number violating interactions. Even after an orbifold 
is introduced, however, a subgroup of KK-number conservation known as KK-parity can re- 
main present which insures the conservation of the "evenness" or "oddness" of KK number in 
an interaction [Sj. For instance, in the case of one extra dimension compactified on an S 1 /Z% 
orbifold, KK-parity is a Z2 symmetry under which only KK modes with odd KK number are 
charged. 

Phenomenologically, KK-parity acts similarly to how R-parity conserves the "evenness" or 
"oddness" of the number of superpartners within the context of supersymmetry. In general, 
KK-parity is conserved in UED if no explicit KK-parity violating interactions are introduced 
on the orbifold fixed points. To determine whether this is the case, one would have to consider 
a specific compactification scheme and have a full UV completion of the theory. Although KK- 
parity appears to be fairly natural within the context of UED models, without a UV completion 
we are unable to address with certainty whether KK-parity will be conserved. Throughout this 
review, however, we will proceed under the assumption that KK-parity is conserved. 
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KK-parity can be written simply as P = (— l) n , where n denotes the nth KK mode. There 
are several important phenomenological consequences of KK-parity conservation. Most impor- 
tantly for dark matter, KK-parity insures that the Lightest Kaluza-Klein Particle (the LKP, as 
it is often called) is stable [251 US]- Of consequence for collider experiments is that odd level 
KK states can only be pair produced [9j ETJ [28l [29] . Finally, all direct couplings to even number 
KK modes are loop-suppressed. 



2.2.3 The Particle Spectrum 



In this section, we describe the spectrum of Kaluza-Klein (KK) states in a 5-dimensional UED 
model, including radiative corrections and the effect of boundary terms. Beginning at tree level, 
all of the Standard Model fields appear as towers of Kaluza-Klein states with masses of: 

n 2 



m 2 x(n) 



n- 



(7) 



where X^ is the rath Kaluza-Klein excitation of the Standard Model field, X, and R ~ TeV -1 
is the size of the extra dimension. X^> denotes the ordinary Standard Model particle (known 
as the zero mode). Assuming that is considerably larger than any of the Standard Model 
zero mode masses, this leads us to expect a highly degenerate spectrum of Kaluza-Klein states 
at each level, n. This picture is somewhat modified, however, when the effects of radiative 
corrections and boundary terms are considered. 

Corrections to KK masses are generated by loop diagrams traversing around the extra 
dimension, called bulk loops, and by brane-localized kinetic terms which appear on the orbifold 
boundaries. Together, these contributions (at one loop) are given by [30J: 
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(8) 



Here 03) ~ 1.2020... is the third zeta function, and A is the cutoff scale of the theory, defined as 
the scale where the effective 5-dimensional theory breaks down, that is where the 5-dimensional 
couplings become strong and the theory is no longer perturbative. A is a parameter of the 
theory; to estimate its largest possible value, one can consider the loop expansion parameters 



(El) 



where the on are the 4-dimensional standard model gauge couplings, the index i = 1,2,3 labels 
the corresponding gauge groups, and iVj indicates the corresponding number of colors. Finally, 
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2.2. The UED Model 



AR counts the number of KK modes below A. Although the 4-dimensional SU(3) coupling 
becomes more asymptotically free above each KK level, the 5-dimensional SU (3) interaction 
becomes non-perturbative in the ultraviolet before the other gauge interactions. In particular, 
£3 becomes of order one for values of the cutoff scale typically of the order of A ~ 10 i? -1 . 

The brane kinetic term contributions in ([8]) are those proportional to the logarithm of the 
cutoff. In the case of KK scalars and spin-1 fields, the (n/R) 2 term in Eq. ([7]) gets replaced by 
(n/R) 2 + S(m 2 ), while in the case of KK fermions (n/R) 2 — > (n/R + 5(m)) 2 . KK gravitons do 
not receive appreciable radiative corrections, and their masses are simply given by n/R. 

After the KK modes of the W and Z bosons acquire masses by eating the fifth components 
of the gauge fields and Higgs KK modes, four scalar states remain at each KK level. These 
modes have masses given by: 

2 n 2 2 , r 2 
m H n = -ft2+ m h + Sm Hn 

2 n 2 2 . r 2 
m iT± = ^2 + m W + 5m H n 

m A n = fi2+ m Z + t™ 2 H n > (10) 

where the radiative and boundary term corrections are given by 

Sm 2 Hn = -£^2 (^g 2 + \g' 2 - 2A^) In AR + m 2 H . (11) 

Here, \h is the Higgs quartic coupling and rhn is the boundary mass term for the Higgs mode. 
Lastly, the additional contribution from the top quark Yukawa coupling yields: 

<W( m t(™)) 

In Fig. [H an example of the first level Kaluza-Klein spectrum is shown. The choice of cutoff 
used (A = 20 R" 1 ) leads to approximately 20% corrections to the KK quark and KK gluon 
masses, and corrections of a few percent or less for the rest of the first level states. 

The above expressions have been calculated under the assumption that the matching con- 
tributions from the brane-localized term cancel at the scale of the cutoff, A. Such contributions 
could, in principle, compete or even dominate over those corrections shown in Eqs. ([8]), (|lip 
and (I12p . With this in mind, it is possible that the KK spectrum might not resemble the re- 
sults described here. If such non-symmetric boundary term contributions are relatively small, 
however, the spectrum should closely reflect that described in this section [301 129] . 

The observation that radiative corrections to the mass parameters of the UED theory are 
naturally expected to lie at the cutoff scale - as defined in the usual sense of naive dimensional 
analysis as the scale for which loops are unsuppressed - leads to a little hierarchy problem in 
the context of UED [31J. For example, the zero mode Higgs field would naturally have mass 
at the cutoff A S> 1/R. Possible solutions to this fine-tuning problem between the cutoff scale 
and the compactification scale arise in the context of extra-dimensional versions of Little Higgs, 
or Twin Higgs models, where a symmetry in the bulk forbids the potential for the Higgs field. 
The symmetry is then broken in such a way as to generate a potential leading to a mass of 
the order of the weak scale. In the Little Higgs model case, the phenomenology of the theory 
might differ significantly from that of the simplest UED model, in that the compactification 
scale approximately coincides with the scale where the masses of the new states (the "Little 
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Figure 1: The spectrum of first level Kaluza-Klein states, including the effects of radiative 
corrections and boundary terms. A compactification radius of -R -1 = 500 GeV, Higgs mass of 
120 GeV and cutoff of A = 20 R' 1 have been used. From Ref. [29]. 



Higgs partners") appear. The latter scale is instead considerably higher for the Twin Higgs 
model, and the expected phenomenology (unless both scales are significantly low) is similar to 
that of UED 01]. 

As we discussed in Sec. 12.2,21 the lightest of the first level KK states can be stable in UED. 
Since we are particularly interested in dark matter within the context of these models, we would 
like to determine which of the first level KK states is the lightest and, therefore, potentially 
stable. The results of Ref. [30J would lead us to conclude that the KK photon will be the 
lightest KK particle (LKP), due to the small (and negative) correction to its mass. Things 
are complicated somewhat, however, by the fact that the KK photon is a mixture of the 
and VF 3 W fields, just as the zero-mode of the photon is a mixture of and W 3( -°\ The 
corresponding mass matrix, in the B^ n \ W 3 ^ basis, is as follows: 

^ + <K (n) + ^V h'av 2 

Here v » 174 GeV is the Higgs vacuum expectation value. In the zero-mode case (n = 0), the 
off-diagonal terms are comparable in magnitude to the diagonal terms, leading to significant 
mixing between and 1U 3 ^ ^ (sin 2 6\y ~ 0.23). At the first KK level, however, the diagonal 
terms are much larger than the off-diagonal terms, leading to only a small degree of mixing 
between the B^ and VF 3 *- 1 -* states. For R~ l =\ TeV, for example, the effective first KK level 
Weinberg angle is approximately sin 2 #14/1 ~ 10~ 3 . Thus the mass eigenstate often called the 
"KK photon" is nearly identical to the state B^\ 

If the corrections given by Eq. ([8]) do not fully describe the contributions to the KK spectrum, 
it is plausible that, instead of B^\ another of the first level KK states could be the LKP. 
Limiting our choices to electrically neutral and non-strongly interacting states, we are left with 
the possibilities of z/W, TU 3 ^ 1 ) or the first level KK graviton for the LKP. KK neutrinos were 
considered as a possible choice for the LKP in Ref. [25], but were found to generate unacceptably 
large rates in direct detection experiments in Ref. [32J. As for the KK mode of the neutral SU(2) 
gauge boson, the SM contribution to its mass makes it difficult to envision a scenario where it 
is lighter than the U(l) KK gauge boson. 
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A latter possibility is that the LKP is the KK graviton (for a discussion of when this is 
the case in the context of the minimal UED model, see e.g. the analysis of Ref. [33] )• The 
cosmological consequences of a KK graviton LKP might turn out to be rather daunting. In 
particular, KK gravitons can be overproduced in the early universe, and can also distort the 
process of light elements nucleosynthesis, the cosmic microwave background, as well as the 
diffuse gamma-ray background. Even if the KK graviton is not the LKP, the presence of a KK 
graviton tower can potentially lead to similar problems. A range of reheating temperatures 
does exist, however, for which these problems can be avoided [34J . We discuss the physics of 
KK gravitons in some detail in Sec. 13.7.21 

Throughout the majority of this review, following most of the phenomenological analyses 
carried out so far in the literature, we will assume that the is the LKP. 

2.3 Electroweak Precision Constraints 

Very precise measurements have been made of a number of Standard Model observables, in 
particular in the electroweak sector. These measurements agree remarkably well with with the 
predictions of the Standard Model. Any model of physics beyond the Standard Model must be 
able to conform to this very restrictive body of data. 

To illustrate how restrictive electroweak precision data can be, consider an extra-dimensional 
model in which only the gauge bosons are free to propagate in the bulk. In such a model, LEP 
data restricts these Kaluza-Klein gauge bosons to have masses far above the electroweak scale, 
at least as heavy as a few TeV [201 EI] ■ 

Such constraints are weakened, however, in models in which all of the fields of the Standard 
Model are allowed to propagate in the bulk (UED models) [91 [351 E3 EZ] . In this class of models, 
the strongest constraints come from measurements at the Z-pole, including the relationship 
between the Z and masses. 

The full independent set of electroweak precision observables (EWPO) can be defined by: 
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where II denote the new-physics contributions to the transverse gauge boson vacuum polariza- 
tion amplitudes, with n'(0) = dH(q 2 )/dq 2 \ q 2 =0 , etc. The parameters S, T and U are related to 
the usual S, T and U by S = 4 sin 2 9wS/a, T = T/a and U = —4 sin 2 6v/U /a. Of these seven 
quantities, all but V have been well determined by LEP data, and provide stringent constraints 
on new physics models. 

Within the context of UED, the strongest constraints come from the parameter T, and 
to a somewhat lesser extent from S. In particular, these parameters acquire leading order 
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Figure 2: Left: The magnitude of the contributions to the parameter T from each of the first 
three KK levels (dashed) for txikk = 400 GeV as a function of Higgs mass, as well as the sum 
over the first 10 KK modes (solid). The magnitude of the Higgs-dependent correction (the 
second term in the expression for T in Eq. (|15p ) is shown as the dotted line (first from above). 
Right: The 95% (dashed) and 99% (dotted) confidence limit exclusion regions, as a function of 
Higgs mass and mass ttikk = 1/-R. From Ref. |36j . 



contributions from the presence of KK modes given by: 

g 2 m\ hg 2 s\y? d w m 2 h 



96m 2 v mj <K 1152 cos 2 6\vfn 2 KK ' 
* _ g 2 m 2 g 2 m 2 

576m 2 KK 2304m 2 KK ' { ' 

Comparing these two expressions, we see that the contributions to T are larger than those to S 
by a factor of ~ Qmf/m^. Given that the experimental constraints on T and S are comparable, 
it is clear that the parameter T will dominate the electroweak precision constraints on UED 
models. 

In the left frame of Fig. [21 the size of the contributions to T are given for the first three KK 
modes for the case of tukk = 400 GeV, as a function of the Higgs mass. Also shown are the 
contributions from the zero modes alone, which depend significantly on the Higgs mass. The sign 
of these two classes of contributions are opposite enabling a cancellation to occur. Summing the 
contributions, the 95% and 99% confidence limit exclusion regions are shown in the mKK- m h 
plane in the right frame of Fig. [2j For a light Higgs mass, limits as strong as ra^K ^ 700 — 800 
GeV can be obtained. For a heavy Higgs (m^ > 300 GeV), electroweak precision observables 
constrain the viable range of from below and from above: for instance, at 95% C.L., 
400 < J?- 1 /GeV < 600 for m h = 500 GeV, and 300 < R^/GeV < 400 for m h = 800 GeV [36]. 



2.4 Accelerator Searches: the Current Status 

Accelerator searches for extra dimensions have been performed for a variety of theoretical setups. 
The most recent and strongest constraints on models with large extra dimensions and warped 
extra dimensions were obtained at the Tevatron. Those results are summarized in Refs. [38} 139] . 
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The only existing analysis of direct searches for KK excitations in the UED scenario can be found 
in Ref. |40j . and is based upon CDF Run IB data taken in the years 1994-1996 1 1 1 j . This analysis 
relies on the multi-lepton channel only, and is restricted to the case of one extra dimension. Also, 
it is assumed that the spectrum of KK particles follows the radiative corrections we summarized 
in Sec. 12.2.31 with rh^j, the boundary mass term for the Higgs mode, set to zercH. 

The heavier n = 1 states cascade decay into the stable LKP, B^\ emitting soft SM particles. 
The LKP then escapes detection, leading to a missing energy signature. The KK modes which 
should be most abundantly produced at a hadron collider are KK quarks and gluons. KK parity 
dictates that KK states are pair produced (see Sec. [3]). The production cross sections for KK 
modes at the Tevatron were first computed in Refs. \27\ I28) 4 I . We show the cross section for 
the production of a pair of KK quarks and/or gluons at the Tevatron Run I in Fig. O Once 
produced, the heavier strongly interacting KK states follow a homogeneous pattern in their 
decay chains (unlike, e.g., the wide range of possibilities in the case of generic supersymmetric 
spectra). We show a typical instance in Fig. [H In the figure, the dominant decay channels are 
represented by solid lines, and the sub-dominant by dotted lines. 

The heaviest KK mode, the KK gluon, decays with almost 100% branching ratio (BR) to 
a KK quark and an ordinary quark. In particular, BR(g^ — » Q^Q) — 0.5, and BR(g^ — > 
u^u) + BR(g^ -> dWd) ~ 0.5 [29]. The SU(2) singlet KK quarks then decay into the zero- 
mode (ordinary) quarks and an LKP, B^\ resulting in a jets plus missing transverse energy 
signature. The SU(2) doublet KK quarks, on the other hand, will mostly decay to the 

zero-mode quarks and a KK or a gauge boson. The latter subsequently decay to 

KK leptons plus ordinary leptons. In their turn, KK leptons will decay to an ordinary lepton 
and a B^ l \ The resulting signature would then be 3 leptons plus 2 jets plus missing energy 
(from the W^'Z^' channel), or 4 leptons plus 2 jets plus missing energy (from the Z^Z^ 
channel). Fig. [5] illustrates the Feynman diagram leading to the multi-lepton plus 2 jets plus 
missing energy signature in KK SU(2) double associated production. 

The analysis carried out in Ref. [40] only considers the multi-lepton channels which con- 
tribute to the signal search in the Tevatron data, effectively disregarding final states which 
include taus. Furthermore, it was assumed in the Monte Carlo simulations for the signal that 
the relative branching ratio of KK SU(2) doublets into versus that into Z^ is 4 : 1. The 
resulting KK gauge bosons are then forced to decay into those channels included in the signal 
analysis (see Tab. 4.1 of Ref. [40]). 

Although the largest final state signature for UED at the Tevatron would be jets plus 
missing energy, the SM background suppression requires one to resort to the much cleaner multi- 
lepton plus missing energy signature. In this latter case, the main sources of background are 
bb/cc, tt, WZ and ZZ production, as well as Drell-Yan and WW production with misidentified 
leptons. The blind analysis procedure followed in Ref. [40] consists of identifying the data 
sample containing three leptons passing certain criteria in the CDF dataset, generating both 
the SM and the signal events through Monte Carlo simulation, and passing it through the CDF 
detector simulator. Selection criteria are then adopted to enhance the signal and to suppress 
the background, eventually obtaining an upper limit on the number of UED signal events at a 
given confidence level from the number of observed data events. Finally, this is translated into 
an upper limit on the production cross section, and on the size of the compactification radius. 

The search for three or more charged leptons and missing transverse energy was performed 
with a total integrated luminosity of 87.5 pb -1 collected during CDF's Run IB period, resulting 
in a production cross section limit for SU(2) KK quark pairs at 95% (90%) Confidence Level 
(CL) of 3.3 pb (2.5 pb). This translates into a lower limit on the inverse extra-dimensional 
radius of 1/R >270 GeV (280 GeV). A limit for the total KK quark pair plus KK gluon-quark 

3 A non-zero value for fh 2 H should not significantly affect the results of the analysis, though. 

4 We collect in App. [rOl the production cross sections for strongly interacting KK states at hadron colliders. 
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Figure 3: The production cross section of KK quarks and gluons at the Tevatron Run I (upper 
panel) and Run II (lower panel). The solid black curve represents the total production cross 
section, while the other lines show the separate contributions from KK quark pairs (red dotted 
line) KK quark-gluon (green dashed line) and KK gluon pairs (blue dot-dashed line). The KK 
top production cross section is indicated by the orange double-dot-dashed line. Adapted from 
Ref. [28] • 
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Figure 5: The Feynman diagram leading to the multi-lepton plus 2 jets plus missing energy 
signature in KK SU(2) double associated production. 
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plus KK gluon pair production cross section of 7.9 pb (6.0 pb) was found, which translates into 
a limit of 1/R >280 GeV (290 GeV). 

2.5 Other Phenomenological Constraints 

Physics beyond the SM can manifest itself not only via direct production of new particles 
at colliders, but also indirectly, e.g. contributing to precision electroweak data, as discussed 
in Sec. 12.31 or to other processes or observables such as rare decays, flavor physics, or the 
anomalous magnetic moment of the muorJl. In this section, we summarize the effect this latter 
class of phenomenological constraints has on the UED model. 

The inclusive b —* branching ratio is known to provide extremely stringent bounds on 
physics beyond the SM, for instance low energy supersymmetry [33]. The contribution of UED 
KK states to b — > was first evaluated in Refs. [444 145] . The effective Hamiltonian for AS = 1 
B meson decays is given by 

n eS = ^=-v tb v t * s ewe), (16) 

where, at leading order, the operator relevant for the transition b — > S7 is 

07 = J^mbSLaV^bRaF^. (17) 

In UED, the contributions to the relevant coefficient C-j come both from loops involving the 
KK top quark and (1) the KK tower and (2) the KK states of charged (physical) would- 
be-Goldstone bosons. Contribution (2) is dominant, and gives 

r UED r . ™? ( t> ( m 2 t +{n/Rf \ 1 ( m* + {n/Rf \\ 
7 ^\ J + Wl l HR? ) 6 V (n/Rf ))> 

where the loop functions are given by 




Combining theory and experimental la errors in quadrature, the 95% CL constraint on the 
contribution of KK states, with [/,(, — > mj, can be cast as 

|[^7° tal (^)] 2 /[C 7 SM (^)] 2 -l| <0.36 (21) 

In the limit of mt <C (1/R), Ref. [45J finds the approximate expression 

Y,m t /{m t + (n/R) 2 ) <0.5. (22) 

n 

Using the exact formulas, the resulting bound on the size of the inverse compactification radius 
from b — > 57 quoted in Ref. [15] reads 1/R > 280 GeV. A refined analysis of the UED contri- 
butions to the processes was then carried out in Ref. [46J . The latter includes a next-to-leading 
order (NLO) treatment of the QCD corrections, which has two important consequences: first, 



3 See e.g. the recent review (35] on low energy precision tests of supersymmetric models. 
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Figure 6: The branching ratio for b — > 57 and E~ > 1.6 GeV as a function of the inverse 
compactification radius, 1/R. The grey band represents the experimental result, the dashed 
lines the SM values, and the solid lines the UED results. The upper (lower) curves make use of 
a value for the ratio m c /mb = 0.22 (m c /mb = 0.29), and all of the curves feature an estimated 
theoretical uncertainty of order 10%. From Ref. [46] . 



depending on the scale \ib the enhancement of the NLO result with respect to the LO one can 
be larger than a factor 2; and, secondly, and perhaps most importantly, the dependence upon 
the scale is vastly suppressed with respect to the LO result, and amounts, for reasonable 
values of to not more than 1.5%. The final result of the analysis is shown in Fig. [U where 
the grey band represents the experimental result, the dashed lines the SM result, and the solid 
lines the UED result; the upper (lower) curves make use of a value for the ratio m c /mb = 0.22 
(m c /mb = 0.29), and all of the curves feature an estimated theoretical uncertainty of order 10%. 
As can be deduced from the figure, taking into account the theory uncertainty, one can conclude 
that consistency with b — > 57 implies 1/R > 250 GeV |46j . a value slightly more conservative 
than the one quoted in Ref. [45] . The more recent analysis carried out in [47], making use of 
the exclusive branching ratio B — > K*~f, confirmed the importance of the detailed knowledge 
of the relevant form factors, and shows that the experimental constraint, under conservative 
assumptions, can be translated into a bound of 1/R > 250 GeV. 

Another precision electroweak observable which is sensitive to physics beyond the SM is the 
Zbb vertex [48j [49j [50] . In the context of UED, corrections to the Zbb vertex were first considered 
in the seminal paper of Ref. [9], and then re-examined in Refs. [51] and [52]. Shifts in the Zbb 
coupling due to radiative corrections, either from the SM or from new physics, affect observables 
such as the branching ratio Rb = Tb/Th, where Tb = T(Z — ► bb) and Th = T(Z — > hadrons), or 
the left-right asymmetry, Aj,, defined as 

^fl^ff, (23) 



where, in turn, the couplings 



9L = -\ + S f + ^ + 5gT (24) 
9R = S f + Sg s R u +S g r (25) 
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enter in the Zbb vertex as 

—lrf{gLPL+9RPR)bZp. (26) 

Both within the SM and within most of its extensions, the largest corrections, scaling with the 
square of the top quark mass, affect the gi coupling. A shift Sg^ p then entails a shift in R b 
given by 

SR b = 2R b (l - R b )^kl_ 5g f^ (27) 
9l + 9r 

and a shift in A b of 

SAb = A 9 2 R 9L Sg xp (2g) 
(9l+9 R ) 2 

The dominant contributions to 5g^ p within UED come from loops containing charged KK 
Higgses and top quarks. Summing over all possible loop corrections yields the result [51] 

^ued „ V2Gft^R\ (29) 

Using the SM and experimental input values for R b , one gets from Eq. f|29|) a bound of 1/R > 
230 GeV from R b , and a much weaker bound from A b ; however, sub-leading contributions in 
m^/mf, stemming from loops involving W/t and , give an extra contribution which 

pushes the limit on the compactification inverse radius from R b up to 1/R > 300 GeV, according 
to the estimate of Ref . [51] • 

Loops of KK quarks and KK gauge and Higgs bosons also contribute to observables related 
to the unitarity triangle, as well as to rare K and B decays. The impact of UED on the — B® 
mixing was first addressed in Ref. [S3], and re-evaluated, together with the UED corrections to 
other quantities relevant to the unitarity triangle (including e.g. ex, AM S , AMk), in Ref. [52]. 
The conclusion of the analysis of Ref. was that with the current la experimental range on 
AMrf, the mass difference between the B® mesons, one can derive a constraint of 1/R > 165 
GeV. If future measurements of the B meson decay constant, and the QCD correction 
parameter, Bb, were able to reduce the errors on these quantities by one third (and yielded 
the same central value for the combination /bVPb), this could lead to a la bound on the 
inverse compactification radius as large as 1/R > 740 GeV. Ref. [52], however, pointed out 
that analyzing all UED contributions to the unitarity triangle (UT), values of 1/R ~ 200 
GeV are perfectly compatible with the present fits of the UT, and that an improvement in the 
accuracy of the determination of /bV^b is not enough to derive more stringent bounds on 1/R. 
Rather, an improvement on both Jb^/Bb and \Vtd\, if determined through AM^/ AM S , can yield 
an improved maximal or minimal bound on 1/R sensitively depending upon the maximal or 
minimal value of /bV^b- 

Ref. [52] pointed out that the short distance one-loop functions relevant for the rare decay 
modes K + — ► 7t + vv, Kl — ► iftvv and B — * X s ^vv are larger than the SM value by around 
10%, when the compactification scale is set to 1/R ~ 300 GeV, due to KK contributions to 
the Z penguin diagrams. For the decay modes K + — > tt + uu, Kl — > 7r°i/P and B — ► X^vv this 
enhancement is partly compensated by the fact that |V^|ued < |T^td|sM> an d that r/uED < ^sm> 
resulting in the enhancement of the branching ratios for the aforemenitoned rare decays with 
respect to the SM expectation of 9%, 10% and 12% respectively, for l/Rm 300 GeV [52J. As 
B — > X s uv is governed by the CKMM matrix element \Vt s \, that is common to the SM and 
UED, the projected enhancement of the BR over the SM prediction is around 21% [52]. 
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Analogously, the short distance one-loop function relevant for the B^ s and Kl — > 

decay is enhanced by around 15% in UED. This results in a 33% enhancement of BR(B S —* 
M + M~)> 23% of BR(B d -> and 20 % of BR(K L -» for the UED case over the 

SM. The analysis of the exclusive decay modes 5 — ► Kl + l~ and 5 — > K*l + l~ carried out 
in Ref. [47] confirms the enhancement with respect to the SM, and point out that in general 
no constraints can be drawn on 1/R from those decay modes. The subsequent analysis of 
Ref. [54] considered the branching fractions and the lepton polarization asymmetries for the 
rare decay modes B — ► X s t + t~ - and B -> t + t , focusing in particular on the transverse 
asymmetries. They find that for the exclusive transitions the uncertainties due to the hadronic 
matrix elements are small, and in the large energy limit for the light hadron in the final state, 
the polarization asymmetries are basically free of hadronic uncertainties. 

The branching ratio B — ► A s /x + /i~ has also been investigated in Ref [46] . Interestingly 
enough, UED enhance the SM result by a factor of up to 12% (for 1/R ~ 300 GeV), which 
goes in the right direction, as the data from BELLE [55J indicate that the measured branching 
ratio is larger than the SM expectation. In view of the current experimental and theoretical 
uncertainties, however, this is certainly not enough to draw any significant constraint on the 
UED model. However, it was pointed out in Ref. [16] that a quantity of particular interest 
for the discrimination of a signature from UED is the forward-backward asymmetry in the 
B — > X s fj, + fi~ decay, defined as 



and 9i is the angle between the fi + and B meson momenta in the center of mass frame. The 
location of sq, the zero of Ap^(s) is a particularly good observable, as it is sensitive to short 
distance physics, and it is subject to very small non-perturbative uncertainties. This motivated 
the authors of Ref. [46J to state that future experimental results on Afb(s) could lead to 
particularly strong tests of the UED model. This point of view has been confirmed by the 
recent analysis of Ref. [IT]- Furthermore, the experimental results of the BELLE collaboration 
indicate that the relative signs of the Wilson coefficients CV and Cg is negative, implying that 
Afb should have a zero [56J. An analogous analysis focusing on the exclusive A decay was 
carried out in Ref. [57], where it was shown that for 1/R ~ 300 GeV the branching ratio 
is enhanced by around 20% with respect to the SM result, and where the importance of the 
forward-backward asymmetry was also emphasized. Ref. [58] studied double-lepton polarization 
asymmetries in A& — ► M + l~ decays. Ref. |58j showed that, unlike the case of single-lepton 
polarization, where a discrimination between the SM and UED does not seem feasible, various 
double-lepton polarization asymmetries in the above mentioned decays can be useful tools for 
establishing new physics signals as predicted in the UED model, as the predictions for the SM 
and for UED differ substantially, and to a level that might be experimentally accessible in the 
foreseeable future. 

The branching ratio for the rare decay mode Kl — > 7r°e + e _ is enhanced in UED up to 
10% more than its SM value [46j ; the latter is however still more than 2 orders of magnitude 
smaller than the current experimental sensitivity [59|. More promisingly, the e'/e ratio, which 
parameterizes the size of the direct CP violation with respect to the indirect CP violation 
in Kl — > tttt decays, is typically of the same order, or larger, than the current experimental 
uncertainty for 1/R < 400 GeV. However, Ref. [M] pointed out that, even fixing a specific value 




(30) 



where 



s = 



(31) 
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1 /2 3 /2 

for m s (m c ), the uncertainties on the main non-perturbative parameters B & and B s is such 
that no constraints can be drawn model-independently from the experimental results and the 
UED expectations. 

An interesting by-product of the analysis of Ref. [52] is that B and K physics can, in 
principle, distinguish between UED and super symmetry. For instance, the expectation for the 
MSSM in the large tan/3 regime is a suppression of AA1 S with respect to the SM prediction, 
while a generic prediction of UED is an enhancement: the ratio AM| xp /AMf M can thus indicate 
a preference for one option rather than the other. Similarly, in the regime of small tan /3, the 
MSSM predicts a suppression of the branching ratios for K + — > -k + vv, Kl — > ifivv, B — > X^yv 
and Bd — > while an enhancement is predicted in the UED scenario, as pointed out above. 

Among other flavor changing neutral current effects considered in Ref [46] , the UED predic- 
tions for the process b — > sg show very significant departures from the SM result as well as, once 
again, a strong enhancement due to NLO QCD corrections [60]. However, both the SM and the 
UED results strongly depend upon //& and on the ratio m c /mfo. In principle, a good determina- 
tion of the two quantities could allow for a discrimination of UED effects up to 1/R ~ 600 GeV. 
Unfortunately, however, a significant reduction of the uncertainties related to //& and to m c /m& 
looks problematic, and the extraction of the branching ratio from the experimental data also 
appears very difficult, if not impossible [60| IB"!]. 

The role of the muon {g — 2) for UED was first considered in Refs. [44j . |62j and |27| . 
The contributions to the anomalous magnetic moment of the muon in UED arise from loops 
involving KK gauge bosons or KK scalars and KK muons or muon neutrinos. The contribution 
of the n-th KK level to the muon {g — 2) reads 

f 5e 2 m 2 3e 2 (24 - l)m 2 f\ p-* 6x 2 + Uxy - Ux 

7777? 



128wHn/Ry 32^cP w (n/Ry J J \ + {x + y) 



leading to the approximate result, expanding in powers of the small parameter e = (mzR) 2 , of 

[23 

- 439 x 10~ n (e - 0.23e 2 + ...) . (33) 

The first thing to notice is that the overall sign is opposite to the potential difference between 
the experimental data and the SM prediction. Secondly, the absolute value of the contribution 
of UED is at most (g — 2)tjed — — 4 x 10~ 10 which is only one quarter as large as the SM 
electroweak contribution, and is smaller than the final expected sensitivity of the muon {g — 2) 
experiment at BNL |63| . In addition, the independent assessment of (g — 2)tjed carried out in 
Ref. [62j finds a final result that is smaller by a factor of approximately 3 than that found in 
Ref. [27J, corroborating the conclusion that the muon {g — 2) does not yet provide any useful 
constraint on the UED scenario. 

As a final remark, we wish to comment on the interesting recent result of Ref. [B3], where 
it was pointed out that the charged lepton flavor violating (CLFV) decays Zj — > lj"f in UED 

are suppressed relative to the SM result by factors of (ffliy/m^ ) 8 , making them an absolutely 
hopeless channel to observe effects of the presence of UED. Five-dimensional UED can therefore 
produce observable lepton flavor violation only if explicitly flavor- violating boundary terms are 
assumed to be present. Finally, CLFV decays in the context of UED with more than one 
extra-dimension were discussed in Ref. [65]. 

We give a summary of the constraints on the size of the extra-dimensional radius from flavor 
physics and other low-energy probes of new physics discussed in this section in Tab. [H 
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Limit on 1/R 


Ref. 


BR(6 -► s 7 ) 
BR(5 -> K*i) 
R ( r(z->bb) \ 

Ji-fe ^r(Z-^hadrons) J 

Rb (plus sub-lead, contrib.) 


> 250 GeV 

> 250 GeV 

> 230 GeV 

> 300 GeV 


m 
m\ 


AMd (m B o — mgo ) 

\ d d J 


> 165 GeV 


[53] 



Table 1: A summary of the constraints on the 5-dimensional UED model from rare decays and 
flavor physics. 



3. Kaluza-Klein Dark Matter 



23 



3 Kaluza-Klein Dark Matter 

Particle dark matter candidates must fulfill an array of phenomenological constraints pertain- 
ing to their interactions with ordinary matter, dictated by the results of dark matter search 
experiments and other considerations. In particular, in the range of particle masses of relevance 
here (above ~ 300 GeV and below or around the TeV scale) a dark matter candidate which is 
charged under strong or electromagnetic interactions would become bound with ordinary mat- 
ter, forming anomalously heavy isotopes. Negative results of searches for anomalously heavy 
isotopes lead to the exclusion of any strong or electromagnetically interacting dark matter par- 
ticles in the mass range of relevance here [66, 67J. This rules out the KK gluon, all KK quarks 
and the charged KK fermions, as well as the charged KK Higgs and gauge bosons, as the LKP. 

Neutral KK SU(2) gauge bosons and Higgs bosons are unlikely to be the lightest KK particle 
since their masses squared are shifted by their large zero mode masses with respect to KK 
neutrinos and the KK photon (see Sec. l2.2.3l) . Although departures from the spectrum described 
in Sec. 12.2.31 could potentially change this conclusion, it is somewhat difficult to envision a 
theoretically motivated setup where the bulk corrections and the boundary terms overturn this 
natural hierarchy. As shown in Ref. |33| . one such scenario is possible if the SM Higgs mass is 
assumed to be very heavy (m^ > 230 — 250 GeV). In this case, depending upon the value of 
the compactification radius, and assuming, again, a minimal UED setup for the radiative mass 
corrections, the LKP can be the charged Higgs boson, iJ^M. This latter option, as explained 
above, is phenomenologically excluded, however, unless special, very low reheating temperature 
ranges, Trh < 1 GeV, and large compactification radii, R > 1 TeV, are assumed [33]. In any 
case, the H^ 1 ^ relic density would be insufficient to be a significant contributor to the overall 
dark matter budget. 

The only KK fermion which could play the role of a dark matter particle is the KK neutrino. 
Relaxing the restriction that SU(2) charged fields are heavier than their singlet counterparts 
(as it is the case in the minimal setup [30]), one can consider the possibility that one of the KK 
neutrinos is the LKP [25j . The relic abundance of KK neutrinos was computed in Ref. |25j in 
four different setups: with one and three neutrino flavors, and with or without coannihilations 
with the left handed KK leptons, taken to be degenerate in mass with the KK neutrinos. 
While most coannihilations play a relatively small role here, it turns out that the number of 
flavors is a critical input, as extra Z^> t-channel exchanges can occur between neutrinos with 
different flavors. As a result, the overall effective cross section is increased, and the range 
of masses compatible with the observed cold dark matter abundance is shifted toward larger 
values. Namely, the acceptable mass range for one KK neutrino flavor is between 0.8 and 1.0 
TeV, while with three flavors the favored mass range lies between 1.1 and 1.3 TeV. 

As discussed in Sec. 12.2.31 the constraints on a KK neutrino LKP from the negative re- 
sults of elastic scattering experiments such as CDMS and EDELWEISS are particularly strong. 
Quantitatively, we find that the current CDMS limits [68J and the results of Ref. [32] on the KK 
neutrino-nucleon scattering cross section imply that a KK neutrino must be heavier than around 
1200 TeV if it is to be a viable particle candidate to make up the dark matter of our universe. 
This is clearly at odds with the range favored by thermal relic abundance considerations. As 
is well known, thermal production in a standard cosmological setup is not the only mechanism 
available to produce the dark matter in the early universe. Alternatives include non-thermal 
production and late time entropy injection. Such alternatives can in principle enable a very 
heavy relic, which would be produced in the standard scenario in excessive amounts, to be 
brought into accord with the desired cold dark matter abundance (for example, see the recent 
analysis of Ref. [69]). Nevertheless, an inverse compactification radius of O(10 3 ) TeV, although 
phenomenologically viable, is certainly not theoretically favored, nor welcome from the point of 
view of detecting KK states at colliders or at dark matter searches other than elastic scattering 
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Figure 7: The process of the thermal freeze-out of a stable, weakly interacting particle in the 
early universe. As the temperature drops below the particle's mass (x > 1), the number density 
of such particles becomes Boltzmann suppressed. As the universe expands further, eventually 
these particles encounter no others of their species with which to self-annihilate, leading their 
density to "freeze-out" . 



experiments. We therefore conclude that a KK neutrino LKP is at present not a favored option. 

The choice of LKP which has received the greatest interest is by all means the B^, some- 
times called the KK photon. In this section, we discuss at length the phenomenology of a 
LKP from the point of view of its implications for cosmology and dark matter searches. It is 
worth bearing in mind, however, that another type of LKP may also constitute a viable dark 
matter candidate, the KK graviton. We devote Sec. 13.71 to a discussion of the cosmology and 
phenomenology of KK gravitona^). 



3.1 Relic Abundance 



In UED models, KK dark matter (KKDM) states are abundant in the early universe (T ^> 
R^ 1 ~ TeV), being freely created and annihilated in pairs. As the universe expands and the 
temperature drops below that needed to produce such states in chemical equilibrium, however, 
the number density becomes rapidly suppressed. A certain density of stable KK states (here, 
the B^s) fr eezes out, and remains in the form of a thermal relic of the universe's hot youth. 
In this section, we review the calculation which is performed to determine the thermal relic 
abundance of KK states in the universe today [71] . 

The number density of B^s evolves according to the Boltzmann equation: 
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(34) 



where H = ^y8Trp/3M p \ is the Hubble rate and (av ) is the B^s self-annihilation cross section. 
The equilibrium number density of B^s is given by: 
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6 Another possibility for dark matter within the context of extra dimensional models in which scalar fields are 
allowed to propagate in the bulk is that of soliton-type states. The existence and stability of such states was 
demonstrated in Ref. [70] . The lightest such state might feature a mass around a TeV, be electrically neutral, 
stable, and in principle be a viable dark matter candidate. 



3. Kaluza-Klein Dark Matter 



25 



where g = 3 is the number of degrees of freedom of the LKP. At T > rn B m , the number density 
of B^s was very close to its thermal equilibrium value. As the temperature dropped below 
m B (i), the number density became exponentially suppressed, until eventually the annihilation 
rate was overcome by the effects of Hubble expansion (see Fig. [7]). 

Numerical solutions of the Boltzmann equation yield a relic density of |71j : 



2 1.04 x 10 9 x F 

BW ~ M P1 ^(a + 36/x F )' (36) 

where = m B (i)/Tp, Tp is the temperature at freeze-out, 5* is the number of relativistic 
degrees of freedom available at freeze-out (g* ~ 92 for the case at hand) , and a and b are terms 
in the partial wave expansion of the LKP annihilation cross section, av = a + bv 2 + $(v ). 
Evaluation of xp leads to: 



xp = In 



,,. / 45 g m B( i)Mpi(a + 6b/xp) 



(37) 



where c is an order 1 parameter determined numerically. Note that since xf appears in the 
logarithm as well as on the left hand side of the equation, this expression must be solved by 
iteration. WIMPs with masses near the electroweak scale generically freeze-out at temperatures 
in the range of approximately xp ~ 20 to 30. 

The B^ self-annihilation cross section is roughly constant as a function of temperature, 
and is approximately given by (see Appendix IB.2jl : 

(ov) - J 5g { ■ (38) 
3247rm^ {1) 

This is unlike what is generically expected in the case of supersymmetry: when the lightest 
neutralino is bino-like, the pair-annihilation into gauge or higgs bosons pairs final states is 
suppressed, and the main annihilation channel is into fermion pairs through sfermion t-channel 
exchange. As a consequence of the Majorana nature of neutralinos, fermion pair (//) final 
states feature an s-wave amplitude for the pair annihilation cross section suppressed by a factor 
TO^/m^, where m x is the neutralino mass [13]. Therefore one expects a sizable variation of the 
thermally averaged pair annihilation cross section with temperature. In addition, the mentioned 
suppression affects the pattern of pair annihilation final states branching ratios at T = 0, 
relevant for indirect detection: namely, independently of the composition of the neutralino, 
prompt annihilation in lepton pairs is always greatly suppressed. This is unlike the case 
LKP, where // final states are not helicity suppressed, but, rather, feature a sizable branching 
ratio. In the approximation that all heavier level one KK modes have the same mass, the 
relative annihilation fraction can be determined from simply the hypercharge of the final state 
fermions. We collect the branching ratios for B^> pair annihilations in Tab. [2j in the first 
column, Aj = 0, we assume a completely degenerate KK mass spectrum, while in the second 
one, A (1) = 0.14 the relative mass splitting between the LKP and KK quarks is assumed to 
be 0.14. Clearly the relative annihilation fractions are not particularly sensitive to the details 
of the spectrum so long as KK leptons are the same mass or lighter than KK quarks as the 
one-loop radiative corrections suggest. 

When the cross section in Eq. (f38|) is inserted into Eqs. ([36j) and ([37]) . a relic abundance 
within the range of the cold dark matter density measured by WMAP (0.095 < Qh 2 < 0.129) 
[72] can be obtained for m B (i) approximately in the range of 850 to 900 GeV. Lighter values 
lead to the LKP over-annihilating and thus under-producing the abundance of dark matter. 
Heavier values generate more dark matter than is observed (see the solid line in Fig. [8|). 
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process 


annihilation fraction 
A, = A qW = 0.14 




-> "eVe, V^V^, V T V T 


0.012 


0.014 




— ► e + e~, t + t~ 


0.20 


0.23 




— > uu, cc, tt 


0.11 


0.077 




— > dd, ss, bb 


0.007 


0.005 




- 


0.023 


0.027 



Table 2: The relative annihilation fraction into various final states. The numbers shown are 
not summed over generations, and the Higgs mass was assumed to be lighter than m B a) /2. 

This conclusion can be substantially modified if other KK modes have masses quasi-degenerate 
to the LKP and thus contribute to the freeze-out process [73]. This is expected to be especially 
important in UED models due to the quasi-degenerate nature of the KK spectrum, leading to 
many particle species being present at the time of freeze-out. To include the effects of other 
KK modes in the freeze-out process, we adopt the formalism outlined in Ref. [74] to account for 
co-annihilation processes in the process of thermal freeze-out of species. In Eqs. ([36]) and ([57]) . 
we replace the cross section (0, denoting the appropriate combinations of a and b) with an 
effective quantity which accounts for all particle species involved: 

a eS = ^ (1 + A,) 3 / 2 (1 + Aj) 3 / 2 e -(Ai+A,). (39) 

Similarly, we replace the number of degrees of freedom, g, with the effective quantity: 

3eff = ;>>(! + A *) 3/2e_lAl - ( 4 0) 

i 

In these expressions, the sums are over KK species, Uij denotes the coannihilation cross section 
between species % and j, and, finally, the As denote the fractional mass splitting between that 
state and the LKP. 

To illustrate how the presence of multiple KK species can affect the freeze-out process, we 
describe below two illustrative cases. First, consider a case in which the coannihilation cross 
section between the two species, 01,2, is large compared to the LKP's self-annihilation cross 
section, 01.1. If the second state is not much heavier than the LKP (A2 is small), then a e s 
may be considerably larger than 01 1, and thus the residual relic density of the LKP will be 
reduced. Physically, this case represents a second particle species depleting the WIMP's density 
through coannihilations. This effect is often found in the case of supersymmetry models in which 
coannihilations between the lightest neutralino and another superpartner, such as a chargino, 
stau, stop, gluino or heavier neutralino, can substantially reduce the abundance of neutralino 
dark matter [75] . 

The second illustrative case is quite different. If 0^2 is comparatively small, then the effective 
cross section for two species will tend toward e g- rj a\^g\j{g\ + 02) 2 + 02,252/(51 + 9%)' '■ ^ 
02,2 is not too large, a e g may be smaller than the LKP's self-annihilation cross section alone. 
Physically, this scenario corresponds to two species freezing out quasi-independently, followed 
by the heavier species decaying into the LKP, thus enhancing its relic density. An example of 
this behavior was found in Ref. [25] , in which the impact of KK leptons only slightly heavier 
than the LKP were considered in the freeze-out process. The resulting enhancement to the relic 
abundance of B^s is shown in Fig. [8] Notice that a similar situation can occur in the context 
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R" 1 (GeV) 



Figure 8: Left: The thermal relic abundance of £?Ws without the effects from any other KK 
species (solid line) and including the effects of KK leptons 5% and 1% heavier than the LKP 
(dashed and dotted lines). Shown as a horizontal band is the measured dark matter abundance 
from the WMAP experiment [72]. Adapted from Ref. [25]. Right: The change in the cosmo- 
logically preferred value for as a result of varying the different KK masses away from their 
nominal values in the minimal UED scenario (indicated by the circles on each line). Adapted 
from Ref. [73] . 



of supersymmetry, for instance with slepton coannihilations when the neutralino is higgsino- or 
wino-like, see e.g. Ref. [76] . 

Recently, the LKP freeze-out calculation, including all coannihilation channels, has been 
performed by two independent groups [771 I73j . We will briefly summarize their conclusions 
here. 

As expected, the effects of coannihilations on the LKP relic abundance depend critically on 
the KK spectrum considered. If strongly interacting KK states are less than roughly ~ 10% 
heavier than the LKP, the effective LKP annihilation cross section can be considerably enhanced, 
thus reducing the relic abundance. KK quarks which are between 5% and 1% heavier than the 
LKP lead to a with the measured dark matter abundance over the range of masses of 
m BW « 1500 to 2000 GeV (instead of 850-900 GeV for the case of the B^ alone). If KK 
gluons are also present with similar masses, m B (i) as heavy as 2100 to 2700 GeV is required to 
generate the observed relic abundance. 

On the other hand, if the strongly interacting KK modes are considerably heavier than 
the LKP, other KK states may still affect the LKP's relic abundance. If, for example, all 
three families of KK leptons are each 1% more massive than the LKP, then the observed relic 
abundance is generated for m B (i) between approximately 550 and 650 GeV. If the KK leptons 
are instead 5% more massive than the LKP, the observed abundance is found for m B (i) ~ 670 
to 730 GeV (see Fig. El- 

The discussion above is quantitatively illustrated in the right panel of fig. The plane indi- 
cates the value of the compactification radius R^ 1 and Aj, the relative mass splitting between 
the flW LKP and a given class of coannihilating KK particles: one (red dotted) or three (red 
solid) generations of SU (2)w-s'mglet KK leptons; three generations of 5f7(2)vK-doublet leptons 
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l/R (GeV) 

Figure 9: Regions, in the (R~ 1 ,mf l ) plane, where the relic abundance is consistent with 
the CDM abundance as inferred from the WMAP results, at the 1 and 2-a level (dark and light 
grey, respectively). The relic abundance computation takes here into account the effect of the 
second KK level resonances. Adapted from Ref. [78] . 

(magenta); three generations of SU{2) ^-singlet quarks (blue) (the result for three generations 
of .Scf (2)vy-doublet quarks is almost identical); KK gluons (cyan) and electroweak KK gauge 
bosons (green). Each line indicates the the value of R^ 1 needed for a given mass splitting Aj 
to get a B^ relic abundance £l B (t)h 2 = 0.1. 

Additionally, it was pointed out in Refs. [7_9l [80] that second level KK modes can also play 
a significant role in the freeze-out process, in particular in the case that va h (i) is very close 
to 2m B (i) . For ttikk =1 TeV, the cross section given in Eq. f|38|) can underestimate the full 
value by up to ~15%, thus increasing the mass range which generates the observed dark matter 
abundance. If (m h ( 2 ) — 2m B (i))/2m B (i) is larger than a few percent, however, second KK level 
particles have a negligible effect on the relic abundance calculation. 

Finally, assuming a minimal UED model where the matching conditions from brane-localized 
terms cancel at the cutoff scale, A, barring second-level KK modes resonances and setting the 
Higgs mass = 120 GeV and the cutoff scale AR = 20, Ref. [73 J finds that the preferred range 
for the compactification scale (and the LKP mass) is 500 < R-i/GeV < 600. In this range, 
the full computation of the thermal relic density of BW's, including all coannihilating particles, 
falls in the measured range of the cold dark matter abundance. Ref. [81] showed, again for a 
minimal UED setup, that for large Higgs masses, as the mass splitting between the B^ and 
the KK Higgses becomes small, coannihilation effects shift the viable region of compactification 
radii to larger values. For instance, at = 230 GeV, Ref. [81] quotes the range 800 < 
R^/GeY < 1200. Including the effect of the second KK level resonances, taking into account 
all coannihilation processes and varying the Higgs mass (but fixing AR = 20), according to the 
analysis of Ref. [78], the preferred range of R~ l for the minimal UED model consistent with 
the observed cold dark matter abundance is 600 < R^ 1 /GeV < 1400, as illustrated in Fig. [9] 

3.2 KKDM Protohalos 

While freeze-out signals the departure of WIMPs from chemical equilibrium, it does not signal 
the end of WIMP interactions in the early universe. Elastic and inelastic scattering processes 
of the form Xf — > Xf or Xf — > X' f keep the dark matter particle, X, in kinetic equilibrium 
until later times (lower temperatures) |82|, 183"], 184] . Here / and /' are SM particles in the 
thermal bath (leptons, quarks, gauge bosons) and X' is an unstable particle that carries the 
same conserved quantum number as X. The temperature of kinetic decoupling, 2\d, sets the 
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distance scale at which linear density perturbations in the dark-matter distribution get washed 
out — the small-scale cutoff in the matter power spectrum. In turn, this small-scale cutoff sets 
the mass, M c , of the smallest protohalos that form when these very small-scales go nonlinear at 
a redshift z ~ 70. There may be implications of this small-scale cutoff for direct |85j and indirect 
[86] dark-matter detection (see also Ref. [87]), although this depends critically upon the degree 
at which these structures are tidally disrupted through encounters with stars [881 [89"1 [87] . For 
instance, the dumpiness of dark matter on small scales enhances the rates for indirect detection 
by a factor that can be approximately cast as / • 5 [90] ; where / is the fraction of dark matter 
mass in substructures, and 5 is the dark matter density contrast, 



At present, however, no clear-cut, first-principle predictions are available for the size of the boost 
factor resulting from a given small-scale cutoff in the dark matter power spectrum. Sufficiently 
large substructures can be in principle also directly detected with gamma-ray telescopes, as we 
discuss in Sec. 13.41 

Ref. [91] showed that the WIMP's kinetic decoupling temperatures can vary over a large 
range, from several MeV to a few GeV. According to the estimate of Ref. [Q2] , 



which leads to a range of masses for the smallest protohalos (in generic WIMP models thermally 
producing the right amount of relic cold WIMPs) from 10 -6 M ffl to 10 2 M^ (M® denotes the 
mass of the Earth). Eq. (|42p accounts for both the acoustic oscillations imprinted on the power 
spectrum by the coupling between the dark matter and the relativistic particles in the primordial 
plasma prior to kinetic decoupling and the cutoff due to free-streaming of dark matter after 
kinetic decoupling. In particular, Ref. [91] addressed, together with supersymmetric models, 
the case of KK dark matter within the context of UED models. 

The kinetic decoupling temperature Tkd is defined by T r {T^) = i/ _1 (Tkd) [84], where H(T) 
is the Hubble expansion rate, and the relaxation time r r is given by 



Here, ni(T, mi) ~ T 3 is the equilibrium number density of the relativistic particle species I (the 
true mass dependence can be crucial here for some of the species under consideration such as 
the fjL and r leptons), aix(T) is the thermally averaged scattering cross section of the WIMP X 
off I's, and the factor (raj/T)" 1 counts the number of scatters needed to keep the WIMPs in 
kinetic equilibrium. Here, we consider I £ {v e ,n,r, e, fx, r, u, d, s, c} and neglect the scattering 
off light quarks below Aqcd — 150 MeV (our results are insensitive to the precise value taken 
for Aqcd and the detailed nature of the QCD phase transition). We neglect the scattering of 
WIMPs off mesons and baryons below Aqcd because this process is suppressed with respect to 
their scattering off light leptons by the relative abundance of the species in the thermal bath 
and by hadronic form factors. 

Within the context of UED, the scattering cross section of off light fermions proceeds 
through the Feynman diagrams shown in Fig. [Tol a). Unlike the case of neutralino dark mat- 
ter [83], the intrinsically degenerate nature of the KK spectrum, where m B (i) « m L (i), clearly 
enforces a resonant enhancement. Ref. [91] finds, to leading order in Ei/mx, and in the rela- 
tivistic limit for I and non-relativistic limit for the LKP particle, that 



J d 3 r d {p c i{rZi)) 2 
Po/d 3 r cl /Oci(rci)' 



(41) 



M c ~ 30 (r k d/10 MeV)" 3 M, 



(42) 




(43) 




(44) 
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(a.l) (a.2) (b) 



Figure 10: Feynman diagrams contributing to the scattering of B^ (a.l and a.2) and of z/W 
(b) off leptons. 




Figure 11: Left: The kinetic-decoupling temperature, T^, as a function of the LKP mass 
for minimal UED models featuring a B^ (with 2 < Ai? < 200, grey shaded region, and for 
AR = 20, solid black line) and a iA 1 ' LKP. Right: The WIMP protohalo characteristic comoving 
wavenumber, k c , (left axis) and mass, M c , (right axis) as a function of the LKP mass. 



where Yrl denote the hypercharge quantum number. In analogy with the case of neutralino 
dark matter, Ref. [91] pointed out that the oxi oc Ef scaling [83] holds for a B^ KK dark 
matter WIMP as well. The case of a KK neutrino was also studied in Ref. [SI], even though 
the KK neutrino is strongly disfavored as a dark matter particle by direct detection searches. 
KK neutrinos belong to a category of dark matter candidates which feature the same couplings 
to the Z gauge bosons as ordinary heavy neutrinos. In this case, one gets 

where gn,L stand for the couplings of the left and right handed lepton, I, to the gauge boson, 
Z, and g v (i) u (i) z = e/(sin 20yy). Again, the axi oc Ef scaling holds, and, remarkably, o v (i) [ does 
not depend on the LKP mass. 

In Fig. [11] we show the results for (left) and for M c (right), respectively, versus the 
LKP mass. We set the KK spectrum according to the minimal UED prescription for radiative 
corrections to the KK masses [30] (see Sec. I2.2.3p . setting the cutoff scale 2 < Ai? < 200, and 
showing the Ai? = 20 case with a black solid line. 

We conclude that KK dark matter protohalos feature a kinetic decoupling temperature lying 
between 10 MeV and a few hundred MeV. The resulting size of the smallest KK dark matter 
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protohalos varies between 10 -3 and 10 3 Earth masses. Recently, Ref. |93] also addressed the 
computation of kinetic decoupling of WIMPs in the early universe, including the case of the 
LKP (in particular, see their App. C for a complete derivation of the LKP scattering amplitude 
off light SM fermions). In the relevant range of masses, the resulting estimates for the kinetic 
decoupling temperatures agree with those of Ref. [91] . 



3.3 Direct KKDM Searches 

A large number of experiments have been designed and developed in the hope of observing 
the elastic scattering of dark matter particles with nuclei. Some of these experiments include 
CDMS [681 El], Edelweiss [95], ZEPLIN [96], CRESST [97], CLEAN, COUPP, DEAP, DRIFT, 
EURECA, SIGN, XENON [98], WARP [99], KIAS, NalAD [TOO], Picasso [TOT] . Majorana [102] . 
DUSEL, IGEX [103], ROSEBUD pE], ANAIS, KIMS, Genius [105], DAMA [106] and LIBRA. 
This collection of experimental programs is collectively known as direct detection. 



3.3.1 Spin-Independent Scattering 

The elastic scattering of a WIMP with nuclei can be separated into spin-independent and spin- 
dependent contributions. Spin-independent scattering can take place coherently with all of 
the nucleons in a nucleus, leading to a cross section proportional to the square of the nuclei 
mass. As a result, the current direct detection constraints on spin-independent scattering are 
considerably stronger than for the spin-dependent component. 

The spin-independent .B^-nuclei elastic scattering cross section is given by: 

^MiV.SI « B(1 \ N ,2 WP + ( A ~ Z )fn\\ (46) 

n{m B (V) +m N ) z 

where m^, Z and A are the mass, atomic number and atomic mass of a target nucleus. f p and 
/„ are the BW's couplings to protons and neutrons, given by: 



q=u,d,s H q=c,b,t 



where a q are the tfW-quark couplings and f^j « 0.020 ± 0.004, 0.026 ± 0.005, fff ~ 

0.118 ± 0.062, f^> ps 0.014 ± 0.003, f^> ps 0.036 ± 0.008 and fff ps 0.118 ± 0.062 [1071 
11081 |109j. The first term in Eq. (|47|) corresponds to interactions with the quarks in a target 
nucleon, whereas the second term corresponds to interactions with gluons in the target through 
a quark/KK-quark loop diagram, is given by 1 — — fjp — ffj ps 0.84, and analogously, 

f$, « 0.83. 

The coupling a q receives contributions from two classes of diagrams: the s-channel exchange 
of KK quarks and t-channel Higgs boson exchange (see Fig. fl~2]) . This leads to [26j [32]: 

,,2 ( V 2 >v2\(^2 , ^2 



m q g{ (Y q 2 R + Y q \) (m| (1) + m* (1) ) mq g2 



> ■-) ■> n o i (48) 

Hd)K(i)- m ,(i)) 8m BW m z h 

The first term in this expression should only be included for the light quarks (q = u, d, s), while 
the second term contributes to both light and heavy species. This leads to a contribution to 
the cross section from Higgs exchange which is proportional to l/(m^ (1) m 4 .) and a contribution 
from KK-quark exchange which is approximately proportional to l/(m^ (1) A 4 ), where A = 
( m g (i) — m Bw)/ m BW- Numerically, this leads to a .B^-nucleon cross section approximately 
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Figure 12: Leading Feynman diagrams for .B^-quark elastic scattering. IJ^-gluon scattering 
can also occur through analogous diagrams with a quark/KK-quark loop. From Ref. [32J. 
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Figure 13: The spin-independent .B^-nudeon elastic scattering cross section. In the left frame, 
rrih has been fixed to 120 GeV. In the right frame, m B (i) has been fixed to 1 TeV. The quantity 
A is defined as A = (m qW — m BW )/m B (i). From Ref. 



given by: 

1n , /lTeV\ 2 r/lOOGeVV / 1 TeV\ 2 /0.1 \ 2 1 2 

In Fig. [TBI this cross section is plotted as a function of m B (i) and m^. 

The range of elastic scattering cross sections predicted in this model are well below the 
current experimental limits. The strongest constraint on spin-independent scattering has been 
reported by the CDMS collaboration [68J, which finds roughly asi > 10~ 6 pb in the mass range 
of interest here. A number of other experiments have placed limits which are only mildly less 
restrictive, including ZEPLIN, Edelweiss and CRESST. 

CDMS and other direct detection experiments are currently working toward improving their 
sensitivity. Although the current constraints are around three orders of magnitude or more 
below the cross sections predicted for KKDM, the advanced phases of Super-CDMS should 
be able to reach the needed sensitivity to directly observe dark matter in these models. In 
addition to Super-CDMS, ton-scale experiments using liquid elements such as Argon or Xenon 
appear very promising. Although such noble liquid detectors may very well reach the level of 
sensitivity needed to observe KKDM, these technologies are currently less well understood than 
the cryogenic Germanium and Silicon detectors used by CDMS and Edelweiss. 

We illustrate the prospects for direct KKDM detection in future direct detection experi- 
ments in Figs. and [151 focusing on three minimal UED benchmark models, and exploring 
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Figure 14: Contours of constant spin-independent .B^-proton scattering cross sections in the 
(1/R, nth) plane, for two choices of AR = 4 and 20. The reach of the future direct detection 
experiments "Xenon-1 ton" and "Super-CDMS C" approximately correspond to the black solid 
line and to the red dashed line. We also indicate the location of three of the benchmark models 
of Appendix |Aj 



various slices of the minimal UED parameter space (for a detailed description of the motiva- 
tions and of the particle spectra for these models see Appendix |A|) . For all benchmark models, 
a vanishing Higgs boundary mass term Model is assumed, and the models are therefore fully 
defined (through the radiative mass corrections outlined in Sec. 12.2.3]) by the values of 1/R, AR 
and m/j. All benchmark models feature a LKP with a thermal relic abundance in accord 
with the cold dark matter content of the universe. Model UED1 has 1/R = 550, AR = 20 and 
m h = 120 GeV, model UED2 has 1/R = 850, AR = 4 and m h = 120 GeV, and model UED3 
1/R = 1000, AR = 20 and m h = 220 GeV. 

Specifically, we explore the .B^-proton scattering cross section in Fig. [T4"l where we scan 
the (1/R, mh) plane, for two representative values of AR = 4, 20 (which determine the mass 
splitting, A). The grey shaded regions are excluded by the negative results of the LEP2 Higgs 
searches. We also plot curves representing the reach of representative future direct dark matter 
search experiments (the projected sensitivity at a WIMP mass of 1 TeV is that indicated on each 
curve) . The black diamonds locate the benchmark models of Appendix [A] The approximate 
sensitivity reach of the future Xenon-1 ton and Super-CDMS C experiments approximately cor- 
respond to the solid black and to the dashed red curves, respectively. While models UED1 and 
UED2 will both escape detection at Xenon-It, they lie within reach of the planned sensitivity 
for the phase C of Super-CDMS. Model UED3, on the other hand, will not be detectable with 
foreseeable direct detection devices. Fig. [15] shows the same curves, and two of the benchmark 
models, this time in the (AR,mf r ) plane. 

These two figures demonstrate how the three parameters defining the minimal UED setup 
affect direct detection rates. We notice that in the low-l/i? regime the role of mh is less critical 
than with a heavier KK spectrum. Also, the parameter AR plays a role (affecting the LKP-KK 
quark spectrum) only when it takes low values < 10; at larger AR, the direct detection rates 
become largely independent of AR. 
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Figure 15: As in Fig. [T4"l but in the (AR, m/J plane. 



3.3.2 Spin-Dependent Scattering 

In addition to spin-independent elastic scattering, axial- vector couplings can lead to .B^-nuclei 
spin-dependent scattering. As we mentioned earlier, the current experimental constraints on 
spin-dependent scattering are far weaker than for spin-independent processes. New experimental 
technologies, however, such as those being developed by the COUPP collaboration, may enable 
greater sensitivity to spin-dependent couplings. Spin-dependent scattering can also play an 
important role in the capture of WIMPs in the Sun (see Section [375]) . 
The cross section for spin-dependent scattering is given by: 

Hm 2 Bil) m 2 N A 2 J(J + 1) 
^ (1) ^ SD ~ Mm BW +m N Y (ml (1) -ml w r ^ 

where J is the spin of the nuclei species and the quantity A is given by: 

^ — a p(S~p) + O-n(Sn) 



where 
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Here the Ags parameterize the fraction of spin carried in a proton by each quark species: 
Au = 0.78 ± 0.02, Ad = -0.48 ± 0.02, and As = -0.15 ± 0.07 [TTUl fTTT] . 
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For spin-dependent scattering with a proton, the cross section is well approximated by: 
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And with a neutron: 
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The strongest experimental limits on spin-dependent scattering in the mass range we are 
interested in have been set by the CDMS [94J and NalAD [100J experiments, for scattering 
with neutrons and protons, respectively. These experiments only constrain <tsd ^ 0.1 — lpb, 
however, and thus do not have the sensitivity needed to test KKDM models. 

3.4 Gamma Ray Searches 

In addition to detecting WIMPs directly, efforts are underway to observe the stable products 
of dark matter annihilations, including gamma rays, energetic neutrinos and antimatter. The 
relevant particle model input quantities for indirect detection are the pair annihilation cross 
section and the relative frequency (branching ratio) of a given pair- annihilation final state. We 
outlined both in Sec. 13.11 and in Tab. [2] above. 

In the case of KKDM, gamma rays can be produced as continuum photons from final state 
radiation and the cascades of other annihilation products, and as line emission from loop- 
diagrams to 77, jZ or 7/1 final states. 

The locations most likely to produce an observable flux of gamma rays from dark matter 
annihilation are those which contain a very high density of dark matter and are relatively 
nearby. The center of our Galaxy has long been considered one of the most promising regions 
of the sky for observing dark matter annihilation radiation [112} 11131 II 14] . 

The flux of gamma rays from the Galactic center is found by integrating the density squared 
along the line-of-sight to the observer [115] : 

$ 7 (V,£ 7 ) = av <i Z\ \ [ dsp 2 (r(s,ij)), (55) 

dE 1 47rm^ (1) ./line-of-sight 

where av is the WIMP's annihilation cross section, dN^/dE~ f is the gamma ray spectrum pro- 
duced per annihilation, and the coordinate s runs along the line-of-sight, in a direction making 
an angle, ift, from the direction of the Galactic center. p(r) is the density of dark matter at a 
distance r from the Galactic center. 

This expression can be separated into factors which depend on the dark matter distribution, 
and which depend on the characteristics of the dark matter species. The first of these is described 
by the quantity J(tp): 

J M = 5-^— f n [ dsp 2 (r(s^)). (56) 

8.5 kpc V0.3GeV/cmV Jlme-of-sight 



We further define J(AO) as the average of J(tp) over the solid angle, AQ, centered on 1(1 = 0. 
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P 


7 


R (kpc) 


J(l0- 3 str) 


NFW 


1.0 


3.0 


1.0 


20 


1.352 x 10 3 


Moore 


1.5 


3.0 


1.5 


28.0 


1.544 x 10 5 


Iso 


2.0 


2.0 





3.5 


2.868 x 10 1 


Kra 


2.0 


3.0 


0.4 


10.0 


2.166 x 10 1 



Table 3: Parameters, and values of J averaged over 10 3 steradians, for some commonly used 
galactic halo profiles, as used in Eq. (|58p . The models shown are those of Navarro, Frenk and 
White (NFW) [116] . Moore et al. (Moore) [TTT] . modified isothermal (Iso) [TT5] and Kravtsov 
et al. (Kra) [Tig] . 

We can now express the flux of gamma rays over a solid angle AO, around the Galactic 
center: 

« 2.8 , 10- ^ ( 3xl0 - cm3/8 ) (^) 2 7(A«) AO c m -V (57) 



There are considerable uncertainties involved in the distribution of dark matter in the Galactic 
center region. N-body simulations suggest that there exists a universal dark matter profile, with 
the same shape for all masses, epochs and input power spectra [116J. Such a profile is often 
parameterized as 

PM = (r/Rp[l + {r/RY}^)/ a ' (58) 

In Table 13.4} the parameters used in the parameterization of Eq. (|58p are given for several 
representative halo profiles. Of these profiles, NFW [116] and Moore et al. [117] each predict 
very dense concentrations of dark matter in the inner regions of the Galactic halo (p(r) oc r _1 
or r -1 ' 5 ). Such a cusped halo profile is typically required if gamma rays from dark matter 
annihilations are to be observed from the Galactic center. 

None of above profiles derived from N-body simulations account for the effects of baryonic 
matter, however. As the baryonic density dominates over that of dark matter in the inner 
regions of the Milky Way, it is possible that N-body simulations will not accurately predict the 
dark matter distribution in this region. The effects of baryons on the dark matter profile are, 
unfortunately, not easily predicted. One possibility is that the effects of adiabatic compression 
on the dark matter profile near the Galactic center could play an important role, potentially 
enhancing the dark matter density by up to an order of magnitude in the inner parsecs of 
the Milky Way [119J. Baryonic matter could also flatten the inner cusp of the dark matter 
distribution. 

At the center of the Milky Way is a 2.6 x 10 6 solar mass black hole. It is possible that the 
adiabatic accretion of dark matter onto this object could lead to extremely high densities of 
dark matter [120J. An initial density profile with a slope, p(r) oc r -7 , would lead to a density 
"spike" with a slope of 7 sp = (9 — 2j) /(4 — 7) [121 j . Such a feature would dramatically increase 
the rate of dark matter annihilations near the Galactic center. 

A spectrum of gamma rays extending from at least 200 GeV to 10 TeV (and above) has 
recently been observed from the Galactic center by the HESS [122] and MAGIC [123] Atmo- 
spheric Cerenkov Telescopes (ACTs). While initially it was thought to be possible that this 
signal was the result of dark matter annihilations |124|, 1125} 1126} 1127] . the spectrum now ap- 
pears incompatible with that predicted from dark matter and is likely to be of an astrophysical 
nature. This spectrum represents a challenging background to future searches for dark matter 
in the Galactic center pL28J. 
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Figure 16: Contributions to B^sW £+^-7. From Ref. [127] . 



With this potential obstacle involving searches for dark matter near the Galactic center in 
mind, other regions are beginning to appear more attractive. In particular, a class of companion 
galaxies to the Milky Way, called dwarf spheroidal galaxies, may provide an observable flux of 
dark matter annihilation radiation from low background regions of the sky [1291 11301 H31| I132j . 
Galaxies external to the Milky Way have also been considered |133j . Alternatively, if density 
spikes form around Galactic intermediate mass black holes, a number of point sources of dark 
matter annihilation radiation may become observable in the near future [134} 1135] . 

3.4.1 The Gamma Ray Spectrum 

The majority of KKDM annihilations produce pairs of quarks or leptons. Many of these particles 
then proceed to decay and fragment, producing a continuous spectrum of gamma rays [136]. At 
modest and high energies, this spectrum is dominated by photons produced in the semi-leptonic 
decays of tau leptons and by radiative processes (bremsstrahlung) associated to and e + e~ 

final states. All these charged lepton pairs final states are produced in approximately 20% of 
B^ annihilations |127j (see Tab. [2]). This is clearly in contrast with the case of the lightest 
neutralino, where, as explained in Sec. 13. 1\ leptonic final states are suppressed by factors of the 
order m^/m 2 , and the high energy end of the gamma-ray spectrum is typically dominated by 
gauge bosons pair final states, when these are open. The latter, here, do not, instead, contribute 
appreciably to shape the gamma-ray spectrum. 

Final state radiation from the process '7 (see Fig. [To]) plays, in the context of 

KKDM, a particularly important role. The spectrum of photons produced through this process 
(for a given lepton species) is given by: 



where x = E^/m B (i). Final state radiation is very important in determining the gamma ray 
spectrum from KKDM largely because of the mentioned large fraction of annihilations that 
generate light lepton pairs (approximately 20% to each of e + e~ and /i + /x~). In Fig. [T71 the 
contributions from quark fragmentation, tau decays, and final state radiation are compared. At 
the highest energies, final state radiation dominates. At more modest energies, semi-leptonic 
tau decays dominate. At low energies, most of the gamma rays are produced through quark 
fragmentation. 

In addition to continuum gamma rays, the processes B^'B^ 1 ' —* 77, B^B^ — > jZ 
and B^B^ — > 7/1 result in mono-energetic gamma ray lines with energies of E~ = m B (i), 
m B (!) (1 — m 2 z / Am 2 BW ) and m B (i) (1 — m|/4m^ (1) ), respectively. Such gamma ray lines would 
be a "smoking gun" for dark matter annihilations if they were to be observed. 

The cross section for B^B^ — > 77 involves contributions from twelve Feynman diagrams 
containing a fermion/KK-fermion loop and 22 diagrams involving the scalar sector of the the- 
ory [137J. The fermion diagrams dominate the cross section, which was found to be approxi- 
mated by cr B (i) B (i) v ~ 2 x 10~ 30 cm 3 /s (1 TeV/m B (i)) 2 . The cross section also grows as the 




dN 7 _ d(ai+i- 7 v)/dx _ a (x 2 - 2x + 2) \ m B m 



(59) 
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Figure 17: The spectrum of photons per B^B^ annihilation (solid line). Shown as a dashed 
line is the contribution from quark fragmentation alone. The dotted line is the quark fragmen- 
tation contribution plus the contribution from r leptons. From Ref. [127J. 

mass of the KK-fermions is moved closer to the LKP mass. Varying the (?7if(i) — m B p)) /itt-bW 
between 0.3 and 0.01 increases the cross section by a factor of roughly 3 [137j . 

Calculations of the cross sections of B^B^ —* 7Z and B^B^ — > 7/1 have not yet 
appeared in the literature, although preliminary work indicates that they are not expected to 
be as favorable channels as B^B^ — > 77 for the monochromatic gamma ray line emission 
[T38] . 

3.4.2 Observational Prospects 

Both ground and space-based gamma ray detectors are potentially sensitive to dark matter 
annihilation radiation. A number of Atmospheric Cerenkov Telescopes (ACTs), including HESS 
|139j . MAGIC |140j and VERITAS [141 j are currently operational. These experiments are 
sensitive to gamma rays above 50-500 GeV, depending on the telescope and the location of the 
source in the sky (its zenith angle). Additionally, the GLAST satellite is scheduled to begin its 
mission in 2007 [142, 143J. GLAST will be sensitive to gamma rays between sub-GeV energies 
and about 300 GeV. Unlike ACTs, GLAST will be able to observe large regions of the sky at 
a time, although with a much smaller effective area (approximately one square meter) than 
ACTs. 

The prospects for the detection of gamma rays from KKDM annihilations depend critically 
on the distribution of dark matter in the center of our Galaxy, in dwarf spheroidals, in dark 
matter clumps and in other such regions. In most cases, these distributions are very poorly 
known, making the prospects for future gamma ray dark matter searches difficult to predict. 

For the case of the Galactic center, considering an NFW halo profile and m B (i) ~ 1 TeV, we 
expect GLAST to detect ~ 10 photons above a few GeV over its first three years (assuming the 
Galactic center is in its field-of-view about one third of the time). This rate could potentially 
represent a detection of dark matter if it were not for the backgrounds which have been observed 
by HESS and MAGIC |128j . A rate roughly an order of magnitude larger would be needed to 
overcome these backgrounds, assuming that the power-law spectrum observed by HESS and 
MAGIC extends to lower energies [128] . A lighter dark matter particle would be somewhat less 
difficult to detect with GLAST. 

The annihilation rate and corresponding gamma ray flux from dwarf spheroidals is gener- 
ally expected to be somewhat lower than is predicted from cusped halo models of the Galactic 
center, although predictions can vary considerably. The prospects for detecting dark matter 
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annihilation radiation from such objects is difficult to assess at this time. The prospects for 
the detection of a signal from dark matter mini-halos or clumps (see Sec. I3.2|) strongly depends 
on the survival probability of such objects and on their specific dark matter distribution. Po- 
tentially GLAST might be able to detect gamma-rays from annihilations of KKDM in such 
structures |144j . or even the proper motion of the sources themselves [145J. 



3.5 Neutrino Searches 

As WIMPs annihilate in the halo, they do not only produce gamma rays, but also other particles 
including neutrinos. Neutrinos are far more difficult to detect than gamma rays, however. For 
this reason, another approach must be taken if neutrinos from dark matter annihilations are 
to be observed. In particular, WIMPs can be captured in the Sun, where they annihilate 
efficiently. Of the various annihilation products, only neutrinos can escape the solar medium 
and potentially be observed. 

WIMPs become captured in the gravitational potential of the Sun at a rate given by [146] : 



is -i f Piocal \ ( 270km/s\ 3 / oh,sd + 2.6<th,si + 0.175 a H e,sA / lTeV\ 



2 



C « 3.35x10 s-\ 03GeV/cm3 ){— {— ~ ) ) 

(60) 

where pi oca ,\ is the local dark matter density and V\ oca \ is the local RMS velocity of halo dark 
matter particles. cth.sDj ch,si and ctie.si are the spin-dependent, .B^-proton (hydrogen), spin- 
independent, .B^-proton and spin- independent, .B^-helium elastic scattering cross sections, 
respectively. The factors of 2.6 and 0.175 include information on the solar abundances of 
elements, dynamical factors and form factor suppression. 

As shown in Sec. 13.31 B^'s elastically scatter with protons far more efficiently through spin- 
dependent than spin-independent couplings (note that this is overturn by coherence effects once 
one considers the scattering off nuclei rather than nucleons). The dominant elastic scattering 
cross section for capture in the Sun is thus cr B (i) p SD ~ 1-8 x 1CT 6 pb {m B (i) /TeV)~ 4 (A/0.1)~ 2 . 

The evolution of the number of B^s in the Sun, N, is given by: 

N = C Q - AqN 2 , (61) 

where Aq is the annihilation cross section multiplied by the relative WIMP velocity per unit 
volume, given by: 

Here, V e s is the effective volume of the Sun's core, which is determined by matching the core 
temperature to a WIMP's gravitational potential energy at the core radius. This was found to 
be [1471 [148] 

V cS = 1.8 x 10 26 cm 3 ( ±™) . (63) 



m BW J 

This leads to a present annihilation rate of B^'s in the Sun of: 

r = ^A Q N 2 = X - C Q tanh 2 [^C^t & ) (64) 

where t & ~ 4.5 billion years is the age of the solar system. If ^/Co-Aoi© S> 1, then the anni- 
hilation and capture rates reach equilibrium, maximizing the resulting neutrino flux. Capture- 
annihilation equilibrium is reached if the following condition is met: 

<* « v ( M A V2 ^H,SD + 2.6a H ,si + 0.175aHe,si A 1/2 / 1 TeV\ 1/4 

3 ' 35X UxlO- 26 cm3/s J { lO^vb J l^rj >>L (65) 
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Inserting the appropriate annihilation and elastic scattering cross sections, the left-hand side of 
this condition reduces to: 

leading us to conclude that capture-annihilation equilibrium will generally be reached in the 
case of KKDM. 

Annihilations of B^'s then proceed to generate neutrinos either directly or through the 
cascades of other annihilation products. The resulting muon neutrino spectrum at Earth is 
given by: 

dN v C Q F Eq fdN u \ lni (g7) 



dE v 4vrL»2 \ dE 



ES 

Here, F E(i = tanh 2 (^/ U^X^t®) is the non-equilibrium suppression factor (ss 1 for capture- 
annihilation equilibrium), Z?es is the Earth-Sun distance, and (^^) In ^ is the neutrino spectrum 
injected in the Sun per annihilating B^\ Most of the energetic neutrinos produced come from 
annihilations to either neutrino or tau pairs. For the 2% to 3% of annihilations which 
produce tau neutrino or muon neutrino pairs, this injected spectrum is simply a delta function 
at E v = m B (i) (neglecting the effects of neutrino scattering in the solar medium). More than 
20% of B^ annihilations produce tau pairs, which generate a hard spectrum of neutrinos 
through their leptonic decays r~ — > /i _ z/ T z^ and a softer component through charged pions 
from hadronic cascades. 

The spectrum of Eq. (|67p can be modified due to a number of factors. Firstly, energetic 
neutrinos can scatter with nucleons in the Sun. The probability of a neutrino escaping from 
the Sun without interacting is approximately given by [1491 1150] : 

P = ff-Zx/Bk (Q 8 ) 

where ~ 130, 160, 200, 230 GeV for u^, u T , V^, V T , respectively. Most of these interactions 
are through charged current, and therefore result in the absorption of the neutrino. Charged 
current interactions of tau neutrinos, however, generate a tau-lepton which can decay before 
losing the majority of its energy to the solar medium. This leads to an effect known as tau- 
regeneration. Finally, neutrino oscillations must also be taken into account. Muon and tau 
neutrinos mix into equal flavors as they travel from the Sun to the Earth. 

Upon reaching the Earth, a small fraction of the muon neutrinos interact via charged current 
in the ice or water near a neutrino telescope, generating energetic muon tracks at a rate of: 



dN Ufi do v 



Events =i / / -j^^iEu^y) - y)E v ) A eS dE u Jy, (69) 



where R^{E^) is the distance a muon travels below falling below the energy threshold of the 
detector, called the muon range, and A e g is the effective area of the detector. In Fig. [HI we 
show the spectrum of muons produced in this way. These muons lead to an event rate in a 
kilometer-scale neutrino telescope such as IceCube |151| 1152] or KM3 [153] , which is shown in 
Fig. [H 

For these neutrinos to be identified, they must overcome the background of atmospheric 
neutrinos. Above 100 GeV, this corresponds to roughly 80 background muons in the Sun's 
angular window per square kilometer, per year. Over a decade of observation, a 3a detection 
at a kilometer-scale experiment would, therefore, require a rate of ~ 3^80/10 ~ 8 per square 
kilometer, per year. 
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Figure 18: The spectrum of muons generated through charged current neutrino interactions at 
Earth, for the various annihilation modes of the B^K In the left and right frames, m B (i) = 600 
GeV and 1000 GeV, respectively. The spectrum is dominated by annihilations to neutrino and 
tau pairs. From Ref. [154J. 




Figure 19: The event rate in a kilometer-scale neutrino telescope resulting from B^ annihila- 
tions in the Sun. Values of 0.1, 0.2 and 0.3 for the quantity A = (m q (i) — Tn B {V))/m B (i) have 
been used. From Ref. |154j . 
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We see from Fig. [TUJthat observable event rates can potentially be generated for a combina- 
tion of a light B^ and small mass splittings between the LKP and KK-quarks. For KK-quarks 
10% heavier than the LKP, an observable rate is expected for m B (i) < 900 GeV, for example. 

Currently, constraints on the energetic neutrino flux from the Sun have been placed by 
the Super-Kamiokande pj>5], BAKSAN [T5S]. AMANDA [HU and MACRO [158] experiments. 
These experiments only limit the rate to less than ~2000 events per square kilometer, per year, 
however, so do not yet constrain these models. Future experiments with considerably larger 
effective areas will improve upon these limits considerably, potentially reaching the level of 
sensitivity needed to observe on the order of ten events per square kilometer, per year. 

Comparing the case of KKDM to neutralinos in supersymmetric models, we point out that 
larger neutrino rates are typically expected for KKDM for a number of reasons. Firstly, the flux 
of neutrinos produced is enhanced by the relatively large fraction of B^ annihilations which 
produce t + t~ (20-25%) and Vpp^ or v T v T (2-3%). Neutralinos, in contrast, annihilate largely 
to heavy quarks and gauge or Higgs bosons, which produce fewer energetic neutrinos in their 
decays. Secondly, B^'s elastically scatter in the Sun primarily through spin-dependent cou- 
plings, and thus are not limited by the results of CDMS and other direct detection experiments. 
The elastic scattering of neutralinos in the Sun can be dominated by either spin-dependent or 
spin-independent contributions. Spin-independent scatterings typically lead to either a direct 
detection rate in excess of the CDMS limit, or a neutrino rate below the reach of kilometer-scale 
neutrino telescopes [159J. 

3.6 Anti-Matter Searches 

The possibility of revealing the presence of an exotic particle population in our Galaxy through 
cosmic-ray searches has been considered for a long time [160J. It was soon realized that if the 
"missing matter" was made up of WIMPs, one could hope to indirectly detect it through the 
detection of its stable annihilation products. The relative under-abundance of antimatter in 
cosmic rays, and the fact that WIMP annihilations typically yield as much matter as anti- 
matter, motivated the possibility of disentangling an exotic signature originating from WIMP 
annihilations in the Galactic halo from the secondary and tertiary cosmic-ray background. 

The possibility of using positrons [161} 1162} 1163] and low-energy antiprotons [161j in this 
way has been discussed for some time. Although early studies mainly focused on the possibility 
of tracing anomalies in the observed antimatter spectra back to an exotic contribution from 
neutralino annihilations in the Galactic halo |163|, 11641 I165j . the possibility of constraining 
dark matter models (particularly neutralino dark matter scenarios) through their antimatter 
yields was also outlined |166l 1167] 1168] . Following the subsequent release of experimental data 
[169} 1170] , it became clear that a standard cosmic ray background could in most cases provide 
an excellent fit to the available antimatter fluxes [167] . What is sometimes dubbed "HEAT- 
anomaly" [171] , an excess in the high energy end of the positron fraction (the ratio of positrons 
to positrons plus electrons), is a possible exception, although the measured positron spectrum 
itself (as opposed to the positron fraction) can be satisfactorily fitted with a pure astrophysical 
background, and the associated statistics are rather poor. The possible interpretation of this 
excess in terms of positrons produced through Galactic dark matter annihilations can be found, 
for example, in Refs. \WE \ \T72 \ [173] . 

Antimatter searches are entering an exciting era, where space based experiments such as 
PAMELA and AMS-02 [174\ I175| will collect an unprecedented wealth of data covering a sig- 
nificantly larger energy range than has been explored thus far by balloon-borne experiments 
[1691 IT70j . Furthermore, these new experiments will help in constraining and appropriately mod- 
eling the cosmic-ray background as well. Recent analysis addressing the question of whether 
these next-generation experiments will have the ability to statistically discriminate a WIMP- 
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annihilation-induced positron and antiproton component for various particle physics setups and 
dark matter halo profiles can be found in Refs. |168} 1176} 1177] . 

The possibility of detecting a dark matter annihilation signature in the antideuteron flux was 
envisioned only much more recently [1781 H79j. Compared to the flux of antiprotons produced in 
WIMP annihilations, the number of antideuterons is significantly suppressed. Since the thresh- 
old energy for the production of an antideuteron in proton-proton collisions is rather large, 
however, the impinging proton must feature a large kinetic energy. Consequently, the center- 
of-mass velocity of the proton-proton system is typically large, and the resulting antideuteron 
kinetic energy sizable. Therefore, the background of low kinetic energy antideuterons is ex- 
pected to be highly suppressed. The exciting experimental possibility of selectively capturing 
slow antideuterons, and of identifying them with a remarkably high accuracy through atoms 
X-ray de-excitations, has lead to the design of a dedicated device, the general antimatter par- 
ticle spectrometer, GAPS [180] . In some setups, even the detection of a single low-energy 
antideuteron could constitute a signature of exotic physics, including that of dark matter anni- 
hilations in the Galactic halo |181j . 

In the following sub-section, we discuss the latter case of KKDM detection by low-energy 
antideuteron search experiments. The case of antiproton and positron searches are discussed in 
Sec. 13.6.21 and 13.6.31 respectively. In Sec. 13.6.41 we discuss the prospects for KKDM detection 
by space-based antimatter exploration experiments. 



3.6.1 Antideuterons 

Antideuteron searches can be performed either with magnetic spectrometers mounted on balloon- 
borne (BESS/BESS-Polar) or space-borne (AMS) missions, or through GAPS-like devices, 
based on the radiative emissions of antiparticles captured into exotic atoms. The latter can 
be installed either on balloons or on satellites, and are specifically designed to look for low- 
energy antiparticles. 

For both antideuterons and other antimatter species (see the following two sections) the 
effects of solar modulation can be accounted for within the framework of the Gleeson-Axford 
analytical force-field approximation [182J , where the interstellar flux at the heliospheric bound- 
ary, d^t/dTfe, and at the Earth, d^/dT^, are related by 

d$®, , _ Pi 
dT e pi dT b 

where the energy at the heliospheric boundary is given by E b = E^ + \Ze\4>p, and p b and p^ 
stand for the momenta at the boundary and at the Earth, and (f)p is the above mentioned solar 
modulation Fisk [183J parameter, which is taken to be approximately charge-sign independent. 
Physically, the Fisk parameter corresponds to the effective electric potential associated to the 
presence of solar wind. 

The flux of low-energy antideuterons is currently constrained by the recent results of Ref . |184j . 
The BESS experiment looked for low-energy antideuterons during four flights (1997, 1998, 1999, 
and 2000) in the kinetic energy interval 0.174- 1.15 GeV/n. The upper and lower kinetic energy 
limits come, respectively, from the particle identification procedure and from the decrease of 
the geometrical acceptance and mean free path through the detector. Without assumptions on 
the D spectral shape, the BESS collaboration, by combining all four missions, derived an upper 
limit on the D flux, at 95% C.L., of 

^BESS < L g x 1Q -4 m -2 s -l sr -l( GeV / n )-l. ( 71 ) 

The Fisk solar modulation parameter &p was derived from the p data from the same experiment, 
and is set to 500, 610, 648, 1334 MV for the four years of data-taking. 
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The future space-based spectrometer AMS-02 will be sensitive to low-energy antideuterons [185J. 
For a total data-taking time of 3 years, and for an antideuteron kinetic energy band extend- 
ing from the AMS threshold of 100 MeV/n to ~ 2.7 GeV/n, the inferred acceptance reads 
5.5 x 10 7 m 2 s sr GeV [175] . 

The GAPS experiment has recently undergone a rich phase of R&D, carried out at the KEK 
accelerator in Japan [180J. It has been realized that solid and liquid targets can greatly simplify 
the needed payload mass (thanks to the removal of the dead mass of the gas handling system) 
and the complexity of the apparatus, yielding an increased background rejection capability 
enabling the capture of more than 3 X-rays, as initially conceived. Furthermore, pion showers 
(it*) and nuclear X-rays from the antiparticle annihilation in the target nuclei, neglected in the 
original sensitivity calculations |179j . have been shown to significantly increase the antiparticle 
identification capability. The GAPS collaboration plans to test the finalized payload with a 
prototype as early as 2009, and to achieve a long duration balloon (LDB) flight from Antarctica, 
or an ultra-LDB (ULDB) flight from Australia as early as 2011 [180, 186J. A preliminary 
evaluation of the two balloon-borne options (with the LDB sensitivity based on 3 flights) found 
projected sensitivities of [186J 

<^ DB « 1.5 x 10~ 7 m -2 s- 1 sr~ 1 GeV~ 1 and c/>^ LDB ps 3.0 x 10~ 8 nrV^^GeV -1 (72) 

over a bandwidth of 0.1 < T-£j/(GeV/n) < 0.25. The sensitivity of the GAPS prototype mounted 
on a satellite was assessed in Ref. |179j . For a 3 year satellite mission, the projected sensitivity 
is 

^gaps/s _ 2 g x 1() _9 m -2 s -i sr -l GeV -l ; for 0.1 < %/(GeV/n) < 0.4. (73) 

A last, very optimistic option, mentioned in Ref. [179] , is to send GAPS on a deep-space 
probe, where eventually solar modulation effects can be significantly reduced. Depending on 
the spectral shape of the differential D yield from WIMP pair annihilation, solar modulation 
can deplete the low energy antideuteron flux, hence this interplanetary GAPS setup might 
represent the ultimate probe for dark matter searches via the detection of low energy Ds [177] . 

The computation of the differential flux of Ds per kinetic energy per nucleon interval induced 
by WIMP pair- annihilations involves a number of steps, which we briefly review below, referring 
the reader to Refs. [1781 11771 H81] for more details. 

The computation of the source spectrum for the primary antideuteron flux originating from 
WIMP pair annihilation is based on three hypothesis: (1) correlations in the associated produc- 
tion of two antinucleons are neglected, i.e. the probability of producing a pair of antinucleons 
is assumed to be equal to the product of the probabilities of producing two single antinucleons 
(factorization) (2) the antineutron production cross section is equal to the antiproton production 
cross section (isospin invariance) and (3) the formation of an antideuteron can be described by 
the coalescence model. We refer the reader to Refs. [1871 11781, 1188] 1189] for a through discussion 
of the validity of these assumptions. The factorization assumption is conservative, in that the 
probability of pair producing antinucleons in the same jet is presumably not factorized, since 
their momenta will not be isotropically distributed. The main idea of hypothesis (3) is that 
whenever the difference of the momenta of an antiproton and an antineutron produced in a 
jet resulting from a WIMP pair annihilation is less than a phenomenologically given value 2po, 
where po indicates the coalescence momentum, then an antideuteron is formed. 

The differential energy spectrum of primary antideuterons produced in the pair annihilation 
of KKDM can then be expressed by [178] 
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Figure 20: Left: The differential antideuteron flux from the four benchmark models (UED1-4) 
of Appendix[A](the results for UED2 and UED3 fall almost on top of each other). The shaded 
regions correspond to the sensitivities of various existing and proposed experiments featuring 
antideuteron searches. Right: The sensitivity reach of antideuteron search experiments for the 
UED model, with a B^ LKP (dashed black lines) in the (m B (i), (av)o/in 2 BW ) plane. We also 
indicate the location of three of the benchmark models of Appendix lAl 



where = + k^j, f indicates any final state of the WIMP pair annihilation process occur- 
ring with a branching ratio BR(B^ B^ 1 — > /), and dNp^ /dEp is the antiproton differential 
yield for the final state /. The latter can be computed using the Pythia Monte Carlo event 
generator [1901 1191] . as implemented in the DarkSUSY package |192j . for example. The reference 
value we assume for the coalescence momentum is 58 MeV, the same choice as in Refs. j!87j and 
[178) . not too far from what is expected from the antideuteron binding energy, \fm p B ~ 46 
MeV. For a discussion of the effects of the propagation of antideuterons through the Galactic 
magnetic fields, we refer the reader to Ref. [181J. 

The source spectrum of antideuterons is specified at every point in the Galactic halo once 
the shape of the dark matter halo itself is given. We take here as an (optimistic) reference 
model the adiabatic contraction [193 J of the N03 halo profile [194J (see Ref. [195J for details), 
which closely resembles the profile proposed by Moore et al., [196J. We consider here a smooth 
halo profile, but we discuss in Sec. 13.6.41 the potential effects of dark matter halo substructures. 

We collect in the left frame of Fig. [20] the results of the computation of the antideuteron 
flux at Earth for the four benchmark models found in Appendix [A[ We assume here a Fisk 
parameter <j)p = 800 MV. In the figure, the shaded regions correspond to the sensitivities 
and kinetic energy range of various existing and proposed experiments featuring antideuteron 
searches, as outlined above. Of the benchmark models we consider here, models UED1, UED2 
and UED3 are likely to generate at least one low-energy D at a space-based GAPS experiment 
(grey-shaded area). Model UED1 might even produce some events at a balloon-borne GAPS 
experiment or at AMS-02. 

In the case of highly sensitive space-based antideuteron search experiments, it has been 
pointed out in Ref. [181) that the role of the cosmic ray background cannot be neglected (unlike 
the case of, e.g. balloon-borne setups). Previously neglected antideuteron production and 
energy loss processes, including secondary antideuteron production from antiproton scattering 
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Figure 21: The local interstellar (left) and solar-modulated (right) antiproton flux, as a function 
of the antiproton kinetic energy, for the four benchmark models (UED1-4) of AppendixlAl (the 
results for UED2 and UED3 fall almost on top of each other). In the right frame, we also 
include an estimate of the cosmic ray antiproton background [177] and the current antiproton 
flux measurements [170] 



off the ISM, and a tertiary antideuteron component originating from the previously neglected 
non-annihilating inelastic scattering processes, have been shown to largely populate the low- 
energy end of the antideuteron spectrum }189j . As a result, the detection of a single low-energy 
antideuteron would not be sufficient to claim an exotic signal, as it would more likely originate 
from the standard cosmic ray background |181j . A statistical approach showed that the number 
of low-energy D events needed to claim, at the 95% confidence level, a D component of exotic 
origin corresponds to one event for ULDB GAPS/balloon, five for GAPS/satellite, three for a 
GAPS detector mounted on an interplanetary probe, and six for AMS. 

In the right frame of Fig. 1201 the 95% confidence level sensitivies of the various D search 
experiments are shown in the (m B (i), (av)o/m Bil) ) plane. The dashed black line in this figure 
indicates the standard value of the pair annihilation cross section of KKDM neglecting second- 
level KK mode resonances (notice that model UED3 does not lie on the dashed black line 
precisely because n = 2 resonances do contribute). UED will be probed at balloon-borne GAPS 
missions for KKDM masses below ~ 400 GeV in the present halo model setup. Analogously, 
satellite-based experiments extend the reach in the D channel up to masses around « 600 GeV, 
and an interplanetary probe up to m B{ i) rj 800 GeV. 



3.6.2 Antiprotons 

The number of antiprotons (the case of positron is completely analogous) per unit time, energy 
and volume element produced in WIMP annihilations is called the source function, Q, and is 
given by [197J 

Q(T, 2) = Mo ( e*2W) ' £ ^BR(Bm B m I) (75 ) 
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where T stands for the antiproton (positron) kinetic energy, and for a given annihilation channel 
final state, /, dN^ jdT is the fragmentation function into antiprotons (positrons). The effects 
of antiproton propagation through the Galaxy is described in Ref. [197J (see also Ref. |177| ). The 
effects of solar modulation are as described in the previous section for the case of antideuterons. 

In Fig. [2TI the spectra of antiprotons expected for the four benchmark models (UED1- 
4) of Appendix [A] are shown. In the left frame, the spectrum does not include the effects 
of solar modulation, while in the right frame, those effects are included. Also shown in the 
right frame is the expected cosmic ray background |177j and the spectrum of antiprotons as 
measured by balloon-borne experiments [170J. In the right frame, for consistency with the time 
of experimental data-taking for the measured fluxes shown, we make use of an averaged Fisk 
parameter (ftp = 610 MV. 

The spectrum of antiprotons in UED is rather soft and mostly originates from quark- 
antiquark final states. Compared to expectations for neutralino dark matter, the large branching 
fraction of pair annihilations in leptonic states suppresses the antiproton yield. Ref. |198j 
also computed the contribution from the radiative pair annihilation channel into two gluons, 
and found that the resulting branching ratio is smaller than 0.5%. We notice however that 
the ratio of the expected signal over background can be large, especially at very large energies 
(Tp > 50 GeV) , which will be probed at future space-based antimatter exploration experiments 
(see Sec. I3.6.4|) . The low energy region (0.1 < Tp/GeV < 10 GeV), instead, is not only plagued 
with large uncertainties in the background computation and in the effects of solar modulation, 
but also typically features a very suppressed signal-to-background ratio. With smooth dark 
matter halos, currently available data on antiproton spectra (that are limited to relatively low 
antiproton kinetic energies) do not provide significant constraints on the UED model, even in 
the case of a relatively optimistic halo model setup (such as the one used here). 

Ref. [198] . where the antiproton yield from KKDM was analyzed in detail, gets to similar 
conclusions to those sketched above, pointing out that the most promising energy range for the 
detection of a signal over the secondary background in the antiproton spectra to be measured 
in forthcoming space-based experiments lies between 10 GeV and a few 100 GeV. The presence 
of dumpiness in the dark matter distribution in the galactic halo can significantly enhance 
the predicted flux of antiprotons [198]. The resulting increase in the signal depends upon the 
fraction of dark matter mass assumed to lie in clumps, the assumed density contrast, and 
the antiproton energy [198] . and might be crucial for the detection of a WIMP annihilation 
signature, possibly to be correlated with a signal in the positron channel as well (see the following 
Sec. 13.6.3]) . Ref. [199J, which also addressed the flux of primary antiprotons in B^ annihilations 
in the galactic halo, arrives at similar conclusions, emphasizing that the best energy range to 
detect a signal appears to be the high-energy tail, and the necessity of an increase in the 
antiproton flux from a clumped dark matter distribution. The uncertainties in the secondary 
antiproton background in that energy range were analyzed in Ref. |200j and compared to the 
expected primary antiproton fluxes from WIMP pair annihilations (including the case of B^ 
dark matter). 

3.6.3 Positrons 

B^ pairs annihilate very frequently (around 20% of the time) into monoenergetic electron- 
positron pairs with energies at production corresponding to the B^ mass. An additional 40% 
of time they annihilate into either muon or tau pairs which, again, yield energetic electrons at 
production. As pointed out in Refs. [261 1173] . this implies a sizable positron flux from KKDM 
pair annihilations in the local Galactic halo. 

Shown in the left frame of Fig. [22] is the interstellar positron flux expected from the four 
benchmark models (UED 1-4) found in AppendixO These spectra exhibit, especially at large 
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Figure 22: The local interstellar (left) and solar-modulated (right) positron flux, as a function 
of the positron kinetic energy, for the four benchmark models (UED1-4) of Appendix [A] (the 
results for UED2 and UED3 fall almost on top of each other). In the right frame, we also 
include an estimate of the cosmic ray positron background |177j and the current positron flux 
measurements |169| . 



energies, a spectral shape close to a power law, following (for this particular choice of a halo 
model and propagation setup) the approximate scaling 

dN e + „ in7 /lGeV\ 6 /m B( i)\ 3 , _ 2 _ x _, 

, « 3 X 10' u GeV 1 cm 1 s 1 sr \ (76) 

dE e + V m BW J V E e+ J 

The naive scaling with mass of the combination (av) /m 2 B ^ oc fnZ^x) is compensated somewhat 
by the differential number of positrons produced per LKP annihilation, which grows approxi- 
mately linearly with the mass. 

The effects of solar modulation on the interstellar positron fluxes are included in the right 
frame of Fig. [221 Also shown is the expected background |177j and the current experimental 
data on the positron spectrum. We use an averaged solar modulation parameter <j> ~ 490 
MV, appropriate for the period when the data shown were taken. We notice that the positron 
background has a significantly softer spectrum than the signal from KKDM pair annihilations, 
particularly at large energies. At high energies, the signal-to-background can be very large, 
even close to unity. Unfortunately, this occurs at energies which will probably be beyond the 
sensitivity of even space-based antimatter experiments. 

To minimize the effects of solar modulation, a useful quantity to consider is the so-called 
positron fraction, i. e. the ratio, at a given energy, of the positron flux to the sum of the positron 
and electron fluxes. In this context, the background is given by the standard expectation for the 
positron fraction, and the dark matter contribution shows up as an extra component, enhancing 
the positron fraction. 

In 1994 and 1995, the High-Energy Antimatter Telescope (HEAT) reported an excess in 
cosmic positron fraction, peaking in the range of 7-10 GeV, and continuing to higher energies 
|171j . In 2000, an additional HEAT flight confirmed this observation |201j. Many previous 
experiments, although less precise, have also recorded a larger than expected positron flux above 
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Figure 23: Left: The positron fraction from KKDM annihilations as a function of positron 
energy. The solid and dashed lines represent 300 and 600 GeV B^'s, respectively. The annihi- 
lation rate was treated as a free parameter, used for normalization. The dotted line represents 
the background predicted with no contribution from dark matter annihilation. The error bars 
shown are from the 1994-95 and 2000 HEAT flights. Right: The positron fraction from an- 
nihilation of KKDM for several choices of propagation parameters. The solid line represents 
the model with the same propagation parameters as in the left frame. The dashed line is for 
a model with an energy loss rate smaller by a factor of two. The models represented by dot- 
dashed and dotted lines use the full energy loss rate but diffusion constants that are 80% and 
50%, respectively, of the value used in the left frame. Lastly, the spectrum with both half the 
diffusion constant and half the energy loss rate falls on top of the solid line. In the right frame, 
m BW = 300 GeV. 



about 10 GeV (see Ref. [171 J and references therein). The excess of positrons at large energies 
has also been recently confirmed by AMS-01 data, although with a rather poor statistics [202J. 
The study of the astrophysical production of positrons |203j has been thoroughly investigated 
[204J, with the conclusion that the ratio of cosmic positrons to electrons above about 10 GeV 
is higher than is suggested by secondary production in a model of a diffusive halo. 

This excess can be, and indeed was, tentatively interpreted as originating from the produc- 
tion of positrons in Galactic dark matter annihilations (see e.g. Refs. |165t I172j ). Ref. |173| 
explored this interpretation in the context of KKDM. Shown in Fig. [23] is the HEAT data com- 
pared to the expectation from pair annihilations. In the left frame, two values for the B^ 
mass are shown (300 and 600 GeV). In each case, the normalization of the annihilation rate was 
chosen to obtain the best possible fit to the HEAT data. This procedure implies annihilate rates 
larger than those expected from a smooth dark matter halo distribution. To fit the HEAT data 
with KKDM, annihilation rates boosted by a factor on the order of 10 2 are typically required. 

The spectrum shown in Fig. [23] can be modified somewhat through the choice of diffusion 
parameters which are adopted. In the left frame, the following diffusion parameters were used: 
a diffusion constant of K(E e +) = 3.3 x 10 28 (£' e +/l GeV) 0,47 cm 2 /s p05], an energy loss rate 
of b{E e +) = 10- 16 (£ e+ /lGeV) 2 GeV/s and a 2L = 8 kpc thick slab for the diffusion zone 
[205J. In the right frame of Fig. [23] these parameters are varied somewhat to illustrate the 
uncertainties involved. 

Space-based antimatter search experiments, discussed in the following section, will collect far 
better statistics and extend the explored range of positron energies considerably. It is reasonable 



50 



3.6. Anti-Matter Searches 



to expect that PAMELA and AMS-02 will definitively confirm, or rule out, the above mentioned 
HEAT positron excess and its interpretation in terms of dark matter annihilation. 



3.6.4 Prospects For Space-Based Antimatter Experiments 

Upcoming space-based antimatter experiments will tremendously enhance the resolution and 
accuracy of positron and antiproton spectra measurements in comparison to existing balloon 
borne results. The PAMELA experiment [206], on board the Resurs DK1 satellite, was launched 
from the Baikonur cosmodrome on June 15, 2006. Also, AMS-02 |185j is expected to be launched 
and installed on board the International Space Station within the next few years. With the 
purpose of assessing the capabilities of these experimental facilities to identify signatures of 
dark matter annihilations, we follow Ref. |177j . To this extent, we will implement a statistical 
X 2 analysis to compare the case of a pure background measurement to that of a signal from 
dark matter annihilations. 

The relevant experimental parameters entering the estimate are given by: 

- The geometrical factor of the experimental facility, i.e. its effective area, A 

- The time of data acquisition, T 

- The energy coverage of experiment, as defined by the number of and size of the energy 
bins 

We will declare that a given model is discriminable at a certain future experiment at given 
confidence level, X%, if the x 2 induced by the dark matter model is larger than the (x 2 )n^° 
corresponding to degrees of freedom. Letting Nf = Nf + Nf be the number of projected 
signal plus backgroud events in a given bin, i, and Nf be the number of observed events, with 
a standard deviation, A N o , the x 2 is defined as 



"b 



2 



X 



E , V . Nf = N? + Nf. (77) 



A N p 



Supposing that the standard deviation has a Gaussian distribution, one finds 

A N o ~ sJn?. (78) 

We are interested in finding the limiting cases, i.e. those cases for which the signal is a small 
component with respect to the background, as is the case for all KKDM models considered here 
(see Figs. [2T1 and [22]) . We therefore make the assumption that: 

Nf « Nf , or Nf ~ Nf . (79) 



Eq. (|77j) then reduces to 



™ b (N s ) 2 



Since the number of events in an energy bin AE is given, as a function of the flux of particles 
4>, by 

N = (AE) ■ (j) ■ A ■ T, (81) 
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Figure 24: The projected "Visibility Ratio" for the PAMELA experiment, defined as the ratio 
Itj,/ (I<f))o (see Eq. ([85]) ) as a function of the LKP mass, m B (i) , for antiprotons (left) and positrons 
(right). The visibility ratio was computed using a reference halo model (the adiabatically 
contracted N03 profile), and two other halo models (the Burkert profile, and the Moore et al. 
profile with a boost factor of 5 resulting from dark matter substructure). The reach of the 
PAMELA experiment is indicated in each case by a horizontal red line. We also indicate the 
location of three of the benchmark models of Appendix lAl 



and indicating with <p s and 4>b the signal and background fluxes, respectively, Eq. (|80p can be 
recast as 

4>b 



X~ = E ^f--(AE)i-A-T. (82) 



i=l 

The quantity given in Eq. (|82[) is what we use to asses the future sensitivity at antimatter 
experiments. We declare that a model is within discrimination capabilities of a given future 
experiment at the X% confidence level if it satisfies the relation 

X 2 > U 2 )l % - (83) 

We focus here, for definiteness, on the case of the PAMELA detector |206j . and compute the 
reduced x 2 f° r an effective area of 24.5 cm 2 sr, an exposure time of 3 years, and resorting to a 
trial energy binning as sketched in Ref. [207J. The results we show are in the limit of "known 
background", i.e. in the (optimistic) scenario in which degeneracies in the parameters used 
to model the propagation of charged cosmic rays in the Galaxy are resolved, for example, by 
precision measurements of ratios of secondaries to primaries for several light cosmic-ray nuclei 
species. As a rule of thumb, the sensitivity of AMS-02 is expected to be one order of magnitude 
better than that of PAMELA [208] . 

Eq. (|82p can be recast in the limit of a large number of energy bins approximating the 
discrete sum with an integral. The region of integration will be given by the lowest and the 
largest energies, T m ; n and T max , accessible by the experiment. We henceforth define the quantity 



Tmax (4> s ) 2 



T 



dE. (84) 
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The visibility condition, which reflects the continuum version of Eq. (|82|) . is then given by 

( Y 2)X% 

h > (85) 

In the parameter Ij,, the dependence on the extreme of integration is small as long as the peak of 
the signal-to-background ratio falls within them (which is always the case for the models we are 
considering and for the values expected for PAMELA and AMS). The effective area is assumed 
to be independent of energy. If this is not a good approximation, a weighting function should 
be introduced accordingly. For the PAMELA setup already exploited, supposing a number of 
bins given by rej = 55 -f- 60, and taking (x 2 )^" ~ ^5, we find the following critical values for 
J<£, after one and three years of data taking: 

I 1 ^ = 9.7. l(T 8 cm- 2 sr _1 s~ 1 , 1^ = 3.2 • 10 -8 cm _2 sr _1 s _1 . (86) 

In Fig. [2H the reach of PAMELA (after 3 years of data) to antiprotons (left) and positrons 
(right) from KKDM annihilations is shown. The solar modulation Fisk parameter has been 
set to 4>f = 600 MV, a value appropriate for the projected solar activity over PAMELA'S 
data taking time. The reference dark matter halo model is, again, the adiabatically contracted 
N03 halo profile described in Sec. 13.6.11 However, also shown are the results for two other halo 
setups, namely a (more conservative) cored Burkert profile [209] . and a (more optimistic) Moore 
et al. profile featuring a degree of dark matter substructures responsible for a boost factor of 5 
in the antimatter fluxes. 

We find that it will be modestly more easy to disentangle an exotic signature in the an- 
tiproton spectrum than in the positron channel. With antiprotons, KKDM models with masses 
m B (i) < 550 GeV will give a significant signal for the case of the reference halo model (adiabati- 
cally contracted N03). In the case of positrons, this figure is reduced down to m B (i) < 400 GeV. 
The comparison between these two channels is similar when other halo profiles are considered. 

3.7 KK Gravitons 

The inclusion of the gravitational sector within the framework of a model with universal extra 
dimensions yields the existence of a KK graviton tower. Depending upon the KK spectrum 
adopted, the lightest KK graviton, which we shall hereafter indicate with the symbol (not 
to be confused with the KK gluon, g^'), could possibly constitute the LKP. If the LKP is 
the first-level KK graviton, the resulting dark matter phenomenology will be that of a super- 
Weakly Interacting Massive Particle (a "SuperWIMP", see e.g. Refs. EH]). If instead 
KK gravitons are not the LKP, they will be meta-stable, in which case the production of these 
relatively long-lived states can significantly affect the cosmology of the early universe, including 
the relic abundance of LKPs. In Sec. 13, 7. H we consider the phenomenology of models where the 
is the LKP, while in Sec. 13.7.21 we address the consequences of the existence of a tower of 
(unstable) KK gravitons in scenarios where the LKP is not the 

3.7.1 The KK Graviton as the LKP 

As in the case of super symmetry, where the LSP need not be a SM partner, within the context 
of UED the LKP could in principle be the first KK excitation of the graviton. Since radiative 
corrections to the KK masses |30| are typically positive, and they are negligible for the G^ 
(being Planck-scale suppressed), this is a fairly natural option. Assuming the boundary condi- 
tions found in the minimal UED model (see Sec. I2.2.3p . and assuming m G (i) = -R -1 , the mass 
splitting between the and the G^ can be approximated, neglecting electroweak mixing 
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between the KK neutral gauge bosons, as [M] 

2 2 ^39 aiC(3) la! 2 \ /1\ 2 g\v 2 

m Bm ~ m GW \Y^3~ 64^ HAR) ) [r) + ~ (8?) 

Making use of the full expression, the mass splitting is positive (hence the is the LKP) for 
R' 1 < 809.1 GeV, and negative elsewhere (hence the is the LKP for R' 1 > 809.1 GeV). 
Also, again within the context of minimal UED, one finds that the maximal mass splitting, for 
R ^ 809.1 GeV, amounts to only (m B (i) — m G (i)) < 1.7 GeV. We stress, however, that the 
minimal UED spectrum can be affected by further contributions to the KK masses, e.g. from 
large boundary terms at the cutoff (relaxing the assumption that the boundary kinetic terms 
vanish at the cutoff scale, A). 

The production of KK gravitons in the early universe proceeds in general through two mech- 
anisms: (1) non-thermal production through NLKP decays, and (2) production at reheating 
after the end of the inflationary era. In the first case, the NLKP undergoes the usual thermal 
freeze-out when its pair annihilation rate drops below the Hubble expansion rate, and subse- 
quently decays into KK gravitons. As the NLKP-LKP interactions are suppressed, the NLKP 
freeze-out occurs at much earlier times than the NLKP decays into the KK graviton. In this 
scenario, since every NLKP decays into one KK graviton, the KK graviton relic abundance is 
simply given by 

m G( i) 



$7 G (i) = Onlkp 



"1NLKP 



Since the KK graviton mass and the NLKP mass are generically predicted to be rather close to 
one another, if the NLKP relic abundance is similar to the CDM abundance the KK gravitons 
will inherit the correct relic abundance through NLKP decays. As long as the reheating tem- 
perature is low enough to avoid regenerating large numbers of KK gravitons, this is the main 
production mechanism for LKP KK gravitons. 

In any case, KK gravitons are also produced at reheating. If reheating occurs at a temper- 
ature, 1rh> a population of KK gravitons is produced through the collisions of particles in the 
thermal bath. The evolution of the number density of nth level KK gravitons, n G ( n ), follows 
the Boltzmann equation, 

— + 3Hn G ( n) = C G ( n ) , (89) 
where C G n is the collision operator and can be parameterized as follows: 

C G( n) =Ca(g*(T)n ) 2 , (90) 

where 

a = iSf' ^-(levrGV)- 1 (91) 

and where is the strong coupling constant. In the equation above, C is the graviton pro- 
duction parameter that can be understood as the fraction of all possible collisions which will 
interact strongly to produce gravitons, g*(T) ~ g^ K TR is the number of effective degrees of 
freedom (g^ K is a model dependent quantity that for the case of one extra-dimension compact- 
ified on the S 1 jZi orbifold equals approximately 197.5), and no = C(3)T 3 /7r 2 . The cumulative 
effect of all KK levels up to a KK mass of the same order of Trh gives, again for the S 1 /Z% 
UED model and for m GW = R' 1 [212] . 

agUo^xCxfCy'V— V(^y /2 . (92) 
V 200 / \m G (i)J \m G mJ 
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Figure 25: The lifetime of NLKP in the KK graviton LKP scenario, as a function of the mass 
splitting Am between the NLKP and the G^ 1 ' . In the inset we show the relation between R^ 1 
and Am in the minimal UED scenario. 



The value of C is subject to some controversy in the literature. Ref. |212j estimates C ~ 
0(0. 01) based upon the analogy with the supersymmetric gravitino case worked out in detail in 
Refs. [213} 1214] . The independent estimate of Ref. [53] . in contrast, yields C ~ 0(1). As a rule 
of thumb, when Trh ^ 10 2 TeV, LKP production from the KK graviton tower is expected to be 
the dominant production mechanism over the G^ production from long-lived NLKP decays, 
and vice- versa. A lower limit on Trh ; in the case of production from long-lived NLKP 
decays, stems from the requirement that the universe is not reheated after thermal freeze-out, 
i.e. T RH >iT 725. 

Constraints on the KK graviton LKP scenario come from limits on the amount of energy 
injected through NLKP decays into G^s. NLKP decays at early times can affect the predictions 
of Big Bang Nucleosynthesis (BBN), the Cosmic Microwave Background (CMB) or conflict with 
observations of the diffuse photon background [215] . The first necessary step to determine these 
constraints is to compute the lifetime of the NLKP [2T0| [2T2] . If the NLKP is the the 
resulting decay width is given by 



(i) 



cos 2 6 W rn 7 B (i) 



m 



144vrM 2 mi 



1 + 3- 
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+ 6- 
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G(i) 



G<!) 
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(93) 



BW 



while if the NLKP is a KK chiral fermion, /j 1 ^, one has [212] 



(i) 

L,R 
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L,R) 
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GW 
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/(i) 



(94) 



We show in Fig. [25] the predictions for the NLKP lifetime as a function of the NLKP-C7W mass 
difference, Am, in the case of a NLKP, at three different values for the G^ mass (two of 
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Figure 26: The constraints on the R ^Am parameter space of the KK graviton LKP mass 
versus the mass splitting from BBN (the excluded region is shaded in light blue), 

CMB (light grey) and the diffuse photon background (yellow). The black line indicates the 
minimal UED prediction for the B^-G^ mass splitting as a function of R" 1 . 



them feature the same value of -R -1 as that of models UED1 and UED4), and for R^ 1 = 1 
TeV for the case of a chiral fermion NLKP. We also show the prediction for the B^ lifetime in 
the case of minimal UED, where the B^-G^ mass splitting, as a function of is given in 
the inset. The hatched region near the top of the plot contains an NLKP with a lifetime larger 
than the age of the universe (hence the possibility of a mixed G^-B^ dark matter scenario. 
A KK fermion NLKP in this region is ruled out because either a fraction of the dark matter 
would be charged or, for a KK neutrino NLKP, would be excluded by the negative results from 
WIMP direct detection searches). 

The most relevant bounds on the KK graviton scenario depend upon the total energy release 
into photons in B^ NLKP decays, i.e. on the quantity e 7 Yy, where e 7 is the energy of the 
photons when created and = n 7 /n 7 G is the number density of photons from NLKP decay 
normalized to the number density of background photons n 7 G = 2£(3)T 3 /it 2 , where T is the 
temperature during NLKP decay. Since NLKPs decay essentially at rest, e 7 = (m-NLKP — 
mg (1) )/(2mNLKp), and since one KK graviton is produced in association with each photon, 
Yy = Y G (i) , the photon abundance is given by [210] 



Y GW * 3.0 x ID- 2 r-^L . ( 95) 



-1.2 ( \ TeV \ ( 

m Gm J V 0.23 

Fig. [25] shows that NLKP decays can occur before and during BBN, hence they can in 
principle destroy the successful BBN predictions for the abundances of light elements. Pho- 
tons produced in NLKP decays rapidly thermalize through scattering off background photons, 
77bg — ► e + e~, as well as through Compton scattering. The BBN constraints therefore very 
weakly depend upon the initial energy distribution of the injected photons. We show in Fig. 1261 
(adapted from Ref. 210 ) the region of the i? _1 ,Am parameter space ruled out by the BBN 
constraint [216] with a light blue shading. In the figure, we assume fi G (i) — ^cdm- 
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The injection of photons can also distort the CMB spectrum from the observed blackbody 
distribution as a result of elastic Compton scattering, bremsstrahlung and double Compton 
scattering, e~7 — ► e~77. Before redshifts z ~ 10 5 , the bounds come from corrections to the 
chemical potential [i of the Bose-Einstein distribution to which the photon spectrum relaxes. At 
redshifts smaller than z ~ 10 5 , deviations from the black-body spectrum can be parameterized 
through the Sunyaev-Zeldovich y parameter |217j . The resulting bounds [210J exclude the region 
shaded in light grey in Fig. [26l 

The non-thermal production of soft gamma rays in NLKP late decays contributes to the 
diffuse photon spectrum, with a differential photon flux peaking at an energy 

£max = 490 keV fll^A (96) 
' \ Am J 

and giving a maximal differential photon flux of [210| 

d®i , ™^ ZlTeV\ / Am \ fSl r m\ 911 1 

— l (ET x ) ~ 2.1 x ^7T7 cm" 2 s" 1 sr" 1 MeV" 1 . (97) 

dE^ 7 ' \ m G (i) / \ 1 GeV / \ 0.23 / v ; 

Constraints from the observatories HEAO, OSSE and COMPTEL |218j translate into the ex- 
cluded region shaded in yellow in Fig. [26j 

In the figure, we also show the prediction for the B^'-G^ mass splitting as a function 
of R~ l in the minimal UED model (where, recall, the can be the LKP for m G {i) < 809.1 
GeV) . We conclude that the a graviton LKP is not a viable option within the context of minimal 
UED, as it badly violates both the diffuse photon background and the CMB bounds. A possible 
way out was proposed in Ref. |219j . through the introduction of right-handed neutrinos to the 
model. The B^ then dominantly decays into neutrinos, and does not emit photons. Within 
this context, models with sub-dominant decay modes producing photons have been shown to 
feature branching ratios so small that the CMB and diffuse photon background spectra are not 
distorted at observable levels [219J. 

As a last remark, we wish to shortly comment on the possibility of having detectable signa- 
tures from the KK graviton dark matter scenario. In the case that the NLKP is neutral (which 
appears to be the most natural option), the presence of KK graviton dark matter might reveal 
itself only in small perturbations to the quantities listed above (CMB, BBN) or through spec- 
tral features in the diffuse photon background. If the NLKP is instead charged, other detection 
methods based on the production of the meta-stable charged NLKP at colliders or in high en- 
ergy neutrino-nucleon collisions in the atmosphere could lead to detectable signatures, at least 
in principle. Ref. |220j proposed to build water tanks surrounding the LHC to trap long-lived 
charged particles (the charged NLKP), which could afterwords be studied in their decays in an 
appropriate low-background environment. High energy neutrinos colliding with nucleons in the 
atmosphere can also produce long-lived charged NLKPs, which can in principle travel through 
the Earth long enough to be eventually detected at neutrino telescopes |221^ 1222] . The tracks 
resulting from charged NLKPs can be disentangled from the di-muon background on the basis 
of both their spatial separation (muons must be created close to the detector, unlike charged 
NLKP) and on their energy spectrum (charged NLKPs feature a significantly more energetic 
spectrum) |222j . 

3.7.2 KK Gravitons and the Early Universe 

Even when the LKP corresponds to the B^\ its relic abundance will be, in general, affected 
by B^ production through the decays of the KK tower of gravitons. Indicating with 
the abundance of LKPs without the inclusion of gravitons, the total relic abundance can be 
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Figure 27: Left: Values of the reheating temperature, Trh> obtained by demanding a proper 
dark matter density, Q B {i) = Odm ~ 0.23, for different values of the graviton production 
parameter: C = 1, 0.1, 0.01. Also shown are lines of constant ratios of the reheating temperature 
to the lightest KK mass: Trh = 40 rn-KK; 70 mKK and 100 ttikk- Right: Constraints on the 
value of the mass difference, Am, between and the LKP assuming £l B (i) = S^dm ~ 0.23 
and C = 1, 0.1, 0.01. Possible coannihilation and second KK level effects have been neglected. 
Figures adapted from Ref. [53] . 



expressed as 



Y BW = Y^ + Yg (98) 
n B w = s Y B(1) (99) 

n BW = mBW nBil) (ioo) 



where p c = 5.3 x 10 TeV/cm , and 



Y G = nY G( n) dn 

Jo 



_ !^!2)-^ff (TRH ^. (101) 

7 7T 8 Mi R 

Requiring that £l B (i) ~ 0.23, one gets a prediction for the reheating temperature of 

1 2 /7 

Tw ^ R X { 1.8x10-18 CsoR-z ) ' (1 ° 2) 

where sq is the entropy density today. The left frame of Fig. [27] illustrates, as a function of 
ttikk = -R , the reheating temperature needed to produce B (i) 0.23 for different values of 
the graviton production parameter C. 

The constraints from the diffuse photon background and from BBN mentioned above apply 
also in the case of a meta-stable KK graviton, , decaying into the LKP. In particular, 
since photons from n > 1 KK gravitons are produced before matter domination, they don't 
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affect the diffuse photon background. On the other hand, a line of reasoning analogous to that 
presented above leads to a constraint on the mass difference between the lightest KK graviton 
and the LKP, Am, which reads 

Am > ^2.48 x 10- 3 ((T RH i?) 3/2 - l) (pf^)) ' GeV > ( 103 ) 

where it was assumed that = Hdm ~ 0.23. 

Constraints on BBN come from electromagnetic energy released by the decay of the — * 
B^'j, and from the hadronic decays of the whole KK graviton tower (the G^ n 's, n > 1, decay 
right after BBN). The latter bounds depend strongly on the branching ratio of the decay of 
gravitons into hadrons. Indicating with B^M/Had the electromagnetic/hadronic branching ratio, 
one obtains the following constraint on Am (again assuming = Qdm ~ 0.23) 



Am < — — ^ (104) 

B EM/Had \Pc ~ R SO ^oo) 



where £b is the bound on the energy released per background photon. Conservatively, Ref. 
assumed £ B < 10~ 15 TeV. 

The bounds quoted above are shown in the right frame of Fig. [27] for the same model 
considered in the left frame, for three values of the graviton production parameter, C. As 
pointed out in Ref. [Mj; the prediction of the minimal UED model for Am can be, in principle, 
compatible with the bounds in Eq. (|1 03H104|) within small ranges of the inverse compactification 
radius. 

Setups that mimic the features of a UED scenario without KK gravitons are the so-called 
"deconstruction models", where KK states are merely the manifestation of a chain of gauge 
groups \223\ [224J. For instance, a simple two-site model could resemble the first KK level of 
UED without the introduction of a KK graviton. 

Scenarios where the UED model is encompassed in a higher dimensional setup, and where 
all SM particles are confined to the five-dimensional space-time but not the gravitons (which 
therefore can be significantly heavier than the rest of the KK spectrum), also effectively get rid 
of the implications of KK gravitons discussed above. One example was given in Ref. [225] in 
the context of shape moduli associated to large extra dimensions. 

The cosmology of KK gravitons was also discussed in Ref. |226] , with similar results to those 
outlined above. In that analysis, UED and other models featuring flat extra-dimensions have 
been studied in a cosmological context addressing, in particular, the issue of the cosmology of 
gravi-scalars, or "radions" , i.e. the role of the scalar field associated to the geometrical modulus 
that determines the size of the compactified dimensions. The main parameters entering radion 
cosmology are the compactification radius and the scale Mj = of the inflaton potential. 
Different cosmological scenarios are then discussed, depending on whether inflation took place 
before, after or around the same time of the compactification of the extra spatial dimensions. 
Since, on rather model independent grounds, the radion decay width reads 



647rm|j ' ^ ^ 

the radion is effectively stable (i.e. its lifetime is larger than H^ 1 ~ 4 x 10 17 s) if R^ 1 < 7 x 10 8 
GeV. This implies a stable radion for the UED model discussed here, and, consequently, a generic 
over-closure problem, as the radion can easily dominate the energy density of the universe at 
late times. The radion energy density can be sufficiently damped if the inflation scale Mj lies 
around the TeV scale. As the radion is stable, reheating would have to come from a different 
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field, which would also imprint its fluctuations into the cosmic microwave radiation. Moreover, 
Ref. [226\ points out that there is a narrow window of parameter space where the radion can 
be the dark matter as well. In general, inflation model building for ~ TeV is problematic. 
Specifically, a low reheating temperature, as needed to account for the dilution of the radion 
energy density, necessitates a very efficient reheating, and almost certainly a preheating phase 
as well [226] . 

Along similar lines, Ref. [227J investigated the case of n universal, "small" extra-dimensions 
and p larger spatial dimensions where gravity also propagates. This setup is motivated by the 
idea of explaining the weakness of gravity by the presence of large extra-dimensions orthogonal 
to the 4 + n dimensions of the Standard Model brane. The model features two radions, that can 
be, in turn, identified with a "shape" and a "volume" mode, respectively. Ref. [227j notices that 
models where a radion plays the role of an inflaton, or where the inflaton is a brane scalar field, 
generically present problems; however, it is possible to find setups where bulk scalar fields drive 
successful inflation, although the success of radion cosmology depends upon a more complete 
understanding of the radion potential arising from the fundamental theory governing the low 
energy action. 



60 



4 Universal Extra Dimensions and Colliders 

The KK states of UED can be produced in high energy colliders such as the Tevatron, the Large 
Hadron Collider (LHC), or a future International Linear Collider (ILC). At hadron colliders, 
first level KK modes can be abundantly pair-produced, depending upon the value of the com- 
pactification scale. The colored KK modes then cascade decay according to decay chains similar 
to those shown in our Fig. [5] of Sec. 12.41 Since the B^ is stable and weakly interacting, it will 
escape the detector, leading to the typical signature Z + Z~Z+ 2 jets + missing energy. Given the 
approximate degeneracy of the KK modes, the resulting jets will be rather soft and challenging 
to distinguish from the background. The leptons will also be soft, but usually they can pass 
some reasonably chosen cuts. In the following section (Sec. 14. ip . we review the prospects for 
the discovery of UED at the LHC. 

The advantage of a linear collider featuring y/s > 1/R over an hadronic collider would be 
enormous, and would entail the possibility of a very accurate reconstruction of the KK particle 
spectrum, as well as the unambiguous identification of the spin of the KK modes. We devote 
Sec. 14.21 to the discussion of the collider physics related to UED that could be performed at 
a future ILC. We remind the reader that the current constraints from the Tevatron and the 
discovery potential of Run II were already reviewed in Sec. 12.41 

A crucial issue in the assessment of the potential of collider experiments in exploring various 
beyond-the-Standard-Model scenarios is how to differentiate among them. As testified by the 
catchy title of the seminal paper, "Bosonic supersymmetry? Getting fooled at the LHC" |29j . if 
new physics signals are detected at colliders, a basic point will be how to discriminate between 
supersymmetry and UED (or other TeV-scale physics scenarios). In a nutshell, although the 
expected collider signals from the two setups are strikingly similar, there are three basic handles 
to potentially discriminate between the two: (1) KK first level states in UED have the same 
spin as their SM counterparts, while SUSY partners have opposite spin, (2) the Higgs sector 
of UED carries a different KK parity assignment than the heavy Higgs bosons of the MSSM 
(H, A, H ), making them more similar to the SUSY higgsinos than to the SUSY Higgs sector 
(although the two Higgs sectors in SUSY and UED share exactly the same gauge quantum 
numbers), and (3) UED feature higher level KK modes, unlike the case of supersymmetry. We 
devote the following Sec. [5] for a thorough discussion of these points, as well as of the role, in 
the discrimination between supersymmetry and UED, of dark matter search experiments. 

4.1 The Large Hadron Collider 

The signatures of the production of new heavy particles at collider experiments crucially depend 
upon those particles' interactions and mass spectrum. Here, we focus on the particular case 
of minimal UED [30J, where all boundary couplings are assumed to be flavor-conserving, and 
all the boundary terms are assumed to vanish at the cutoff scale, A (following Ref. |29j . we 
shall adopt here AR = 20). Once the radiative corrections are taken into account, the KK 
mass degeneracy is lifted, and the heavier KK states promptly decay. The resulting collider 
phenomenology is similar to the case of supersymmetric models where all superpartners lie 
relatively close in mass. As in the case of supersymmetry, the decay cascades will terminate in 
UED with the stable lightest KK particle leaving the detector undetected, resulting in generic 
missing energy signatures. 

At hadron colliders, the by far most abundantly produced KK particle species will be KK 
gluons and KK quarks, i.e. the strongly interacting states. Fig. [28] shows the production cross 
sections for strongly interacting first level KK states at the LHC (y/s = 14 TeV). As a first step 
to assess the capability of extracting a signature of UED at the LHC, one needs to consider the 
allowed decays of of the first level KK particles and estimate the relevant branching fractions. 
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Figure 28: The production cross section of KK quarks and KK gluons at the LHC. The solid 
curve represents the total production cross section, while the remaining curves show the separate 
contributions from the production of KK quark pairs, KK quark + KK gluon, KK gluon pairs 
and KK top pairs. Adapted from Ref. |28j . 



The heaviest KK particle in the present setup is the KK gluon, The two-body decays 
to the KK quarks u^ 1 ', are always open and have comparable branching ratios, 



BR(g 



(i) 



BR( 5 



(i) 



0.5, where q = u,d. 



(106) 



The 5'C/(2)-singlet KK quarks, q( \ can only decay to the hypercharge gauge boson, 
their branching ratio into Z^ = W 3 ^ being suppressed by the level-1 Weinberg angle (recall 
that 0$ < 9w, leading to BR(gW -> Z^qW) « sin0$ « 10~ 2 - 10~ 3 , while BR(gW -» 
B^q^) 1). As a consequence, production yields jets plus missing energy. An exception 
to this situation occurs for the KK top quark, , which can decay to W + ^b^ and H+^b^ 
(as noticed in Ref. [29], this is also the dominant mechanism for the production of H^ 1 ^ at 
hadron colliders). 

The decay chain of SU(2) doublet quarks, Q^ 1 ', is somewhat less trivial, as they can decay 
into a SM quark plus or W^. In the sin#^ <C 1 limit, SU(2) symmetry dictates 



BR(Q« _> W^Q'W) « 2 x BR(Q« -> Z^Q^) 
For mn(o) C 1/]J one also has 

BR(Q( 1 ) -> Z0)Q(0)) 9lTl Q {m 



(107) 
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BR(Q« -> B(i)Q(o)) <? 2 l^(m 2 



(108) 



where T3 and Y stand for the weak isospin and hypercharge of the quark. The KK quark decays 
into 577(2) gauge bosons, although suppressed by phase space, are enhanced by the ratio of the 
couplings squared as well as by the quantum numbers. Typically, one numerically obtains 

BR(Q« -> W {x) Q m ) « 65%, BR(Q« -» Z^qM) « 33%, BR(Q^ -» B^qC)) « 2%. 

(109) 
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Once produced, the heavy KK gauge bosons cannot decay hadronically for kinematic reasons, 
and will democratically decay into all lepton flavors (we indicate here with L the charged SU (2) 
doublet lepton of any generation, and with I the SU(2) singlets), 

BR^ 1 ) -> „Wl±(°)) « BRCW^M -> L^V )) « i (110) 
BR(Z« - ^i/C)) « BR^ 1 )^^ 1 ^ ))^! (Ill) 

for each generation. The branching fraction for —* Z^)^ ) is suppressed by sin#^ <C 1. 

The resulting level- 1 KK leptons will directly decay into B^- 1 ' and the corresponding SM 
lepton. As a result, the heavy KK gauge bosons will almost always decay as 

W MD _> B (i) L ±(o) u (o) and ^±(1)^^(1)^(0)^(0) and z ±(i) _ 5(1)^(0)^(0) (112) 

resulting in large e and // yields. Lastly, the KK Higgs will decay according to the available 
channels, dictated by the mass of the SM Higgs. If the KK Higgs is heavier than the KK W, 
Z, t and/or 6, it will decay into those states. Otherwise, the tree level two-body decays will be 
suppressed, leading it to decay to B^ and a virtual zero-mode Higgs boson, or through a loop 
to -E^ 1 ) and a photon. In any case, the production cross section for the KK Higgs is negligible 
with respect to the strongly interacting KK states. 

At the LHC, the signature with the largest overall rate is E^f nss + (N > 2) jets, analogous to 
the well-known traditional squark and gluino searches |228| . Roughly one quarter of the total 
production cross section of KK states at the LHC gives rise (directly or indirectly) to pairs 
which produce the aforementioned signature. However, even though the missing mass in these 
events is rather large, the measured missing energy is rather small, as it is correlated with the 
energy of the relatively soft recoiling jets. The estimate of the LHC reach in this channel is 
around 1/R <1.2 TeV, as can be inferred from the analysis carried out in Ref. [229J. 

A cleaner channel for UED discovery at the LHC is that of multilepton final states arising 
from associated heavy KK gauge boson production. The inclusive pair production cross 
section amounts to roughly one quarter of the total production cross section. The subsequent 
decays of the Q^s yield W^W {1 \ W^Z^ and Z^ZW pairs in the approximate proportion 
of 4:4:1. The subsequent heavy KK gauge boson decays yield final states with missing energy 
and up to four leptons, each of which offer the possibility of discovery. Following Ref. [29], we 
consider here the gold-plated 4/ + E™ 1SS signature. 

Ref. |29] studied the reach of the LHC in this channel. They conservatively neglected direct 
KK gauge boson production and Q^W^ and Q^Z^ associated production processes and 
used for the single lepton trigger the relatively soft cuts of px > {35, 20, 15, 10} GeV with 
\rj(l) < 2.5|, and, in addition, Ef iss > 50 GeV. The SM background from ZZ -> I^Pt+t- -► 
41 + E™ ss , where Z stands for a real or virtual Z or 7 [230] . can be reduced by invariant 
mass cuts for any pair of opposite sign, same flavor leptons. To this extent, Ref. [29J uses 
\mu — Mz\ > 10 GeV and mu > 10 GeV. Additional sources of background, including multiple 
gauge bosons and/or top quark production [230], fakes, and leptons from 6-jets, are estimated 
conservatively to give a background of the order of 50 events after cuts per 100 fb _1 . 

The estimate of the LHC reach in the (1/R, L) plane (where L stands for the total integrated 
luminosity) is shown in Fig. [29l where the two requirements of 5 signal events and of a 5<r 
significance over background are shown separately. The Tevatron Run II performance is also 
estimated, making use of similar cuts and assumptions [29J. The final LHC reach in the 4/+£'™ lss 
channel extends up to inverse compactification radii just below roughly 1.5 TeV. 

Other leptonic channels such as two or three leptons plus E™ 1SS are affected by larger back- 
grounds, but take advantage of the larger branching fraction for — > W^Q'^ and can 
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Figure 29: Discovery reach for the minimal UED scenario at the Tevatron (blue) and at the LHC 
(red) , in the 4£+ $t channel. The reach contours correspond to a 5<r excess or the observation 
of 5 signal events. The y-axis corresponds to the required total integrated luminosity per 
experiment (in fb _1 ) as a function of R~^, for AR = 20. Both the 5<r excess and the 5 signal 
events refer to single LHC experiments, and not to the combination of the two). From Ref. [29J. 

therefore offer higher statistics which can be useful, for instance, to enhance the reach of Teva- 
tron Run II. 

It is important to realize that the results presented above are somewhat model-dependent, 
in that if the cutoff scale was smaller than the assumed value, AR = 20, the mass splitting 
between the first level KK excitations would be reduced, along with the energy of the final 
state leptons, making it harder to differentiate the signal from the background |24j . Moreover, 
should the boundary mass terms be non-zero at the cutoff scale, the phenomenology might also 
be significantly affected [23]. For example, if the KK leptons are heavier that the heavy KK 
gauge bosons, then the cascade decays will be significantly modified from the minimal UED 
case, as will the KK gauge boson decay widths. In such a scenario, new competing decay modes, 
such as — ► W^B^ — ► lvB^\ though suppressed by sin0^ <C 1, would be competitive 
with the off-shell decay process through a KK lepton. 

4.2 The International Linear Collider 

An International Linear Collider (ILC) with a center of mass energy smaller than 3 TeV would 
not be, within the context of UED, a discovery machine. As pointed out before, a signal from 
UED is in fact expected at the LHC when R^ 1 < 1.5 TeV, at least in the multi- lepton plus 
missing energy channel. However, given the cleanness of the environment in e + e~ collisions, an 
ILC would play an invaluable role in accurately measuring the properties of the particles under 
consideration, including their mass, couplings and spin. Such determinations will be critical for 
differentiating between UED and other beyond-the-SM scenarios (see Sec. [5]). 

As alluded to above, a crucial test for the discrimination of UED from other scenarios is 
the possibility of producing KK level- 2 modes. If the ILC is energetic enough to pair produce 
first level KK modes, than single production of second level KK states should also be possible. 
As single production of KK level- 1 modes is not allowed by KK parity, the ILC will face a 
none-or-both situation. We devote Sec. 14.2.11 to the discussion of the production and detection 
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Figure 30: Left: Cross section versus 1/R for the process e + e~ — ► e + e~+ missing energy, in 
minimal UED for various polarization choices and values of AR. Right: The forward-backward 
asymmetry, Aps, versus 1/R for the same process (from Ref. |231j ). 

of level-1 KK modes at the ILC }231| . while Sec. 14.2.21 deals with the physics of singly-produced 
level-2 KK modes |232] . 

If R^ 1 is large and the mass splitting between the heavier KK modes and the LKP is small, 
a serious problem a future high energy ILC will face is that of the huge background from two 
photon processes for a signal consisting of very soft final state leptons or jets [233]. The latter 
are produced by the collisions of two soft photons from initial state radiation, while the e + and 
e~ go down the pipe, resulting in large missing energies. The leptons coming from the decays 
of KK leptons can thus potentially be buried by the two photon background. Larger mass 
splittings than those predicted in the minimal UED model, or the production of KK quarks, 
can alleviate this problem. Another way out is to use e~e~ collisions |233j . In this case the 
signal would be two soft same charge fermions plus large missing energy, and the opposite 
charge soft fermion background from soft initial state radiation would not be as problematic 
as with e + e~ collisions, making it possible to identify and study in detail the properties of KK 
electrons, as well as, in principle, to test the possibility of lepton flavor violation in UED |233j . 

4.2.1 KK n = 1 Pair Production and Detection at the ILC 

Following the analysis presented in Ref. (231) . we assume an ILC featuring a center-of-mass 
energy of -y/i = 1 TeV, and consider the range 250 GeV < R' 1 < 450 GeV within the minimal 
UED framework. The production cross section for SU(2) doublet and singlet KK electrons is 
shown in the left frame of Fig. [30] for different polarizations of the incident beams (including 
the optimal ILC polarization of 80% for e~ and 50% for e + ) and for different values of the 
cutoff parameter, AR. The collider parameters employed in the simulation are specified below. 
The e + e~ — > E^^E L ^ process proceeds through s and t channels, the first mediated by 7 

and Z exchange, and the latter by B^,Z^ and A^> exchange. The KK electrons promptly 
decay into an electron and an LKP. The mass splitting is large enough for the decay to occur 
well within the detector with a nearly 100% branching ratio. The observation of a displaced 
vertex might, in principle, be possible for very small mass splittings. The final state signature 
will then be e + + e~ + 2B^\ the latter resulting in missing energy. 

The same final state can be obtained from e+e" W+^W'W as well. The can then 
decay into a KK neutrino and a lepton, or a KK lepton and a neutrino. Since all three lepton 
flavors are open, and assuming that only the final state electrons are tagged, the e + e _ — > 
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Figure 31: The initial state radiation background from the 2 photon processes (from Ref. [233J). 

p^+Wp^C 1 ) reaction is suppressed by a branching fraction of 1/9 to the final state e + + 
e _ + missing energy and, therefore, provides a subdominant contribution. An even smaller 
contribution comes from the charged KK scalars, A^ 1 ^. 

The dominant SM background, as explained above, comes from initial state radiation 
photon-photon to e + e~ events j233|, 1231] . We illustrate the process in Fig. [31] The 7*7* — * e + e~ 
production cross section at a 1 TeV ILC is expected to be around 10 4 pb, with about half of 
those events resulting in visible particles. The background e + e~ pairs are usually quite soft 
and coplanar with the beam axis [231] . Therefore, an acoplanarity cut significantly helps in 
removing this background, while not significantly reducing the signal [231] . Ref. [231] quotes 
that, for instance, excluding events which deviate from coplanarity within 40 mrad reduces the 
signal by only 7%. Another strategy to estimate and eliminate the aforementioned background 
would be to derive it from observations of \x r events, whose production cross section is sup- 
pressed by a factor 20 with respect to the e + e~ case (due to the s-channel suppression). Once 
measured, this could be then be used to normalize the actual background. Ref. [231] estimates 
that, for an integrated luminosity of 300 fb _1 , the signal events divided by the square root of 
the number of background events would still be around 10 for the range of R^ 1 considered. 
Less significant backgrounds originate from e + e~ — > W + W~ , evW and e + e~ Z, followed by the 
appropriate leptonic decays of the gauge bosons. The cross section shown in the left frame of 
Fig. [30] also assumes kinematic cuts on the lower and upper energies of the final state charged 
leptons of 0.5 and 20 GeV, respectively. The lower cut is the minimum energy resolution for 
identification, while the upper cut eliminates most of the SM backgrounds. A rapidity cut is 
also assumed, admitting only those final state electrons away from the beam pipe by more than 
15 degrees. 

A further handle to deal with the SM background is to make use of the forward-backward 
(FB) asymmetry, defined as Afb = (&F ~ &B)/(c r F + &b)- For the SM background, Afb = 0, 
while the result for KK electrons, as a function of R~ l for various values of AR, is shown in the 
right frame of Fig|301 Beam polarization does not seem to help in enhancing Afb- However, it 
is important to bear in mind that the main background (electron-positron pairs from two initial 
radiation photons, see Fig. I3ip . is forward-backward symmetric, and therefore does not pollute 
the new physics signal from the forward-backward asymmetry. 

The phenomenology of the KK Higgs sector at the ILC was investigated in Ref. [234] , where 
both the importance of soft r leptons detection and the role of backgrounds from pair production 
of KK gauge bosons and KK taus were emphasized. The most promising channel appears to 
be charged scalar pair production. Ref. [234] also concludes that, in the context of the minimal 
UED model, the determination of the value of the parameter rhh from collider data will be 
particularly challenging. 

4.2.2 KK n = 2 Pair Production and Detection at the ILC 

The production and observation of n = 2 KK states would constitute a highly distinctive 
signature for many extra-dimensional scenarios, particularly for the case of UED [29]. While 
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pair production of such states would be challenging at the energy of the LHC (and is surely 
beyond the reach of an LLC), the resonant production of a sing le 5( 2 ) or Z( 2 ) would appear 
as two narrow, closely spaced peaks [235] 1236] . Although this will likely not be resolvable by 
the LHC, particularly because, as we show below, the B^ decay pattern is dominated by two 
jets plus no missing energy (see, however, the discussion in Sec. [5] for prospects with other final 
states), an ILC in the y/s scan mode is potentially the ideal environment to produce and study 
n = 2 KK modes. Due to the structure of the KK spectrum, the production of n = 2 KK states 
at an ILC stands on the same grounds as the pair production of first level KK states. As single 
production of n = 1 modes is forbidden, the ILC will face a none-or-both situation. 

The decay of an n = 2 KK state to two n = (zero mode) states is allowed by KK parity 
conservation, but suppressed by the boundary-to-bulk ratio, as it is KK number forbidden. 
However, there is no phase space suppression, not even if the final state is a tt pair. On the 
other hand, KK number conserving decay modes (2 — ► 2 + 0, 1 + 1) have large couplings, but 
are kinematically suppressed, if kinematically open at all. 

The couplings of n = 2 KK gauge bosons to fermion-antifermion pairs are given by |30j : 

where stand for the U(l)and SU(2) gauge couplings, respectively, T a is the group generator 
(hypercharge and third component of isospin, respectively) and P+ is the ^-even projection 
operator (which is Pl = (1— 7s)/2 for Z^ 2 \ but can be both Pi or Pr for B^). The expressions 
for the boundary corrections, 5, are given by the part of the expressions in Eq. ([8]) proportional 
to ln(Ait!). Eq. (|113p includes the sum of four contributions: (1) one-loop vertices, (2) n = 
2(external)-n = gauge bosons kinetic mixing, (3) (n = 2)-(ra = 0) gauge bosons mass mixing 
and (4) (n = 2)-(n = 0) fermion mass mixing (see Fig. 11 and 12 and App. C in Ref. [30] for 
more details). 

For any level, B^> is the lightest KK state. Thus, B^ cannot decay into other n = 2 states 
plus SM particles. It turns out that the decay into n = 1 pairs is also kinematically forbidden 
for I?( 2 ), leaving as the only possibility the KK number violating decay into zero-mode fermion- 
antifermion pairs. The left frame of Fig. [32] shows the coupling 



y/2 ( Km 2 B(2) z(2) ) _ S(m f . 



X Vf = \\ ~ 2 ^ I > (H4) 



^ (2ft- 1 ) 2 2R 



as a function of AR (notice that for AR = 2, 5 = 0, see Eq. ([8])). Here, V denotes B^ or 
Z( 2 ). As can be seen from the figure, the B^ decays almost entirely into qq pairs, giving as 
a signature two jets with no missing energy. As pointed out in Ref. [232], B&' cannot decay 
into KK-number conserving three or four-body channels, e.g, B^ — > E^E^ 1 '* — > E^eB^ l \ 
at least for the spectrum of the minimal UED type. The reason is that KK-number conserving 
decays must result in two LKPs in the final state, and 2m B (i) > m B {2) over the entire parameter 
space. 

The decay pattern of Z^ 2 ) is more complicated [232J. The is an almost pure VF-like 
state and, therefore, couples only to left-handed doublet fermions. The decay to an n = 1 
pair of doublet leptons (Z2-even) is kinematically allowed, except for very low values of AR < 3. 
The n = 1 states then decay to the corresponding n = leptons, plus the LKP B^\ so that 
the resulting collider signature is a pair of soft leptons (for charged lepton channels) plus a 
huge amount of missing energy (excited neutrinos, of course, will go undetected). These final 
state soft leptons should be detectable [2311 [235J . Similarly, the Z^ can also decay into n = 2 
and n = doublet leptons. Both of these modes are KK-number conserving, but there is an 



4. Universal Extra Dimensions and Colliders 



67 





i 10 



Figure 32: Left: Xyj, the KK number violating couplings, as a function of AR, for = 300 
GeV (the values are independent of R). From top to bottom, the curves correspond to Xzl, 
X ie , XryL, Xzq, Xjd, X yu , and X^q, respectively. Right: Decay widths of (upper pair) 
and i?( 2 ) (lower pair) as a function of AR, for R^ 1 = 450 GeV and 300 GeV (upper and lower 
curves in each pair, respectively). From Ref. [232J. 



important difference: while the coupling is the usual g for the latter channel, it is g/V2 for the 
former ones [232J. 

The also features KK-number violating decay modes, but only to left-handed particles. 
Since the lower limit on R" 1 is about 300 GeV, both of these gauge bosons can decay even to 
a (zero mode) tt pair. 

In the minimal UED model with m| = 0, the Z^ cannot decay through the Bjorken channel 
to ZW-ffW, as it is kinematically forbidden. Moreover, the three-body channels, with a virtual 
ZW or will be even more suppressed. However, if m| < 0, all of the Higgs masses will 

be lowered, and the production of a neutral CP-even Higgs will be possible. The decay pattern 
of the H^ 1 ' is dominantly into right-handed r pairs (assuming the mixing in the n = 1 level 
is small) plus an LKP. If the taus are soft enough, they may escape detection, leading to an 
invisible decay mode. 

The second level KK gauge bosons are produced as s-channel resonances in e + e~ collision 
through KK-number violating couplings. This suppression brings down the peak cross-section 
from an otherwise expected nanobarn level to about 35-45 pb for Z^ and to about 63 pb for 
(for R' 1 = 300 GeV). For Br 1 = 450 GeV, these figures drop down to 16-21 pb and 28 
pb, respectively. The reason for the larger production cross-section is its narrower width 
compared to that of the Z^ 2 \ Shown in the right frame of Fig. [32] are the decay widths of the 
and B^ as a function of AR. 

Since the B^ decays almost entirely to two jets, it will be very challenging to detect at the 
LHC, as the signal would be swamped by the very large QCD background, and because the 
resonance itself is quite narrow. The case of the Z^ is instead more promising, as it features 
various hadronically quiet decay modes, and soft leptons with energy greater than 2 GeV might 
be, in principle, detectable. For a precision study of these resonances, the ILC would be an 
ideal machine. In particular, on the Z^ peak, the ratio, 1Z, of e + e~ to two jets over e + e~ 
to [V fi~ would show a sharp dip. The latter can be even more marked including the missing 
energy events. In fact, the Z( 2 )-width is dominated by decays into n = 1 lepton pairs, and 
quarks can appear only from KK-number violating interactions. On the other hand, 1Z should 
show a sharp peak corresponding to the B^ resonance. Furthermore, the total e + e~ cross- 
section would show a kink between the two peaks, corresponding to the value of where the 
KK-number conserving channels open up. 

With polarized beams, the behavior of the two peaks will also be quite different. Since 
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Figure 33: Projected 95% confidence level exclusion limits for minimal UED in the AR, R^ 1 
plane from combined measurements of leptonic and hadronic final states at an ILC. Deviations 
from the SM expectations can be observed up to the values given. From Ref. 



couples only to the left-handed fermions, with suitable polarization the peak may entirely 
vanish, or it may get enhanced by a factor of 3 (assuming 80% polarization and 60-70% e + 
polarization) |232| . The peak will get enhanced by about a factor of 2 with a left-polarized 
e~ beam, but will never vanish altogether, as the hypercharge gauge boson, B, couples to both 
right and left-handed fermions. 

Although in the minimal UED model the B^ does not have any KK-number conserving 
decay channels, in non-minimal models there are two possibilities for such channels. First, 
non-universal boundary terms can enhance the B^ mass, possibly kinematically opening KK- 
number conserving decay channels. Second, asymmetric boundary terms (different for y = and 
y = ttR) would break KK-parity, resulting in decay modes like B^ — > E^+e. A precision study 
at the ILC would in all cases allow for the discrimination between such different realizations 
of the UED model. As a concluding remark, we point out that the expected background level 
from the continuum at the ILC is less than 10 pb for = 600-900 GeV |237j . and may be 
further reduced by suitable cuts [232J. 

The region of the minimal UED parameter space that can be ruled out at an ILC (at the 
95% confidence level, in the AR, R -1 plane) is shown in Fig. [33j The assumed ILC integrated 
luminosity corresponds to 1 ab _1 with 5Li n t =0.1%. Uncertainties due to the identification of 
leptonic and hadronic final states are assumed to result in systematic errors of 0.1%. In the 
analysis of Ref. [236J a polarization of 80% for the electron beam and 60% for the positron 
beam, and an error of 0.1% for the polarisation measurement were considered. Also, Ref. [236J 
assumes that 50% of the luminosity is spent for each sign combination, (+ , -) and (- , +), of 
the beam polarization. 
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If signals of new physics are detected at the LHC in the form of the production and subsequent 
decay of new heavy particles, a crucial question will be whether it is possible to discriminate 
between different beyond the Standard Model scenarios. As we pointed out before, the phe- 
nomenology of low energy supersymmetry and of UED at colliders are strikingly similar. In 
both models, the lightest new state is neutral and stable, leading to collider signatures includ- 
ing missing transverse momentum (energy) plus a number of jets and leptons. Furthermore, 
both supersymmetry and UED predict the same couplings for the SM particles and their heavier 
counterparts. Distinguishing between supersymmetry, UED and other beyond the SM scenarios 
will likely be a challenging task for the LHC. 

On general grounds, there are a few handles one has which could be used to distinguish 
between supersymmetry and UED. Provided that the LHC measures the mass spectrum of 
the new heavy particles, a quasi-degenerate mass pattern would admittedly favor a UED-like 
scenario, while a split spectrum with sizable mass differences between the heavy states masses 
would certainly be more suggestive of a supersymmetric scenario. Theoretically, in fact, the 
assumption of universality in the soft supersymmetry breaking parameters is thought to hold 
at some high energy scale, e.g. the grand unification scale. As a result, degeneracy at the weak 
scale in supersymmetry is generically rather unnatural. In contrast, the KK structure of the 
UED spectrum predicts an almost completely degenerate spectrum, up to radiative corrections. 

A second distinguishing feature concerns the Higgs structure of the two theories. The particle 
content of the minimal UED setup does not include an analogue of the heavy Higgs bosons of 
the MSSM. More precisely, even though the level-1 KK modes of the Higgs/gauge bosons in 
minimal UED have exactly the same quantum numbers as the MSSM Higgs bosons H°, A 
and H^, the latter do not carry i?-parity. The UED states A 1 ^ and H^ 1 ^, in contrast, 

carry KK parity, and are therefore are more analogous to the higgsino sector of supersymmetry. 
However, the occurrence of additional Higgs bosons is not a robust discriminatory feature of 
supersymmetry, at least at the LHC, as such particles are likely to escape detection over large 
portions of the supersymmetric parameter space. 

One is, therefore, left with the two most striking differences between supersymmetry and 
UED. Firstly, UED predicts the existence of a tower of KK states with almost degenerate masses 
lying around m n ~ n/R for the n-th KK level, rather than the single heavier "copy" of the SM 
particle content found in supersymmetry. Secondly, in UED the spins of the KK states are the 
same as that of their SM counterparts, while in supersymmetry the spins of the superpartners 
differ from their SM counterparts by 1/2. 

In this section, we review the feasibility of discriminating between UED and supersymmetry 
at hadronic (Sec. l5.lj) and leptonic (Sec. 15.2(1 colliders, as well as with astrophysical experiments 
(Sec. 15 .3p . The punchline of this section is that at the LHC, if R^ 1 < 1 TeV, signals of n = 2 
KK states and various spin correlation studies can provide marginal evidence favoring either 
UED or supersymmetry. At a multi-TeV ILC, the two models can be accurately distinguished 
through several processes, including angular distributions of final state leptons and total cross 
section measurements. The particle masses of a UED setup will also be measured to an impres- 
sive degree of precision at such a machine. Finally, direct and indirect dark matter detection 
experiment can, in many cases, provide signatures which would also allow for a discrimina- 
tion between neutralino and KK photon dark matter, or complement and confirm results from 
collider data. 
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5.1 Discrimination of UED and Supersymmetry at Hadron Colliders 

As alluded in the previous section, the detection of second level KK modes would be a provide 
strong evidence in favor of a UED scenario over an accidentally degenerate supersymmetric 
model. The production cross sections and decay patterns of n = 2 KK states at the LHC were 
first addressed in Ref. |238j . It was pointed out that the production cross section of strongly 
interacting n = 2 KK states at the LHC is enhanced, at the same mass scale, with respect 
to the supersymmetric case. This is the case for several reasons. Firstly, the particle content 
in UED is "duplicated" in the KK modes. For example, there are both left-handed and right- 
handed SU(2) doublet KK fermions, while in supersymmetry there are only left-handed-doublet 
squarks. Secondly, the different angular distributions for fermions (1 + cos 2 6) versus scalars 
(1— cos 2 9), when integrated over all angles, account for an extra factor 2. Furthermore, as heavy 
states at the LHC are produced close to threshold, the strong suppression of the cross section 
for the production of scalars in supersymmetry (~ /3 3 ) is replaced by the milder suppression 
for KK quarks (~ (3). However, Ref. |238j points out that any signal from n = 2 KK quarks 
and leptons decaying into n = 1 modes would be likely swamped by the direct production of 
n = 1 KK states. A possible exception might be that of decay chains leading to pairs, 
subsequently decaying into leptons, and in principle giving rise to a very clean and potentially 
observable multilepton signature (Nl+ missing energy, N > 5). The latter would, unfortunately, 
feature very limited statistics, and looks very challenging at the LHC [238J. 

The best prospects for detecting second level KK modes appear to be connected with the 
gauge boson sector. In particular, decays of n = 2 gauge bosons into fermion-antifermion pairs 
can be looked for as a bump in the invariant mass distribution of the decay products, in analogy 
with Z' searches. As discussed in Sec. 14.2.2} the decay widths of KK n = 2 gauge bosons are 
always very suppressed, and much smaller than the typical width of a Z' with SM couplings. 
This implies, experimentally, that the width of the resonance will then be determined by the 
experimental resolution, rather than the intrinsic particle width. 

The production of n = 2 KK gauge bosons proceeds through three main channels: (1) single 
production through KK-number violating operators, (2) indirect production through decays of 
strongly interacting n = 2 states and (3) direct pair production. The first channel was discussed 
in Sec. 14.2.21 and is completely analogous for the case of charged gauge bosons. The production 
of electroweak KK modes in the decays of heavier n = 2 particles is analogous to the dominant 
production of electroweak superpartners in supersymmetry. In particular, the large branching 
fractions of SU(2) doublet, n = 2 KK quarks into and Z^ and of SU(2) singlet, n = 2 KK 
quarks into entail that this channel is a very significant source for the production of n = 2 
gauge bosons at the LHC. Direct pair production through KK-number conserving interactions, 
instead, is kinematically suppressed, since two heavy particles must be produced in the final 
state, and can be neglected for electroweak gauge bosons [238 . 

In analogy to the case of n = 2 KK quarks mentioned above, KK-number conserving decays 
are not very distinctive, since they simply contribute to the inclusive n = 1 sample which is 
dominated by direct n = 1 production. The decays of n = 1 particles will then give relatively 
soft objects, and most of the energy will be lost in the LKP mass. Signatures of second level KK 
modes based on purely KK-number conserving decays are not very promising experimentally 
as the suppressed production cross section for the heavy n = 2 particles is not compensated by 
the benefit of the large mass, since most of the energy is carried away by the invisible LKP. 
KK-number violating channels, in which the n = 2 KK gauge boson decays are fully visible, 
are instead much more promising |238j . In particular, since KK-number conserving decays 
dominate the decays of the n = 2 gluon, searches for the electroweak gauge bosons in leptonic 
channels is predicted to be the best search strategy, even though the branching fraction of 
Z^ 2 \ £?( 2 ) — > l + l~ is always rather suppressed (~2%). 
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Figure 34: 5<r discovery reach for (left) and (right). The y-axis indicates the total 
integrated luminosity L (in fb _1 ) required for a 5a excess of signal over background in the 
dielectron (red, dotted) or dimuon (blue, dashed) channel, as a function of In each plot, 

the upper set of lines labeled "DY" makes use of the single production only, while the lower set 
of lines (labeled "All processes" ) includes indirect production from n = 2 KK quark decays. The 
red dotted line marked "FNAL" in the upper left corner of the left panel reflects the expectations 
for a 72 —> e + e~ discovery at the Tevatron in Run II. The shaded area below R- 1 = 250 GeV 
indicates the region disfavored by precision electroweak data [35J. From Ref. [238J . 



Ref. [238J studied the inclusive production of Z^ 2 \ and looked for a dilepton resonance 
in both the e + e~ and n + n~ channels. As mentioned above, a crucial parameter of the search 
is the width of the reconstructed resonance which, in turn, determines the size of the invariant 
mass window selected by the cuts. Since the intrinsic width of the and B^ resonances 
is so small, the mass window is entirely determined by the mass resolution in the dimuon and 
dielectron channels. For electrons, the resolution in CMS is approximately constant, on the 
order of Am ee /m ee ~ 1%, in the region of interest [238] . On the other hand, the dimuon mass 
resolution is energy dependent and, in preliminary studies based on a full simulation of the 
CMS detector, has been parametrized as [239J 

Am,,,, / m,,,, 
^ = 0.0215 + 0.0128' w 



1 TeV 



The analysis of Ref. |238j used the following cuts: (1) lower cuts on the lepton transverse 
momenta, pr{£) > 20 GeV, (2) a central rapidity cut on the leptons, \rj(£)\ < 2.4, and (3) 
dilepton invariant mass cuts for electrons of 2R~ l — 2Am ee < m ee < my 2 + 2Am ee and, for 
muons, 2R~ l — 2Am w < m MAt < my 2 + 2Am wt . The main SM background to the signal is 
Drell-Yan. 

The discovery reach of the LHC and the Tevatron for the B^ and Z^> resonances, as 
assessed in Ref. [238] . is shown in Fig. [Ml The electron energy resolution assumed for the 
Tevatron was AE/E = 0.01 © 0.1 6/y/E |240j . The Tevatron reach in the dimuon channel is 
worse due to the poorer resolution, and the reach for the Z^ is worse still since m z (2) > m B (2) 
for a fixed R . Even after one year of a low luminosity run, the LHC should be able to test 
for the existence of n = 2 KK states for values of R^ 1 < 750 GeV. Eventually, the LHC should 
be able to probe R^ 1 < 1 TeV. 

As pointed out in Ref. |238j . the discovery of a resonance, while suggestive of a UED setup, 
might be confused with the effect of a Z' gauge boson. A much cleaner signal would be provided 
by multiple quasi-degenerate KK gauge boson resonances, as their is no compelling motivation 
for multiple Z's with approximate mass degeneracy. A double peak structure would be rather 
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Figure 35: The — diresonance structure with = 500 GeV, for the dimuon (left) 
and the dielectron (right) channels at the LHC with L = 100 fb _1 . The SM background is 
shown with the (red) continuous underlying histogram. From Ref. [238J . 



challenging to observe in the dijet channel due to the relatively poor LHC jet energy resolution. 
Again, the dilepton channel from B^ and is the most promising channel. The results for 
the invariant mass distribution of dimuons and dileptons, for R^ 1 = 500 GeV and L = 100 fb _1 , 
and the same cuts described above, is shown in fig. [35j The better mass resolution makes 
the diresonance structure easier to detect in the dielectron channel, although the structure is 
beginning to emerge in the dimuon channel as well. 

The second strong handle on the discrimination between UED and supersymmetry we men- 
tioned at the beginning of this section pertains the determination of the spin of particles at 
the LHC. Spin determinations at hadron colliders appear to be extremely challenging, as the 
center of mass energy of the underlying event is unknown. Furthermore, the momenta of the 
two undetected neutral stable particles in the event are also unknown. 

Ref. [241] suggested that a charge asymmetry in the lepton-jet invariant mass distributions 
from particular cascade decays can be used to discriminate supersymmetry from the case of pure 
phase space decays and is an indirect indication of the superparticle spins. This technique was 
applied in Refs. [242, 238J to the case of the discrimination between UED and supersymmetry. 
Ref. [238] compared the case of UED with = 500 GeV with a supersymmetric model 
featuring a matching particle spectrum and studied the two cascade decays shown in the left 
panel of Fig. [36l The spin of the intermediate particle (Z^ in UED or m supersymmetry) 
governs the shape of the lepton-quark invariant mass distributions, Mi q , for the near lepton. 
In practice, however, one cannot distinguish the near and far lepton, and is forced to include 
the invariant mass combinations with both leptons |241j . This tends to wash out the spin 
correlations but, due to the different quark and anti-quark content of the proton, a residual 
effect remains which, in turn, leads to a difference in the production cross sections for squarks 
and anti-squarks |241j . The spin correlations are encoded in the charge asymmetry [241J 

+ _ = ( dN(qt+) _ dN(q£-) \ I ( dN(qt+) dN(q£~) \ 

{ dM ql dM ql ) J V dM ql dM q i J ' { ' 

where q stands for both a quark and an anti-quark. N{q£ + ) and N(q£~) are the numbers of 
entries with a positively and negatively charged lepton, respectively. The results of Ref. [238] on 
are shown for the UED and supersymmetry cases in the right panel of Fig. [36j Evidently, a 
clean discrimination of the two cases is challenging and probably not completely unambiguous. 
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Figure 36: Left: Twin diagrams in super symmetry and UED. The upper line corresponds to 
the cascade decay q — > qx% ~^ ql ~^ Q^ + £~Xi i n supersymmetry. The lower line corresponds 
to the cascade decay — > qZ^ — > qi^if — > q£ + £~~ji in UED. In either case the observable 
final state is the same: q£ + £~ Ifix- Right: Comparison of the charge asymmetry A" 1 defined in 
eq. ()115p as computed in the case of UED with R~ l = 500 GeV and the case of supersymmetry 
with a matching sparticle spectrum. Adapted from Ref. [238J. 

Further effects (poor statistics near the two ends of the plot, detector effects, backgrounds, etc.) 
can further dilute the cleanness of this discrimination technique. 

Ref. [242J studied, for the same decay chain shown in the left frame of Fig. [36] other spin 
correlations, such as the quark-near lepton mass distribution, the dilepton mass distribution 
and the quark-far lepton mass distribution, and their observable counterparts, arriving at con- 
clusions similar to those of Ref. [238]. In particular, it turns out that it would be much easier 
to rule out a UED setup given a supersymmetry spectrum than vice- versa, as the particle mass 
degeneracy tends to blur the difference in the two cases. 

Recently, Ref. [243] pointed out that a further handle to disentangle UED from supersym- 
metry comes from the determination of the spin of the "gluino" (or KK gluon) from long decay 
chains involving the gluino decay into a sbottom, which decays into the second lightest neu- 
tralino, followed by a slepton, and finally to the LSP (or the corresponding chain of KK modes 
in the case of UED). Using, as an example, the supersymmetric benchmark model SPSla, 
Ref. [243] showed that, using a list of asymmetries constructed from lepton-bottom correlations 
or from pure bottom-bottom correlations, it is possible to distinguish between supersymmetry 
and UED cascade interpretations, and thus possible to identify the fermionic (bosonic) nature 
of the gluino (KK gluon). 

In a scenario where the compactification radius is very small (or, correspondingly, in 
the case of supersymmetry, where the mass of sleptons is small) an observable that could help 
unraveling the spin of KK leptons (sleptons) was proposed in [244] . The analysis considers the 
processes 

qq^Z°/ 7 /+ l~ (116) 

gg->Z°/ 7 - r+T--> X\l + Xir (117) 

and pointed out that since the leptons (e, n) are highly relativistic, their pseudo-rapidities are 
to a very good approximation close to their true rapidities. In particular, making use of a 
function of the pseudopodia difference between the final state leptons, Arji+i-, one no longer 
has to face the problem of determining the center-of-mass frame along the beam direction, still 
being sensitive to the KK leptons (sleptons) angular distribution, that, in turn, depends upon 
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the spin of the particles. Ref. [244] proposed to use the angular variable 



cos fy* = cos (2 tan 1 exp (Ar/i+i- /2)) ; 



(118) 



that corresponds to the cosine of the angle between each lepton and the beam axis in the 
longitudinally boosted frame where the (pseudo-)rapidities of the final state leptons are equal 
and opposite, cos 9f, features the virtue of being longitudinally boost invariant, and is on 
average smaller for SUSY than for UED, providing a potentially spin-sensitive discriminant in 
KK leptons/sleptons pair production at hadron colliders. Ref. |244j showed that with proper 
cuts it would be possible to statistically distinguish between scalar and fermionic heavy leptons. 
In particular, the anticipated LHC reach with 100 (300) fb -1 integrated luminosity extends to 
202 (338) GeV for left-handed sleptons, and to 143 (252) GeV for right-handed sleptons [233] . 

5.2 Discrimination of UED and Supersymmetry at a Multi-TeV e + e~ Col- 



Assuming the LHC observes signals of new physics consistent with either n = 1 KK modes in 
UED or sparticles in supersymmetry, Ref. |235j studied the discrimination of these scenarios at 
a post-LHC facility, namely a linear collider tunable over a center-of-mass energy range between 
1 TeV and 3 TeV. Ref |235j concentrated on a minimal UED model featuring R^ 1 = 500 GeV 
and AR = 20, and adjusted the super symmetric low-energy scale parameters so that the two 
smuon masses and the neutralino mass matched the masses of their KK counterparts (the KK 
muons and the B^ 1 '). The most promising process at an ILC with large center-of-mass energy 
is found to be the pair production of KK muons (smuons), which then decay into a muon 
pair and two B^'s (neutralinos). While, in supersymmetry, the production process of smuon 
pairs proceeds through 7 and Z s-channel exchanges, in UED, KK muons can be also be pair 
produced via s-channel n = 2 KK electroweak gauge boson exchange, a process that can occur 
close to resonance. 

The signal cross section for a center-of-mass energy of 3 TeV is 14.4 fb for the case of 
UED and 2.76 fb for supersymmetry. The SM backgrounds come from n + [i~vv final states, 
mostly due to gauge boson pair production W + W~ — > /i + [i~ v ^, Z°Z° — > /i + [i~vv and from 
e + e _ — > W + W~v e v e , e + e~ — ► Z°Z°f e P e , followed by muonic decays. The total background 
cross section is ~20 fb at yfs = 3 TeV. The background muons are, furthermore, produced 
at small polar angles, therefore potentially biasing the angular distribution. To reduce this 
background, Ref [235J imposed event selection cuts including the requirement of detecting two 
muons, missing energy in excess to 2.5 TeV, transverse energy below 150 GeV and event spheric- 
ity larger than 0.05. The rejection of the background was performed discarding events with 
di-lepton invariant mass compatible with M z o. A further source of background, 77 — * fJ- + [i~ 
from underlying 77 collisions, can be completely suppressed by a cut on the missing transverse 
energy, fi™ lssm 9 > 50 Q e y Finally, events with large beamstrahlung were cut imposing an 
event sphericity smaller than 0.35 and acolinearity smaller than 0.8. The combination of these 
criteria provide a factor of ~ 30 background suppression in the kinematic region of interest, 
while not significantly biasing the lepton momentum distribution [235J . 

A first distinguishing feature in the KK muons / smuons pair production is the angular distri- 
bution of events as a function of the muon scattering angle, 6^. As long as the mass differences 
iW^(i) — M S (i) and M^ — M^o, respectively, remain small, the muon directions are well correlated 
with those of their parents, for which the angular distributions are given by 



lider 




UED 



1 + COS 



(119) 
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Figure 37: Left: The initial state radiation corrected differential cross-section da/dcos9 ll for 
UED (blue, top) and supersymmetry (red, bottom) as a function of the muon scattering angle 
9^. Right: The same, including the effects of event selection, beamstrahlung and detector 
resolution and acceptance , for UED (left) and supersymmetry (right). The data points are 
the combined signal and background events, while the yellow-shaded histogram represents the 
signal only. From Ref. [235J. 



in the UED case, and by 

da 



dcos 9 



~ 1 - cos 2 9. (120) 

SUSY 



in the case of supersymmetry, assuming production occurs well above threshold. The two an- 
gular distributions can be clearly distinguished even after accounting for initial state radiation, 
event selection, beamstrahlung and detector resolution and acceptance, as shown in Fig. [37l A 
X 2 fit to the normalized polar angle distribution shows that UED and supersymmetry can be 
distinguished on the sole basis of the distribution shape, with as little as 350 fb _1 of data at 
= 3 TeV. 

At an e + e~ linear collider, the muon excitation masses can be accurately determined through 
an energy scan across the onset of the pair production threshold. Such a scan, impossible at a 
machine like the LHC (recall that in a hadron the center of mass energy of the partons involved 
in the underlying event producing the new physics massive particle pairs is unknown) will also 
reveal the spin of the particle, as the cross sections for the UED processes rise at threshold oc (3 
while in supersymmetry their threshold onset is oc /? 3 , where f3 is the particle velocity. From the 
estimated sensitivity, da/dM, and the cross section accuracy, Ref. |235j claims that the masses 
of the two UED muon excitations can be determined to ±0.11 GeV and ±0.23 GeV for the 
SU(2) singlet and the doublet states, respectively, with a total luminosity of 1 ab _1 shared in 
three values of the center-of-mass energy, where the particle widths can be disregarded. Fig. [38] 
illustrates the accuracy at which the ILC will be able to tell apart singlet or doublet KK muons 
with mass differences of 2.5 GeV. 

A further handle that will be completely under control at the ILC, unlike the LHC, is the 
determination of the total production cross section for KK muons or smuons as a function of 
the center of mass energy. As we stated above, at a center-of-mass energy of 3 TeV, the two 
cross sections differ by a factor 5, to be contrasted with a projected CLIC absolute luminosity 
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Figure 38: Left: A threshold scan at selected points. The green curve refers to the reference UED 
parameters while for the red (blue) curve the mass of (j&p (a*d ) nas been lowered by 2.5 GeV. 
The points indicate the expected statistical accuracy for the cross section determination at the 
points of maximum mass sensitivity. Right: The muon energy spectrum resulting from KK 
muon production (left) and smuon production (right), including the effects of event selection, 
beamstrahlung and detector resolution and acceptance. The data points are the combined signal 
and background events, while the yellow-shaded histogram is the signal only. From Ref. [235J. 



determination that should be measurable to 0(0.1 %) and the average effective collision energy 
to 0(0.01 %). 

The mass determination also benefits from the study of the characteristic end-points of the 
muon energy spectrum (see Fig. \38\ right panel) . Even though these do not depend on the spin 
of the decaying parent particles, and will therefore not directly distinguish between UED and 
supersymmetry, their study can lead to an impressively accurate mass determination, to the 
level of ±0.19 (stat.) ±0.21 (syst) GeV, where the statistical uncertainty is given for 1 ab _1 of 
data and the systematics reflect the uncertainty in the fii masses. 

In addition to the possibilities discussed above, the distinction between UED and super- 
symmetry can proceed through the study of the sharp peak in the photon energy spectrum 
predicted when one of the mediating s-channel particles is on-shell; in that case, since the decay 
allowed by phase space, there will be a sharp peak in the photon spectrum, 
due to radiative return to the Z^ 2 \ In minimal UED, the peak occurs only for the Z^ 2 ' and not 
for the B( 2 \ since m B (2) /2 < m„(i) , and gives rise to a peak located at a photon energy 

1 ( Ml \ 

E 7 = -E CM (l-^J- (121) 

The same peak cannot occur in supersymmetry, where smuon production is mediated only by 
the Z and the photon, and both particles are always far from being on-shell (recall that the 
lower limit on the smuon mass from LEP is roughly > 100 GeV S> m^/2). 

Ref. [235J also considered the prospects for the pair production of KK leptons and quarks. 
In the case of KK tau pairs, the lower statistics and the inferior jet energy resolution make the 
channel less competitive than the muon channel. For electrons, instead, the extra t-channel 
diagram occurring in the pair production of KK electrons enhances the total production cross 
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section and can be up to two orders of magnitude larger than in the case of muons. However, 
despite the larger event sample, the additional t-channel diagram also distorts the differen- 
tial angular distributions discussed above, creating a forward peak, which causes the cases of 
UED and supersymmetry to look very much alike. Finally, in the case of KK quarks, the jet 
angular distribution will again be indicative of the KK quark spin, and can be used to dis- 
criminate against (right-handed) squark production in supersymmetry. Furthermore, the jet 
energy distribution will again exhibit endpoints which will, in principle, allow for precise mass 
measurements. The jet angular and energy measurements, however, would not be as clean as 
in the case of lepton (muon or electron) final states, which would therefore provide the most 
convincing evidence for UED/supersymmetry discrimination. 

5.3 Discrimination of UED and Supersymmetry With Dark Matter Exper- 
iments 

In addition to using collider experiments, it may be possible to distinguish UED from super- 
symmetry with direct and indirect dark matter detection experiments. Following Ref. [245J , we 
briefly summarize the prospects for this here. 

The spin-independent (scalar) elastic scattering cross section for the with nuclei is quite 
small compared to that found for neutralinos in many supersymmetric models. In particular, the 
WIMP-WIMP-Higgs coupling in supersymmetry is typically, although not always, considerably 
larger than in UED (in part because of the heavy CP-even Higgs in supersymmetry). That 
being said, the elastic scattering cross section (as determined by the direct detection rate) will, 
alone, likely be insufficient to discriminate between these models. 

In contrast, the spin-dependent (axial-vector) elastic scattering cross section for the B^ 
with protons is quite large, & B (i) pSD 1.8 x 10~ 6 pb for a 1 TeV mass and a 10% LKP-KK 
quark mass splitting (see Sec. I3.3.2[) . Neutralinos can possess a similarly large spin-dependent 
cross section, but only if their higgsino fraction is substantial (> 1%), which would also imply 
a fairly large spin-independent cross section [159J. 

The discriminating power of these elastic scattering cross sections come from the ratio of the 
two types. A comparison of the spin-independent cross section, as manifest in direct detection 
experiments, with the spin-dependent cross section, as manifest in neutrino fluxes from the Sun, 
can act as a powerful discriminator of supersymmetry and UED. 

In Fig. 1391 a comparison of these two rates is shown for neutralino and B^ dark matter (see 
Ref. [245J for details on the adopted scan over supersymmetric model parameters). The range 
of supersymmetry and UED modes shown demonstrates that, although neither the direct nor 
the neutrino rate alone will be sufficient to distinguish between these models, the combination 
of the two rates can be very powerful. For a particular value of the spin-independent elastic 
scattering, UED generally produces between one and three orders of magnitude more neutrino 
events. If these rates can each be measured in future experiments, such a comparison would be 
expected to be very useful in discriminating between UED and supersymmetry. 

Other indirect detection channels have the potential to distinguish between UED and super- 
symmetry by observing distinctive features in the WIMP annihilation spectrum. In particular, 
the positron and gamma-ray spectra each have considerably larger fluxes at high energies in 
the case of UED than is predicted for neutralino dark matter. This is a result of the chirality 
suppression of neutralino annihilations at low velocities to light fermions. B^'s, on the other 
hand, annihilate a large fraction of the time to light fermion pairs. 

In the case of the cosmic positron spectrum (see Sec. 13.6.31) . the spectrum is hardened in UED 
as a result of the large fraction of the annihilations to e + e~ (20-25%) and, to a lesser extent, to 
(20-25%) and t + t~ (20-25%). The spectrum resulting from these annihilation modes is 
shown in Fig. l4"0l where it is compared to the spectra representative of neutralino dark matter 
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Figure 39: A comparison of the neutrino and direct detection rates for neutralinos and Kaluza- 
Klein dark matter. The mass of the WIMP in each case shown is in the range of 450-550 GeV. 
The shaded region represents the range of supersymmetric models which satisfy all collider 
constraints, current direct detection constraints, and do not overproduce the abundance of 
neutralino dark matter. The thick solid lines denote the UED case varied across a broad range 
of KK quark masses (m (i) = 1.05 m B (i) to 1.3m B (i)). The two UED contours represent the 
maximum and minimum Higgs masses consistent with electroweak precision observables. From 
Ref. [2i5] . 

(annihilations to bb and IU + IU _ ). As neutralinos annihilate almost entirely to combinations of 
heavy fermions, gauge and Higgs bosons, positrons are largely generated through relatively soft 
channels. Even in the case of neutralino annihilations to gauge bosons, the positron spectrum 
is considerably softer at high energies than is found in UED. 

The gamma-ray spectrum produced in dark matter annihilations in UED is also harder 
at high energies than in the case of super symmetry. This is a result of the large fraction 
of annihilations to e + e~ (through final state radiation) and t + t~. The resulting gamma-ray 
spectrum from annihilations is shown in Fig. [T7] of Sec. 13.41 The spectrum predicted for 
neutralino dark matter, in contrast, is produced through the fragmentation and cascades of 
heavy quarks and other annihilation products, and is described by the softer dashed line in 
that figure. If future gamma-ray telescopes were able to measure the dark matter annihilation 
spectrum with sufficient detail, it would be possible to use this information to distinguish 
between UED and supersymmetry. In some supersymmetric setups, however, the effects of final 
state radiation can become non-negligible, leading to a considerably harder spectrum than what 
expected from fragmentation alone [246, 247J , blurring, to some extent, the possibility of a clear 
distinction between UED and supersymmetry. 

The prospects for each of the methods described in this section depend strongly on the 
precise nature of the WIMP being studied and on a number of astrophysical uncertainties. 
If the -B^ 1 ) is fairly heavy and it has a fairly large mass splitting between it and the KK 
quarks, for example, it will be difficult to observe at future direct detection experiments and 
impossible to detect with planned neutrino telescopes. On the other hand, a lighter UED model 
spectrum, with a small LKP-KK quark mass splitting, could be effectively studied (and likely 
distinguished from a neutralino) by using these techniques. The prospects for anti-matter and 
gamma-ray observations depend on the annihilation rate of dark matter in the galactic halo 
and other astrophysical environments, making the perspectives of detecting a signature with 
such techniques difficult to assess. 

The low-energy phenomenology of "hybrid" effective theory scenarios where a supersym- 
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Figure 40: The ratio of positrons to positrons+electrons as a function of energy, after diffusion 
through the galactic halo, for a 600 GeV WIMPs annihilating to bb (dashes) to gauge bosons 
(dot-dash) and to the modes of Kaluza-Klein (UED) dark matter, which are primarily charged 
leptons (dotted). The errors bars correspond to the measurements of the HEAT experiment 
|171j . To normalize the annihilation rate, a boost factor of approximately 10 has been used. 
The lower solid line represents the background from the secondary production of cosmic ray 
positrons. 



metric setup is embedded in an extra-dimensional framework, and where some or all of the SM 
fields are allowed to propagate in the extra dimensions was also recently addressed from the 
perspective of dark matter physics in Ref. [248]. As the compactification scale in these models 
is much larger than the soft-super symmetry breaking scale, the resulting dark matter candidate 
is the neutralino, and the related phenomenology is very similar to that of supersymmetry in 4 
dimensions [248J . 
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6 Beyond One Extra Dimension 

As compared to the numerous papers devoted to five-dimensional UED models, the phenomenol- 
ogy of models with more than one universal extra dimension has been somewhat rarely discussed 
in the literature so far. In some cases, models with d > 2 have stronger theoretical motivations 
than the d = 1 case within the context of UED (see Sec. 12. 1|) . For instance, the d = 2 case is 
motivated by the cancellation of the global SU(2) anomaly [249] and by the fact that the sim- 
plest chiral compactification of two universal extra dimensions, the "chiral square", preserves 
a discrete symmetry, a subgroup of the six dimensional Lorentz group, which suppresses the 
proton decay rate even when baryon number is maximally violated at the TeV scale [18] (see 
Sees. EI]). 

The Feynman rules for gauge theories in d = 2 UED compactified on the chiral square are 
given in Ref. |250j . The KK excitations are indicated with two KK numbers (j,k) which are 
manifest in the conservation of the quantity (— l) J+fc , somewhat resembling the KK parity of five 
dimensional UED. The KK gauge bosons in this model feature two scalar fields corresponding 
to the polarizations along the two extra dimensions. At every KK level, these two scalars mix 
into a linear combination that is eaten by the spin-1 KK mode, and another which remains 
as a physical real scalar field, referred to as the "spinless adjoint" |19j . As a consequence of 
the presence of the spinless adjoints, the production cross sections and decay modes of other 
KK states are different than are found in the d = 1 case. In particular, the (1,1) states, which 
feature a mass of only y/2/R and can be produced in the s-channel, are particularly relevant to 
collider phenomenology [25H 12501 1252] . 

Several aspects of the phenomenology of six dimensional UED models were recently explored. 
Ref. [253J addressed the questions of neutrino masses and nucleon stability, and studied models 
based on the gauge group SU{2) L x U{l) hR x U(1) B -l compactified on T 2 /Z 2 and left-right 
symmetric extensions featuring a full SU(2)l x SU(2)r x U(1)b-l gauge group, compactified 
on the orbifold T 2 jZ^ x Zj. In both cases, neutrino masses are suppressed by an appropriate 
orbifold parity assignment for the standard model singlet neutrinos, and the proton decay rate is 
suppressed due to a residual discrete symmetry, left over from compactification. An interesting 
phenomenological outcome of these models is that, for low values of the fundamental scale, a 
dominant decay mode of the neutron is three neutrinos. Ref. [253] also pointed out that the 
model possibly container's a two-component dark matter setup, where the stable relics are KK 
right handed neutrinos and photons. 

Ref. [254J elaborated on this latter point, and showed that the left-right symmetric six 
dimensional UED model of Ref. [253] indeed thermally produces the desired amount of two- 
component cold dark matter, for values of the compactification scale between 400 and 650 
GeV, and a mass for the heavy extra gauge boson Z 1 below 1.5 TeV. Ref. [254J finds that a 
generic prediction of this model is a dark matter-nucleon scattering cross section typically in 
excess of 10 -44 cm 2 , accessible to next-generation direct dark matter searches. Furthermore, 
the parameters of the model giving rise to the appropriate amount of thermal relics are such 
that the theory can be fully tested at the LHC. 

In Ref. [255] the two-component dark matter scenario was further investigated, including 
the effect of coannihilations. Also, the possibility that the spinless adjoint photon is lighter than 
the KK photon led the authors to explore this alternate scenario, featuring a the physical scalar 
gauge boson as one of the two dark matter constituents. The range of parameters is similar 
to the one outlined in Ref. [254J , and a discussion of the extra Wr and Z' phenomenology is 
provided as well. 

In Ref. [256J dark matter in a six-dimensional UED model was studied. In this scenario, 
the lightest Kaluza-Klein excitation in this d = 2 model is a neutral scalar, corresponding to 
a six-dimensional photon polarized along the extra dimensions. Annihilations of this "spinless 
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photon", as it is known, proceed largely through s-channel Higgs exchange to gauge bosons, 
Higgs bosons and, to a lesser extent, top quarks. Unlike in the d = 1 case, annihilations of the 
LKP to heavy fermions are chirality suppressed. 

The thermal relic abundance of spinless photons can be consistent with with measured 
cosmological dark matter abundance for masses of approximately 500 GeV or less. The elastic 
scattering cross sections of this dark matter candidate with nuclei are highly suppressed, leading 
to very small rates in direct detection experiments and in neutrino telescopes. 

Many of the phenomenological characteristics of this dark matter candidate in d = 2 UED 
resemble those of a super symmetric neutralino than those of KKDM in d = 1 models. In 
particular, the dominant annihilation modes, indirect (gamma ray and antimatter) detection 
prospects and the preferred mass range are quite similar to those of neutralinos found in typical 
supersymmetric scenarios. 
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7 Conclusions 

In the opinion of the authors, the universal extra dimensional model reviewed here remains 
an extremely fertile field for further theoretical and phenomenological investigation. While the 
UED model possesses several attractive features, numerous aspects of this setup still need to 
be explored or clarified. Not expecting to be exhaustive, we list below those issues we consider 
of particular relevance: 

• As the UED model is one of the most promising new physics scenarios accessible to 
the Tevatron Run-II, we recommend that the CDF and DO collaborations undertake a 
careful data analysis to search for signatures of KK states, extending the analysis of 
Ref. [40j not only to a much larger data set, but also to other search channels, and to UED 
realizations beyond the minimal case. At the same time, it would be useful if dedicated 
collider simulations exploring non-minimal UED setups and novel search channels became 
available. 

• In preparation for the beginning of operation at the LHC, we feel that the issue of dis- 
criminating among different TeV-scale physics scenarios should be more extensively inves- 
tigated, possibly contrasting UED and weak-scale supersymmetry with other beyond the 
Standard Model scenarios. Along these lines, further studies on the spin reconstruction 
of new particles to be produced at the LHC would certainly be particularly interesting. 
Furthermore, it would be highly desirable if collider simulation numerical codes were 
upgraded to include the UED scenario as an option. 

• The realm of UED models beyond one extra dimension certainly deserves the attention of 
both theorists and phenomenologists. For instance, the collider implications and the dark 
matter physics of d > 2 UED models have not been so far thoroughly investigated (with 
the exceptions of the studies mentioned in Sec. [6]), nor was the resulting phenomenology 
systematically contrasted to the five dimensional case. 

• The cosmological and phenomenological issues connected to KK gravitons and to radion 
(gravi-scalar) fields seriously jeopardize the UED model: for instance, the large portion 
of the parameter space of the minimal UED scenario where the n = 1 KK graviton is the 
LKP, as explained in Sec. 13.71 is ruled out. We think that further investigations of the 
physics of KK gravitons in UED and, more in general, of the cosmology of models with 
UED should be undertaken. 

• It would be highly desirable to encompass the UED framework in a theoretically motivated 
high-energy completion, and investigate the candidate theories that might lie beyond the 
UED effective field theory setup. 

• UED triggered a renewed interest in dark matter candidates within the context of extra- 
dimensional scenarios. We think that this line of investigation should be further pursued 
starting, for example, with non-minimal or higher dimensional UED constructions. 

There is a thought that one invariably has after spending many months writing a review 
of this nature: How long will it be before this article is outdated, or even obsolete? Strangely 
enough, we optimistically think that it will not be long at all. 

With the LHC about to begin operation, and numerous dark matter experiments being 
planned and currently underway, it seems very likely that our understanding of the TeV-scale 
will change dramatically in the near future. As this transition takes place, there will be a rare 
and exciting opportunity to witness a new energy scale of physics. The origin of electroweak 
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symmetry breaking, the solution of the hierarchy problem, and many other long shrouded 
mysteries will soon be uncloaked. 

Whether beneath that cloak lies universal extra dimensions or other new physics, it will be 
an exciting time to do physics. 
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A Benchmark Models 

We devote Appendix [A] to the definition of the four benchmark models used throughout this 
review, UED1-4. For all models we assume vanishing Higgs boundary mass term (fhh in the 
notation of Ref. |30j). All of the models produce a thermal relic abundance of B^s which is 
roughly consistent with the cold dark matter abundance as inferred in Ref. [72] for a standard 
ACDM cosmology. 

Model UED1 is motivated by the results of Ref. [73], and assumes a particle spectrum 
given by the radiative mass corrections computed in Ref. [30J, and takes into full account 
coannihilation processes. Second KK level resonances are here irrelevant [78]. It features 1/R = 
550 GeV, AR = 20 and = 120 GeV. In principle, the LKP would be, for this choice, the 
KK graviton, but we assume that somehow the KK graviton mass is lifted above the B^ mass 
and that it doesn't play any role in affecting the relic abundance of LKPs in the Early Universe 

M- 

Model UED2 is motivated by the results of Ref. [77], where the effect of taking a smaller 
value for AR was evaluated, together with full coannihilation processes. We take here 1/R = 850 
GeV, AR = 4 and again m h = 120 GeV. 

Model UED3 is motivated by the results of Ref. [75] , in particular on the effect of the second 
KK level resonances when the Higgs mass is large. It features 1 / R = 1000 GeV, AR = 20 and 
m h = 220 GeV. 

Finally, for model UED4, we neglect the detailed structure of the UED spectrum including 
radiative corrections, and simply assume that all first level KK states, except for the LKP, lie 
at a mass splitting A = 0.02 above the mass of the B^ ~ 1/R = 2250 GeV (i.e. at a mass of 
2295 GeV). The choice of the value of 1/R is here motivated by the analogous models studied 
in Ref. [73]. In principle, this simplified spectrum might occur if suitable boundary kinetic 
terms occur at the cutoff scale, A. The main motivation for considering this case study is to 
understand what occurs when the mass range of the KK particles is in the multi-TeV range 
and, yet, the LKP relic abundance is compatible with the measured density of dark matter. 

We collect in Tab. [3] the details of the mass spectra for the models UED1, UED2 and 
UED3 (the spectrum of model UED4 is trivial, as explained above). All masses are given in 
GeV. 
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Table 4: The input parameters and the mass spectrum for models UED1, UED2 and UED3. 
All quantities with the dimensions of a mass are in units of GeV. Model UED4 features m B (i) = 
2250 GeV, and all the other KK particle masses are set to 2295 GeV. 
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B The UED Feynman Rules and Cross Sections 

We collect in this appendix the relevant Feynman rules (Sec. IB. 1 [) . all of the annihilation and 
coannihilation cross sections relevant for the computation of the relic abundance of the lightest 
Kaluza-Klein particle in the early universe (Sec. IB.2|) . and the amplitudes squared for the 
production of strongly interacting Kaluza-Klein states at hadron colliders (Sec. IB.3[) 



B.l The UED Feynman Rules 

We indicate here only the interactions between SM particles and level-one KK modes allowed 
by KK-parity. Following Ref. [77], we choose to work in the unitary gauge, where the Goldstone 
bosons appear as the longitudinal polarizations of the massive KK gauge bosons, rather than 
as external particles. 

The mass eigenstates, f^ n \ of the n-th KK fermion modes are in general linear combinations 
of the weak eigenstates, /'("). For instance, in the up-type quark sector, one has 



'(") \ / (n) 

- \ / — 75 cos a y ' 



V Q't n) ) V 75 sin a W 
where, neglecting radiative corrections [30] . 




(122) 



tan2o» = for n > 1. (123) 

n/R 

Fermion- Fermion- Gauge Boson vertices. In what follows we consider the vertices involving 
the mass eigenstates. For brevity, however, we will suppress the prime to the fermionic fields 
introduced in Eq. (|122p . Also, we explicitly show only the case of quarks. The Feynman rules for 
leptons can be trivially obtained replacing the appropriate charges. We indicate with s u ^ and 
c Ut d the sine and cosine of the level 1 fermion mixing angles, defined in the previous paragraph. 
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The coefficients Cf 1 f 2 z- 



The other vertices involving two fermions and a gauge boson follow directly from the SM 
vertices in the weak eigenstate basis. To get the right vertices it is therefore sufficient to perform 
the rotation from the weak to the mass eigenstate basis. The same applies for all of the vertices 
involving a gluon or a KK gluon. 
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Fermion-Fermion- Scalar vertices. We include here the non-trivial vertices involving the 
mass eigenstates resulting from the superposition of the charged and neutral KK Higgs bosons 
(or rather, of the zero-level Goldstone bosons which give masses to the Z and the W) and 
of the gauge bosons. We indicate the physical charged and neutral states by H^ 1 ' zk and 
respectively. We neglect here the effect of radiative corrections to the masses. The vertex 
involving the charged KK "Higgs" can be represented as in the figures. The vertices involving 
a scalar KK "Higgs", together with the relevant Cf" R coefficients are also given below. Again, 
the leptonic vertices follow from those reported above via a simple substitution of the relevant 
charges. The vertices involving the KK first-level Higgs are the same as those in the SM once 
the proper rotation of the weak eigenstates into the mass eigenstates is performed. 
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Gauge Boson-Gauge Boson-Scalar vertices. The table below shows the vertices involving 
two gauge bosons and one scalar. We use the sign convention of Ref. [77] and adopt the factor, 
(— which accounts for the e abc from the commutators of the SU(2)\y generators. (3 is even for 
vertices where the electric charge entering the vertex increases in the counterclockwise direction, 
and is odd otherwise. We do not include explicitly H^A^'^A^' , where A = W ± , B,W 3 , 
as they follow exactly from those in the Standard Model. Also, for notational ease (to avoid 
double superscripts), though with a slight abuse of notation, in this appendix we indicate the 
neutral component of the SU(2)w gauge bosons, W 3 , with the symbol Z. 
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B.l. The UED Feynman Rules 



Gauge Boson-Scalar-Scalar vertices. The table below shows the vertices involving two first 
KK level mode scalars and a gauge boson. 
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Gauge Boson-Gauge Boson- Scalar- Scalar vertices. We show in the following table the 
vertices involving a pair of zero-mode gauge bosons and a pair of first KK level mode scalars. 
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Three-Scalar vertices. The following table collects the vertices featuring three scalars, of 
which two are KK first level modes and one is a zero-mode. 



aWaMh 









Four-Scalar vertices. We collect in the following table the vertices involving four scalars, of 
which two are first KK level modes and two are zero-modes. 
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B.2 Annihilation and Coannihilation Cross Sections 

We review here the annihilation and coannihilation cross sections relevant for the computation 
of the LKP relic abundance [25l EH EOl EH EHJ [78] . We collect the various (co-)annihilation 
processes below by the type of (co-) annihilating particles (gauge bosons, fermions and Higgs 
bosons). We neglect in what follows the mass splitting between particles belonging to the same 
KK level, and call their common mass m = 1/R, as well as EWSB. <p indicates generically 
the SM charged and neutral Higgs bosons, and G indicates the zero mode Goldstone bosons. 
Following Ref. [73] . we indicate s as the center-of-mass energy squared, with 



/? = \/l-— (124) 



and with 



L = log (irf ) = - 2tanh-1 ^ ( 125 ) 

The other symbols we employ in the cross sections will be defined later. We also list the relevant 
resonant annihilation cross sections as computed in Refs. [T9" l IBOl 178] . 

I. Gauge Boson-Gauge Boson 
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II. Gauge Boson and Fermion 
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.(w'-^i^^^ (195) 

where Z and Z' belong to different lepton families. 

9 4 Cll 2 «(» + 2m 2 ) 



„ (fl*^ - <zz + z-, + „,) . WWty-w^) (197) 



/ fii— fiT -\ N c g 4 (Y 2 + V|) (a + 2m 2 ) , 



where F = T — 



V ' V 7 24m (a-m| (2) ) +m| (2) r| (2) 



where 



and 



r -- = ^f^ ln2(A2i?2) ' (200) 



2 \ 3/2 



r — r -- +r --- + V[ 1 -^ (201) 
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where 



(203) 



( (1)11) „„\ , 77x7 , A g 4 ((8m 4 - 4m 2 s - s 2 )L - fjajfi + 4m 2 )) 



128tts 3 /3 2 c4 

(204) 



ff 4 (/3s - 2m 2 L) g 4 (/3(s + 4m 2 ) + (s + 2m 2 )L) 



+ - 



32tts 2 /3 2 32tts 2 /? 2 



„ /\,UUi) _ v v \ _ „ (,±W,±W ^ ,±,±\ _ ff 4 (/3 g (2 S -m 2 )-m 2 (4 S -5m 2 )L) 



where §l — g/(2c w ) for neutrinos and ()l = g/V% for charged leptons 

til 
24ns 2 (3 



where g 2 = g 2 (Y f Y ln t 2 w + TfT* ) 



/,+(i),-(i) - ,+,-\ gjg 2 (5/3a + 2(2s + 3m 2 )£) 



(205) 

(206) 



c («' - «-) - , (W - If) - * - «0 - (207, 

(d)- ,_,+\ _ g 4 yLg|(5/3 S + 2(2.s + 3m 2 )£) g 4 (F 4 + Y^js + 2m 2 ) 
" l ^ J ~ 647r S 2 /3 2 c 4 + 967r S 2 /3 c 4 + 

, g 4 (/?(4 s + 9m 2 )+8m 2 L) 



327TS 2 /? 2 

gj (/?(4s + 9m 2 ) + 8m 2 £) ~g\ s + 2m 2 ) 

647rm 2 s/3 2 24tts 2 /3 1 j 
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where T q q is given in Eq. (1200|) , and where 



r , v w^ ,9 2 ™l 4m ff> \ 3g 2 m ( m i™ 



2 



where 



1 ijH.HA - ^ l zl0 J 



^ ^ -H = i&?r (217) 



/,-(!)— T/t ™^ A g 4 ^ ((8m 4 - 4m 2 s - S 2 )£/? S ( S + 4m 2 )) 5g 4 (s + 2m 2 ) 
^ ^ (Z + 7) J = 32^2 A^s 2 B + 



+ 



32tts 3 ^ 2 48tts 2 /3 
g 4 (/3s - 2m 2 L) g 4 (/3(s + 4m 2 ) + (a + 2m 2 )L) g 4 m 2 L 
32tts 2 /3 2 + 64tts 2 /3 + 16tts 2 

(218) 



\_ 9 4 tw (50s + 2(2s + 3m 2 )L) gH 2 w (/3(4s + 9m 2 ) + 8m 2 i) g 4 (s + 2m 2 ) 

CT V L ~* n ) ~ 647rs 2 /3 2 (2s 2 „ - 1) + 647rm 2 s/3 2 (24 - l) 2 + 96tts 2 ^ 

(219) 

/ (!)(!) \_ gX(-2m 2 (4 S -5m 2 )£ + m 2 /? S ) Bs(As-'im 2 ) f gHj g 4 \ 

(J \ L Vl ^ n J 64^mV/3 2 (24 - 1) 647rm 2 s 2 /? 2 ^(24 - l) 2 4 ) 

(220) 

V / 2567rm 2 s/3 V ; 

ff fc' (1 M 1} - rv ; ) = g 4 ^-3m 2 ) 

V L ' V 647rm 2 s/3(24 - l) 2 V ' 

\ ^J- 647rm 2 s /3 2 (24 - 1) 2 [226) 



a 



( (D-mT \ 5 4 (0(4s + 9m 2 ) + 8m 2 L) 

~+ » lW ) = 256,m4 2 ^ <™ 



n (iW) -> /A - ff 4 fe 2 ^(8m 2 £ + /3(9m 2 + 4 S )) foo^ 

r m _ _x g Ht Y 2 Y 2 Us - 3m 2 ) , x 

*(1?IV*^U) = 9 - j (226) 
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Below we collect the KK quark-KK lepton pair annihilation cross sections: 



o („(m> nu) = (4g2 ^ r ^s;^tr " 3m2) (228) 



, ^ _ a) = (4^^ +j ^ + ^ + 4.)) (231) 

. (V 1 ^ - „ u\ W*° Y ^ + g2 ) 2 ( 8m2£ + + 45 )) r2321 

a \v t u L - ^uj - T024W^ (2d2J 

a (i' 1 '^ -^lu)=a - iu) = a - Zu) = <7 - Zz) (233) 

(7 (4 1} «F -» i«) = * ~> ««) = ° ( l R )u R } -» i«) = (7 (z^ZMfl -» il) (234) 

The annihilation cross section involving the annihilation of KK quarks through strong interactions are 
collected below: 

^-^') = 9l & <*»> 

' (««V> - «) = ^(^-W^-^)) (236) 

. - = gl (4m2(4S - " m \ + JS% + 33Sm2 + 12S2)) (23?) 

. - = -^(^ 4 + W+^+^ 2 + 7 a) ) (238) 



cr (V 1 )^ - M ') = a (ugV^ ud) = a - «d) = a (u^u^ - uu) 

_ g|(4 S -3m 2 ) 
2887rm 2 s/3 

ff | (8m 2 L + /3(9m 2 + 4s)) 
288ttto 2 s/3 2 



(239) 



(240) 
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V. Fermion-Higgs 



o - fG) = a (ft A™ - /*) = „ (/<M> - /G±) = ^gg^ 

48tts/3 2 



(241) 



( r(4 1 ^«^tG+) = 2887rm 2 s2/3 2 c 2 (l2fl 2 4^m 2 y? + m 2 £[12c 2 s 2 (m 2 - s)y?g 2 + 

+9<4(m 2 - s)y 2 + 4s ff 4 4] + */?[4« 5 4 *£ - 9c4m 2 y 4 ]) (243) 



_ g 4 (T^c w - 2s w t w Y f ) 2 {m 2 L + s/3) 



1287rm 2 s/? 2 c 2 



(244) 



a (/flffW - /_G+) = a (/IV) - /+G-) = a (/fUW - /_G+) = a (/i 1 ^ - / + G") = 
^ ° - / + *) = * (/?> ^ ^ /-*) = (-5) 

where / is any lepton or quark except and 6^, and /+ (/_) denotes isospin +1/2 (isospin —1/2) 
fermions. 

rr (+W itW +r+\ - ff 4 (c^ - 2s w t w Y f ) 2 {m 2 L + s(3) yt(m 2 L + sf3) 

a{t L H + ^tG )- 1287T m 2 S (3 2 c 2 w + 327rs2/3 2 ( 246 ) 

, - tW*) = 2a (ftW^W - G&) = ^ (247) 

/, (i) „(i) , r _\ .9 4 (c„, - 2g 1fl t, u y / ) 2 (m 2 I ) + S P) , VtL[sg 2 {c 2 w - 2s 2 y 6 ) + cg> - m 2 )y t 2 - cjs(Py?\ 
a \ L ~ ^ ) 128nm 2 s(3 2 c 2 w 32ns 2 (3 2 c 2 w 

(248) 



VI. Higgs-Higgs 

a G+G-") = a (V^ 1 ) G+G^ = 8m2 V £ + A( 2a + m V + 4A g 2r7t2 + 

V /V J 64nm 2 sf3 2 

(249) 

a (hWhW -^HH)=a (a^ A™ -+ Gg) = (250) 
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( H W H (i) ^ GG \ =a (a^A^ - JJJjl = 8m2X 3 2c l L + gKgg + m V + 4Ag 2 c> 2 + 4A 2 m 2 cj] 
V J \ J 647rm 2 s/3 2 c 4 

(251) 

„ k>^> - zz) = „ (a» A* - zz) = g4Ms+ tit;:i (s ' 2m2)£) 



(253) 



, (, 4 <„^> _ tf ) . griHi+g^g) (2M) 



a (Al"A<« -. H>» - If) = c (H»'H<y - ff<« - if) = 1 rf i < 255 > 

V / V + / o-'mfl s-m,,™V+ml,,,rt r ™ 



g 2 m(3 (s - m 2 H(2) ) 2 + m 2 H(2) T 2 H(2) 



( „(D u (i) n+n+ \ -lQg 2 m 2 \cjL + (3 ((2s + m 2 )g 4 - 8Xm 2 g 2 c 2 w + 16A 2 m 2 c 4 ) 
H+ ^° ° J' 1287rm 2 S /? 2 c 4 , 



(256) 



(ttWttW +F\ VciyH-Sscjy 2 - Am 2 g 2 s 2 w )L - 432^4 yf + sp 3 (20sj - I2sj + 9)g 4 - lU S p C 2 w y 2 g 2 s 2 w 

°{ H + H - ^ U )= il52^« 

(257) 



a (h^H™ bb) = 1152 J, W (72c^ 2 (-3 S c^ 2 + . 9 2 m 2 (44 - 3)L) + 



11527rs 2 /3 2 c: 

-432s/3<4</ t 4 + s/3 3 (204 - 24s 2 , + 9)g 4 - 36 S /3c 2 „y 2 5 2 (4 - 7s 2 w + 3) 



(258) 

/ „(i)„(i) ~+~_\ 6m 2 g 2 S (g 2 + 4Ac|)L - 48A 2 m 2 St j + g (g 4 (m 4 - 7m 2 s - 3s 2 ) - 12Am 2 sc|g 2 ) 
" ^ G G J - 192^WZ 



192ttto 2 s 2 /? 2 c 4 

(259) 

a -^HH)=<t (h^H™ GG) = la - G+G") = 

_ 8m 2 g 2 XL + (3 ((2s + m 2 )g 4 + A\m 2 g 2 + 4A 2 m 2 ) 
~~ 1287rm 2 s/3 2 

(rrWrrW r rr\ £_ ^m 2 c 2 w sL ~ (3[4(1 - 2s 2 w ) 2 m 4 + 5 (92s 4 - 140^ + 47) - 24 5 2 C 4 ] ) 
o \H + H_ GH) (260) 

„(hV H V^ 77 \ .9 4 (1 - 24) 4 (4m 2 (, - 2m 2 )L + s£a + 4m 2 )) 
f „(D„(i) 7 A .9 4 4(1 - 24) 2 (4m 2 ( S - 2m 2 )£ + ,s/3(.s + 4m 2 )) 
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a (h^H^ - 77 ) = i > >1 (263) 



V / 247TS 2 /? 2 



(264) 



m m - ^ 4 /3(0/(^(i-24)-24)-^^) 



l + 



/ (!),(!) ^ r \ = 6m 2 g 2 s(g 2 - 2A4)L - g ((m 4 - 7sm 2 - 3s 2 )g 4 + 6Am 2 s44g 4 - 12A 2 m 2 «4) 
CT l ^ J 192vrm 2 s 2 /3 2 c 4 , 



1927TTO 2 S 2 /3 2 4 

(266) 



^(i) A (i) ^ G+G -\ = ff 4 (24m 2 4sL + 0[ - 4(1 - 24) 2 m 4 + S m 2 (-92 S 4 + 1404 - 47) + 24s 2 c 4 ]) 



1927TTO 2 S 2 /3 2 C 4 

(267) 



a (HWAW - W+W-) = .9 4 (l2m 2 (2m 2 + 3 )£ + ^(32m 2 + .)) 

V / 967TS 3 /? 2 



96tTS 3 /? 2 

/ (1) (1) x -54y 2 (g 2 m 2 + 4( S - 2m 2 );/ 2 ) £ + /?[- 54s?/|c| - 27g 2 y 2 S + g 4 S 3 2 (32 S j - 244 + g)/4<j 

CT r A "J = 288^4 

(269) 



(268) 



, . N c 8g 4 (y 2 + Y 2 ) 

\ J 967TSC 4 , 



(270) 



a (h^H™ -> = a (i^fff - GG+) = 

6m 2 g 2 s(g 2 - 2A)L - 12A 2 m 2 s - B[{m A ~ lam 2 ~ 3s 2 )g 4 + 6Am 2 s4g 4 ] 



192ttto 2 s 2 /3 2 

(271) 



■A 



, (h<«h™ GG+ ) . „ (Af - * G+ ) = x 

x (l2m 2 (l - 2*2,) (2 - 3«2,)L - /3[(4m 4 + 47sm 2 - 24s 2 )4 - 24(m 2 - s)ss 2 w c 2 w - 3s(m 2 + 4s)s 4 ] ) 

(272) 

a ( H^H™ - = a f A^ff^ -> = x (273) 

V + J V + / 96tts 3 /3 2 c 4 v 7 

x (6m 2 L[(l - 24)(4m 2 - 2ss 2 w + s) + s] + s/3[4s 4 w (s + Urn 2 ) + (1 - 2s 2 w )(s + 32m 2 )]) (274) 



(275) 



f ff (iWi) +7 A -6y 2 (2g 2 m 2 + sy 2 )L + /?[ S /? 2 g 4 - 6sg 2 4g 2 - 12 S y 4 ] 
a r ^+ ^* 6 J= 6tof (276) 

a {HUH? - //') = a {AV>H? - //') = §g (277) 
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B.3 Production Cross Sections for strongly interacting KK modes at Hadron 
Colliders 

We list below the spin-averaged amplitudes squared, summed over initial for the production of strongly 
interacting KK particles at hadron colliders [242 . The symbols and qn generically indicate up and 
down-type quarks of any generation. 



M(qq -> q ( >q 



4 5 | 



2m 2 (t - to 2 ) 2 + (u - m 2 ) 2 



(278) 



M(qq-*qWqW) 



9i 



2m 1 



1 



+ - 



23 



s (t — to 2 ) 2 t — in 2 
2s 2 8s 6(t-m 2 ) 



8{t-m 2 ) 



2\2T 



6 (t-m 2 ) 2 3(t-m 2 ) 



(279) 



M(qq^ q {1) q (1) ) 



27 



m z 6 



t — to 



6- 



(u — m 2 ) 2 (t — m 2 ) 2 (t — m 2 )(u — m 2 ) 

(t — m 2 ) 2 (u — m 2 ) 2 -s 2 
+2 3^ + 3^ Jrr + 4 



(u — m 2 ) 2 (t — m 2 ) 2 (t — m 2 )(u — to 2 ) 



- 5 



(280) 



M(gg^q^q^) 



9i 



(t — m 2 )(u — m 2 ) \6(t — m 2 )(u — to 2 ) 8 



6(t — m 2 )(u — to 2 ) 8 



2\1 



17 3(t-m 2 )(u — m 



6(t - m 2 ){u - m 2 ) 24 



4s 2 



(281) 



5s 2 



12(t — m 2 ) 2 (t — to 2 ) 2 '{u — m 2 ) 

lls(w — to 2 ) 5(u — to 2 ) 2 (u — m 2 ) 3 



6(t-m 2 ) 2 \2(t-m 2 ) 2 s(t - to 2 ) 2 



(282) 



M(qq' -> 9 (1 V (1) ) 



18 



4m 



(t - TO 2 ) 2 



5 + 4 



(t — TO 2 ) 2 t — TO 2 



(283) 



Miqq'^qVq'W] 



9 



1 



(t-TO 2 ) 2 ' 4 ' (t-TO 2 ) 2 



(284) 



9t 



2s 3 



4s 



(f — m 2 ) 2 (u — m 2 ) 2 (t — m 2 )(u — to 2 ) 



-8- 



(t — to 2 ) 2 (u — to 2 ) 2 (t — m 2 )(u — m 2 ) 



(285) 



2 



M( qq ' - flW^gW^') ' - MM - g g> ff «' 

1 M — to 2 \ 5 4(u — to 2 ) 2(u — to 2 ) 2 



Si 



9 



2to 2 



< — TO 2 (t — TO 2 ) 5 



t — TO 2 (t — TO 2 ) 2 



(286) 
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M(gg^g^g^) 



9$£ 



3m 4 



+ 1 + 



+ (t - TO 2 ) 2 + (U - TO 2 ) 2 2 

(t — m 2 ) 2 (u — m 2 ) 2 
(s 2 + (t- m 2 ) 2 + (u - m 2 ) 2 ) 3 
As 2 {t — m 2 ) 2 (u — m 2 ) 2 



2 + (t - to 2 ) 2 + (u-rt 
s(t — m 2 )(u — in 2 ) 
(t - m 2 )(u - to 2 )" 



- + 



(287) 



M(qq 



9 W 9 



CD) 



M 

27 



4s 3 



57s 



108 s 



+ 



— m 2 ) 2 (u 
20s 2 



l ) 2 (t-m 2 )(u-m 2 ) 



(t — m 2 )(u — to 2 ) 



-93 + 



I08(t-m 2 ) 2 {u-m 2 ) 



(288) 
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